Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,313,472 G→A 6.6% D206D (GAC→GAT rhbC ← possible Iuc family aerobactin synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,313,4720GA 6.6% 230.1 / 4.2 76D206D (GAC→GATrhbCpossible Iuc family aerobactin synthesis protein
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base A (3/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GGATGTACTAAAATCGCTCTATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCCTCTGGTACCAAACTCTCAGA  >  CP000730/2313332‑2313612
                                                                                                                                            |                                                                                                                                            
ggATGTACTAAAATCGCTCTATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATgg                                                                                                                                              <  2:916460/141‑1 (MQ=255)
 gATGTACTAAAATCGCTCTATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGt                                                                                                                                             <  2:912536/141‑1 (MQ=255)
   tGTACTAAAATCGCTCTATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCt                                                                                                                                           <  1:828129/141‑1 (MQ=255)
     tACTAAAATCGCTCTATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCttt                                                                                                                                         >  1:418339/1‑141 (MQ=255)
      aCTAAAATCGCTCTATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCtttt                                                                                                                                        <  2:515951/141‑1 (MQ=255)
      aCTAAAATCGCTCTATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCtttt                                                                                                                                        >  1:495690/1‑141 (MQ=255)
         aaaaTCGCTCTATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCa                                                                                                                                     >  2:358412/1‑141 (MQ=255)
              cGCTCTATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAatcat                                                                                                                                >  1:338065/1‑141 (MQ=255)
               gCTCTATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAatcatc                                                                                                                               >  2:377728/1‑141 (MQ=255)
                ctctATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAAtcatca                                                                                                                              >  1:185193/1‑141 (MQ=255)
                 tctATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATcatcat                                                                                                                             >  2:717479/1‑141 (MQ=255)
                     aaGTCTTCACTTTTCACACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAAtt                                                                                                                         <  2:110665/140‑1 (MQ=255)
                          ttCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCaa                                                                                                                    >  1:354252/1‑141 (MQ=255)
                              cTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGATCTTTTTCAATCATCATAATTTGCAAAggg                                                                                                                >  1:348870/1‑141 (MQ=255)
                                 ttCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGAtt                                                                                                             <  1:969236/141‑1 (MQ=255)
                                       ccTAAACTCTTTAAAAACACACGCATCTGGTTGAAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCtttt                                                                                                       >  1:198357/1‑141 (MQ=255)
                                       ccTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCtttt                                                                                                       <  1:1015177/141‑1 (MQ=255)
                                        cTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCttttt                                                                                                      >  2:148504/1‑141 (MQ=255)
                                         tAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTc                                                                                                     <  2:478410/141‑1 (MQ=255)
                                         tAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTc                                                                                                     <  1:228218/141‑1 (MQ=255)
                                            aCTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCaaa                                                                                                  <  2:757066/141‑1 (MQ=255)
                                            aCTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCaaa                                                                                                  <  2:770665/141‑1 (MQ=255)
                                             ctctTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAAc                                                                                                 >  1:1078861/1‑141 (MQ=255)
                                               ctTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTc                                                                                               >  1:1080350/1‑141 (MQ=255)
                                               ctTTAAAAACACACGAATCTGATAGTAGTAATCGAGCATTAATTCATGCATAATTAATTGATCAGTACCAGCCATAGCTGTGCACACAGCATGGAATTCTTCAATCATCATAAGTCGCTAAGGGATTTCGTTTTCaccctc                                                                                               >  2:979520/1‑136 (MQ=255)
                                                tttAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCa                                                                                              >  2:93040/1‑141 (MQ=255)
                                                   aaaaaCANACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGt                                                                                           <  2:814305/141‑1 (MQ=255)
                                                   aaaaaCACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGt                                                                                           >  2:602513/1‑141 (MQ=255)
                                                   aaaaaCACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGt                                                                                           >  1:301454/1‑141 (MQ=255)
                                                    aaaaCACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTg                                                                                          >  2:34928/1‑141 (MQ=255)
                                                         acacGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCatat                                                                                     >  2:82702/1‑141 (MQ=255)
                                                         acacGAATCTGAGTGTAGTATTCGGGAATTATTTCATCAAGAATAAATTGAGCATGACCAGCCAGAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGTGACGTCTTTTGCAAACTCAGGTACatat                                                                                     >  2:102368/1‑141 (MQ=255)
                                                                  tGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGATCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTAct                                                                            <  1:937643/141‑1 (MQ=255)
                                                                    aTTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGATCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTActtc                                                                          >  2:330112/1‑141 (MQ=255)
                                                                     ttGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTActtct                                                                         <  1:1075941/141‑1 (MQ=255)
                                                                      tGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTActtctt                                                                        <  1:477021/141‑1 (MQ=255)
                                                                         agtaTTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCa                                                                     <  1:3423/141‑1 (MQ=255)
                                                                           taTTCGGGAATTATTTAATCGATAATAAATTGCTCATGACCATATATAGCTGTGGATACAACATGGTCTTTTTCAAATATCATACTCTGCAAAGGGATATCTTTTTTTATCTCAGGTTCATTTGCTCGTACTTCTTCCTTa                                                                   >  2:808745/1‑141 (MQ=255)
                                                                             ttCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGt                                                                 >  1:726501/1‑141 (MQ=255)
                                                                                 ggAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAg                                                             >  2:505665/1‑141 (MQ=255)
                                                                                 ggAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTATAACTCAGTTGCATATGCTGTTACTTCTTCCATAGTTAAg                                                             >  2:560826/1‑141 (MQ=255)
                                                                                  gAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAgg                                                            >  2:635696/1‑141 (MQ=255)
                                                                                     ttattTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCa                                                         <  2:161105/141‑1 (MQ=255)
                                                                                       attTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAAt                                                       <  1:764178/141‑1 (MQ=255)
                                                                                       attTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAAt                                                       <  2:303260/141‑1 (MQ=255)
                                                                                       attTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAAt                                                       <  1:1056269/141‑1 (MQ=255)
                                                                                        tttCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGATCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAAtt                                                      <  2:839513/141‑1 (MQ=255)
                                                                                         ttCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAAttt                                                     <  2:1067369/141‑1 (MQ=255)
                                                                                          tcatcaATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTc                                                    >  2:33792/1‑141 (MQ=255)
                                                                                           catcaATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCg                                                   >  2:242508/1‑141 (MQ=255)
                                                                                              caATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGttt                                                <  2:34420/141‑1 (MQ=255)
                                                                                              caATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGttt                                                <  2:303414/141‑1 (MQ=255)
                                                                                                 taataaATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTgg                                             <  2:1073612/141‑1 (MQ=255)
                                                                                                      aaTTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCNTAGTTAAGGGCANTTTCGTTTTGGTTaaa                                        <  1:1082699/141‑1 (MQ=255)
                                                                                                      aaTTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTaaa                                        <  1:468719/141‑1 (MQ=255)
                                                                                                      aaTTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTaaa                                        <  2:232165/141‑1 (MQ=255)
                                                                                                          gATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGgt                                    >  1:382201/1‑141 (MQ=255)
                                                                                                           aTCATTACCATCCATAGCTGTGCACACAACATGGTCTTGTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGgtg                                   <  1:314088/141‑1 (MQ=255)
                                                                                                           aTCATTACCATCCATAGCGGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGgtg                                   >  2:1024343/1‑141 (MQ=255)
                                                                                                            tCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCAGCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGgtgt                                  <  2:993985/141‑1 (MQ=255)
                                                                                                             cATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGgtgtg                                 <  2:990584/141‑1 (MQ=255)
                                                                                                                    atccatAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATg                          <  2:924994/141‑1 (MQ=255)
                                                                                                                     tccatAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGc                         >  2:22969/1‑141 (MQ=255)
                                                                                                                      ccatAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGcc                        <  2:1064569/141‑1 (MQ=255)
                                                                                                                      ccatAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGcc                        >  1:268829/1‑141 (MQ=255)
                                                                                                                           gCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCATACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCctct                   <  1:214574/141‑1 (MQ=255)
                                                                                                                           gCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCctct                   >  1:442438/1‑141 (MQ=255)
                                                                                                                           gCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCctct                   <  1:765123/141‑1 (MQ=255)
                                                                                                                           gCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCctct                   >  2:489617/1‑141 (MQ=255)
                                                                                                                           gCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCctct                   >  1:44799/1‑141 (MQ=255)
                                                                                                                            cTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTTTTTTTTTTTTTCATATGCTCGTACTTCTTCCATAGTTAAGTGCAATTTGGTTTTGGTTATATGGTGTGTTGGATGCCCCTCTg                  >  2:444766/1‑141 (MQ=255)
                                                                                                                             tgtgCACACAACATGATCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCCTCTgg                 >  1:777942/1‑141 (MQ=255)
                                                                                                                               tgCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCCTCTGGTa               <  2:1065362/141‑1 (MQ=255)
                                                                                                                               tgCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCACTCTGGTa               >  2:77272/1‑141 (MQ=255)
                                                                                                                                gCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCCTCTGGTAc              <  2:318574/141‑1 (MQ=255)
                                                                                                                                 cacacaACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCCTCTGGTAcc             <  2:738138/141‑1 (MQ=255)
                                                                                                                                       acaTGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCCTCTGGTACCAAActc       <  1:1055351/141‑1 (MQ=255)
                                                                                                                                            gTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCCTCTGGTACCAAACTCTCAGa  >  2:327110/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
GGATGTACTAAAATCGCTCTATAGTCTTCACTTTTCAAACCTAAACTCTTTAAAAACACACGAATCTGATTGTAGTATTCGGGAATTATTTCATCAATAATAAATTGATCATTACCATCCATAGCTGTGCACACAACATGGTCTTTTTCAATCATCATAATTTGCAAAGGGATTTCTTTTTCAAACTCAGGTGCATATGCTCGTACTTCTTCCATAGTTAAGGGCAATTTCGTTTTGGTTAAAGGGTGTGTTGGATGCCCCTCTGGTACCAAACTCTCAGA  >  CP000730/2313332‑2313612

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: