Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 580,366 G→T 94.6% R242S (CGC→AGC)  ompT ← DLP12 prophage; outer membrane protease VII (outer membrane protein 3b)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273580,3660GT94.6% 730.3 / 29.2 249R242S (CGC→AGC) ompTDLP12 prophage; outer membrane protease VII (outer membrane protein 3b)
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base T (110/125);  total (117/132)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

TTATTGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTG  >  CP009273/580221‑580510
                                                                                                                                                 |                                                                                                                                                
ttattGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGAt                                                                                                                                                <  2:184620/148‑1 (MQ=255)
 tattGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCtata                                                                                                                                               >  2:2109102/1‑148 (MQ=255)
 tattGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCtata                                                                                                                                               >  2:165541/1‑148 (MQ=255)
 tattGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATa                                                                                                                                               >  2:2553990/1‑148 (MQ=255)
    tgtgATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGt                                                                                                                                            >  2:1705273/1‑148 (MQ=255)
          cATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCt                                                                                                                                      <  1:710943/148‑1 (MQ=255)
           aTAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCtt                                                                                                                                     <  2:1678577/148‑1 (MQ=255)
            tAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCttt                                                                                                                                    >  2:1120114/1‑148 (MQ=255)
             aaaGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCtttt                                                                                                                                   >  1:1662095/1‑148 (MQ=255)
              aaGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCttttt                                                                                                                                  >  1:385501/1‑148 (MQ=255)
               aGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTc                                                                                                                                 >  2:1430354/1‑148 (MQ=255)
               aGTGAAGTATTACCATTTATATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGATAGTTGTGAAGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATCAGTGATTCTTTTTc                                                                                                                                 >  2:3561247/1‑148 (MQ=255)
                gTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTcc                                                                                                                                >  1:671723/1‑148 (MQ=255)
                gTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTcc                                                                                                                                >  1:2550715/1‑148 (MQ=255)
                  gAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCg                                                                                                                              >  1:3010402/1‑148 (MQ=255)
                  gAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCg                                                                                                                              <  1:2873372/148‑1 (MQ=255)
                             ttttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgtgt                                                                                                                   <  2:2981679/148‑1 (MQ=255)
                               ttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTc                                                                                                                 <  2:2296802/148‑1 (MQ=255)
                                tttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCa                                                                                                                <  2:2953104/148‑1 (MQ=255)
                                 ttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCAt                                                                                                               <  2:3550630/148‑1 (MQ=255)
                                 ttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCAt                                                                                                               >  2:3083866/1‑148 (MQ=255)
                                 ttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCAt                                                                                                               <  2:1441990/148‑1 (MQ=255)
                                 ttATTCGTAACCCGAGTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCAt                                                                                                               <  1:2786171/148‑1 (MQ=255)
                                    ttCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGt                                                                                                            >  1:2566882/1‑148 (MQ=255)
                                     tCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGtt                                                                                                           <  1:3087406/148‑1 (MQ=255)
                                     tCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGtt                                                                                                           >  1:1006448/1‑148 (MQ=255)
                                     tCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGtt                                                                                                           >  2:3322155/1‑148 (MQ=255)
                                        tAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGCAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATc                                                                                                        >  1:584352/1‑148 (MQ=255)
                                        tAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGCAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATc                                                                                                        >  1:3133949/1‑148 (MQ=255)
                                          aCCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAg                                                                                                      <  1:934573/148‑1 (MQ=255)
                                          aCCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAg                                                                                                      >  2:2072703/1‑148 (MQ=255)
                                             cGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatg                                                                                                   <  1:2332153/148‑1 (MQ=255)
                                             cGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatg                                                                                                   >  2:1751276/1‑148 (MQ=255)
                                              gATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatga                                                                                                  >  1:2921302/1‑148 (MQ=255)
                                               aTTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatgat                                                                                                 >  2:3453140/1‑148 (MQ=255)
                                                ttCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                                <  1:3277111/148‑1 (MQ=255)
                                                       cgcCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGGTTCCACCCa                                                                                         <  1:1798985/148‑1 (MQ=255)
                                                       cgcCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCa                                                                                         <  1:1324482/148‑1 (MQ=255)
                                                            tCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGc                                                                                    >  1:2439869/1‑148 (MQ=255)
                                                              aaCATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTg                                                                                  <  1:1910510/148‑1 (MQ=255)
                                                                cATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTa                                                                                >  1:2803903/1‑148 (MQ=255)
                                                                cATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTa                                                                                >  1:1893817/1‑148 (MQ=255)
                                                                 aTAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAt                                                                               <  2:3135530/148‑1 (MQ=255)
                                                                  tAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAtt                                                                              >  1:2213533/1‑148 (MQ=255)
                                                                    aaCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTa                                                                            >  1:944380/1‑148 (MQ=255)
                                                                     aCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaa                                                                           >  2:474713/1‑148 (MQ=255)
                                                                     aCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaa                                                                           >  2:1948838/1‑148 (MQ=255)
                                                                     aCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaa                                                                           >  2:862086/1‑148 (MQ=255)
                                                                     aCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaa                                                                           <  1:139054/148‑1 (MQ=255)
                                                                      cTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaaa                                                                          <  1:3410807/148‑1 (MQ=255)
                                                                       ttttGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAt                                                                         <  1:375223/148‑1 (MQ=255)
                                                                        tttGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAtg                                                                        <  2:318599/148‑1 (MQ=255)
                                                                        tttGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCGTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAtg                                                                        <  2:3134014/148‑1 (MQ=255)
                                                                           gCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGc                                                                     >  1:1105782/1‑148 (MQ=255)
                                                                            cGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGcc                                                                    <  1:2109102/148‑1 (MQ=255)
                                                                            cGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGcc                                                                    <  2:1731740/148‑1 (MQ=255)
                                                                              ttAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCAc                                                                  >  1:2817421/1‑148 (MQ=255)
                                                                              ttAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCAc                                                                  >  1:2817991/1‑148 (MQ=255)
                                                                                aGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCg                                                                <  1:137705/148‑1 (MQ=255)
                                                                                  gtgtGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAg                                                              >  1:2282483/1‑148 (MQ=255)
                                                                                  gtgtGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAg                                                              <  2:565497/148‑1 (MQ=255)
                                                                                  gtgtGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAg                                                              <  1:1313396/148‑1 (MQ=255)
                                                                                      gACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCa                                                          >  2:426222/1‑148 (MQ=255)
                                                                                       aCGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCaa                                                         >  1:1725235/1‑148 (MQ=255)
                                                                                          taataaCCTTCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAAt                                                      <  2:2240699/148‑1 (MQ=255)
                                                                                          taataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAAt                                                      <  1:11683/148‑1 (MQ=255)
                                                                                          taataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAAt                                                      <  2:2722547/148‑1 (MQ=255)
                                                                                              aaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTc                                                  <  1:2256814/148‑1 (MQ=255)
                                                                                              aaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTc                                                  <  2:1533896/148‑1 (MQ=255)
                                                                                                 cTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATa                                               <  1:462877/148‑1 (MQ=255)
                                                                                                 cTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATa                                               <  2:2394326/148‑1 (MQ=255)
                                                                                                 cTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATa                                               <  2:3366753/148‑1 (MQ=255)
                                                                                                  tGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTtaataa                                              >  1:407386/1‑148 (MQ=255)
                                                                                                  tGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATaa                                              <  2:2648876/148‑1 (MQ=255)
                                                                                                   gCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAAc                                             >  2:2166219/1‑148 (MQ=255)
                                                                                                    cATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACg                                            <  2:2948398/148‑1 (MQ=255)
                                                                                                     aTTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGa                                           <  1:1992546/148‑1 (MQ=255)
                                                                                                       tGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATa                                         <  2:2527773/148‑1 (MQ=255)
                                                                                                            gCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTc                                    >  1:3129893/1‑148 (MQ=255)
                                                                                                            gCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTc                                    >  2:1150618/1‑148 (MQ=255)
                                                                                                                cAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGt                                <  2:971580/148‑1 (MQ=255)
                                                                                                                 aGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGGTCAAAATCTTAATAACGATAACTTCCAGTc                               >  1:1548899/1‑148 (MQ=255)
                                                                                                                  gAATAGTAATTTTGGTCTTTGACCTTACTGCTATTAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCTAAATCTTCATAACGATAACTTCCAGTCa                              <  2:593323/148‑1 (MQ=255)
                                                                                                                   aaTAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCaa                             >  2:2885274/1‑148 (MQ=255)
                                                                                                                     tAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTTCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGc                           >  1:3451583/1‑148 (MQ=255)
                                                                                                                     tAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGGATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTACAGTCAAga                           >  1:1530860/1‑147 (MQ=255)
                                                                                                                      aGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGcc                          >  2:3264160/1‑148 (MQ=255)
                                                                                                                      aGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGcc                          >  2:799798/1‑148 (MQ=255)
                                                                                                                       gTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCa                         <  1:564527/148‑1 (MQ=255)
                                                                                                                         aaTTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAAt                       <  1:3295447/148‑1 (MQ=255)
                                                                                                                           ttttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCAGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGt                     <  2:2921302/148‑1 (MQ=255)
                                                                                                                            tttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTa                    <  1:1672510/148‑1 (MQ=255)
                                                                                                                             ttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGGCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAg                   <  2:2198437/148‑1 (MQ=255)
                                                                                                                              tGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAgg                  >  1:2013410/1‑148 (MQ=255)
                                                                                                                                   tttGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAtt             <  2:679980/148‑1 (MQ=255)
                                                                                                                                   tttGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAtt             <  2:695828/148‑1 (MQ=255)
                                                                                                                                         cTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCCATGTAGGGCATTTTaaaa       >  2:2946699/1‑148 (MQ=255)
                                                                                                                                         cTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTaaaa       <  2:1252694/148‑1 (MQ=255)
                                                                                                                                             cTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACgtt   <  1:1252795/148‑1 (MQ=255)
                                                                                                                                             cTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACgtt   <  2:2282483/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTATTGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTG  >  CP009273/580221‑580510

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