Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 580,366 G→T 100% 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,3660GT99.1% 768.7 / ‑1.9 234R242S (CGC→AGC) ompTDLP12 prophage; outer membrane protease VII (outer membrane protein 3b)
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (128/104);  total (128/106)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 235 total aligned reads displayed.

TTATTGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTT  >  CP009273/580221‑580509
                                                                                                                                                 |                                                                                                                                               
      tgATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGa                                                                                                                                         <  1:532706/148‑1 (MQ=255)
           aTAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCtt                                                                                                                                    <  1:827156/148‑1 (MQ=255)
            tAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCttt                                                                                                                                   >  1:480672/1‑148 (MQ=255)
              aaGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCttttt                                                                                                                                 <  1:2080592/148‑1 (MQ=255)
              aaGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCttttt                                                                                                                                 <  1:538411/148‑1 (MQ=255)
               aGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTc                                                                                                                                <  2:2269257/148‑1 (MQ=255)
                  gAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCg                                                                                                                             <  1:2860319/148‑1 (MQ=255)
                   aaGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGNCTTCAACATAAACTTTTGCGTTAGNTGNGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTNTGACCTTACTGCTATAAGTGATTCTTTTTCCCgg                                                                                                                            >  1:3401560/1‑148 (MQ=255)
                     gTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGt                                                                                                                          <  2:1473592/148‑1 (MQ=255)
                      tATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTc                                                                                                                         >  1:1465985/1‑148 (MQ=255)
                      tATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTc                                                                                                                         >  2:2378682/1‑148 (MQ=255)
                        ttACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCAt                                                                                                                       <  1:3195658/148‑1 (MQ=255)
                         tACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATa                                                                                                                      >  1:1765424/1‑148 (MQ=255)
                         tACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATa                                                                                                                      <  2:3186600/148‑1 (MQ=255)
                           ccTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgt                                                                                                                    >  1:3145157/1‑148 (MQ=255)
                           ccTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgt                                                                                                                    >  2:35379/1‑148 (MQ=255)
                             ttttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgtgt                                                                                                                  <  1:2133642/148‑1 (MQ=255)
                              tttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGtt                                                                                                                 <  2:1832299/148‑1 (MQ=255)
                              tttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGtt                                                                                                                 <  2:54800/148‑1 (MQ=255)
                               ttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTc                                                                                                                >  2:2509504/1‑148 (MQ=255)
                                tttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCa                                                                                                               <  1:1960369/148‑1 (MQ=255)
                                tttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCa                                                                                                               >  2:114221/1‑148 (MQ=255)
                                 ttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCAt                                                                                                              <  2:3248488/148‑1 (MQ=255)
                                 ttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCAt                                                                                                              <  2:1952071/148‑1 (MQ=255)
                                    ttCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGt                                                                                                           >  2:2238040/1‑148 (MQ=255)
                                         aaCCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCa                                                                                                      >  2:2955809/1‑148 (MQ=255)
                                           cccGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAga                                                                                                    <  1:223508/148‑1 (MQ=255)
                                           cccGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAga                                                                                                    >  2:2760115/1‑148 (MQ=255)
                                            ccGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgat                                                                                                   >  2:848771/1‑148 (MQ=255)
                                              gATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatga                                                                                                 >  1:435371/1‑148 (MQ=255)
                                                ttCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                               <  1:695185/148‑1 (MQ=255)
                                                ttCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                               >  2:1301880/1‑148 (MQ=255)
                                                 tCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACATAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGATCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTc                                                                                              >  1:2998233/1‑148 (MQ=255)
                                                    aTGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCAc                                                                                           >  1:1511618/1‑148 (MQ=255)
                                                     tGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCAcc                                                                                          <  1:855202/148‑1 (MQ=255)
                                                      gcgcCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCAccc                                                                                         >  1:64745/1‑148 (MQ=255)
                                                       cgcCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCa                                                                                        >  1:594016/1‑148 (MQ=255)
                                                        gcCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGGCATAGTGTTCATCGTTATCAGATGATTCCACCCAg                                                                                       <  1:1639185/148‑1 (MQ=255)
                                                          cTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGcc                                                                                     >  1:1177979/1‑148 (MQ=255)
                                                            tCAACATAAACTTTAGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCACGGGTCATAGTGTTCATCGTTATCAGATCATTCCACCCAGCCGc                                                                                   >  1:50162/1‑148 (MQ=255)
                                                               aCATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGt                                                                                <  1:1806221/148‑1 (MQ=255)
                                                                cATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTa                                                                               <  1:2527463/148‑1 (MQ=255)
                                                                 aTAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAt                                                                              <  1:952219/148‑1 (MQ=255)
                                                                   aaaCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAttt                                                                            >  1:2619829/1‑148 (MQ=255)
                                                                    aaCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTa                                                                           >  1:1775806/1‑148 (MQ=255)
                                                                     aCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaa                                                                          <  1:2244948/148‑1 (MQ=255)
                                                                     aCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaa                                                                          <  1:185617/148‑1 (MQ=255)
                                                                      cTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaaa                                                                         <  1:1360271/148‑1 (MQ=255)
                                                                      cTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaaa                                                                         <  2:2578314/148‑1 (MQ=255)
                                                                      cTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGGCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaaa                                                                         <  2:1326751/148‑1 (MQ=255)
                                                                       ttttGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAt                                                                        >  2:1181471/1‑148 (MQ=255)
                                                                              ttAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCAc                                                                 >  1:1289043/1‑148 (MQ=255)
                                                                               tAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCAcc                                                                >  2:2856640/1‑148 (MQ=255)
                                                                                aGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCg                                                               <  1:793837/148‑1 (MQ=255)
                                                                                 ggTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGa                                                              >  2:1779938/1‑148 (MQ=255)
                                                                                    gtgACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGtt                                                           <  2:3323480/148‑1 (MQ=255)
                                                                                      gACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCa                                                         >  1:2467822/1‑148 (MQ=255)
                                                                                          taataaCCTGCATTGACTGCATCAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAAt                                                     <  2:2775917/148‑1 (MQ=255)
                                                                                          taataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCCAGTTCAAAAt                                                     >  2:718150/1‑148 (MQ=255)
                                                                                           aattaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCAGAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATc                                                    <  2:912490/148‑1 (MQ=255)
                                                                                            ataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCt                                                   <  1:1918709/148‑1 (MQ=255)
                                                                                             taaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCaaaaaatt                                                  >  1:613752/1‑144 (MQ=255)
                                                                                              aaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTc                                                 <  1:1055685/148‑1 (MQ=255)
                                                                                                ccTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATATTCAt                                               >  1:1557772/1‑148 (MQ=255)
                                                                                                 cTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATa                                              >  2:143690/1‑148 (MQ=255)
                                                                                                   gCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAAc                                            >  1:807910/1‑148 (MQ=255)
                                                                                                    cATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTGATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACg                                           >  1:2561852/1‑148 (MQ=255)
                                                                                                     aTTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGa                                          <  2:2093200/148‑1 (MQ=255)
                                                                                                     aTTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGa                                          <  2:2816603/148‑1 (MQ=255)
                                                                                                       tGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATa                                        >  1:1899197/1‑148 (MQ=255)
                                                                                                       tGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATa                                        <  2:1108417/148‑1 (MQ=255)
                                                                                                       tGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATa                                        >  1:3282888/1‑148 (MQ=255)
                                                                                                        gACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATaa                                       <  1:1601108/148‑1 (MQ=255)
                                                                                                         aCTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAAc                                      >  2:1577869/1‑148 (MQ=255)
                                                                                                           tGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACtt                                    >  1:2168325/1‑148 (MQ=255)
                                                                                                            gCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACtta                                   >  1:1504548/1‑147 (MQ=255)
                                                                                                            gCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTc                                   <  2:2044766/148‑1 (MQ=255)
                                                                                                             cAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTcc                                  <  1:1429241/148‑1 (MQ=255)
                                                                                                               aCAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAg                                >  1:2279074/1‑148 (MQ=255)
                                                                                                                cAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGt                               >  2:491801/1‑148 (MQ=255)
                                                                                                                cAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGt                               <  2:1989336/148‑1 (MQ=255)
                                                                                                                      aGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGcc                         >  1:291189/1‑148 (MQ=255)
                                                                                                                       gTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCa                        <  2:2859855/148‑1 (MQ=255)
                                                                                                                       gTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCa                        <  2:645635/148‑1 (MQ=255)
                                                                                                                        tAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCaa                       >  2:1426293/1‑148 (MQ=255)
                                                                                                                         aaTTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAAt                      >  1:3504697/1‑148 (MQ=255)
                                                                                                                         aaTTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAAt                      >  2:1043336/1‑148 (MQ=255)
                                                                                                                          aTTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCaaag                     >  2:2672063/1‑146 (MQ=255)
                                                                                                                           ttttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTACAGTCAAGCCAATGt                    >  2:3111099/1‑148 (MQ=255)
                                                                                                                            tttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTa                   <  1:1998516/148‑1 (MQ=255)
                                                                                                                            tttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCCAGTTCAAAATATTCATAACGATAAATTACAGTCAAGCCCATGTa                   >  1:2670141/1‑148 (MQ=255)
                                                                                                                             ttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAg                  <  1:2633120/148‑1 (MQ=255)
                                                                                                                                    ttGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAttt           <  2:3423811/148‑1 (MQ=255)
                                                                                                                                       aCCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTaa        >  2:1525744/1‑148 (MQ=255)
                                                                                                                                       aCCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATATTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATGTTaa        >  2:607940/1‑148 (MQ=255)
                                                                                                                                        ccTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTaaa       <  2:3040871/148‑1 (MQ=255)
                                                                                                                                        ccTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTaaa       >  2:1953551/1‑148 (MQ=255)
                                                                                                                                          ttACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAAc     <  2:2624202/148‑1 (MQ=255)
                                                                                                                                           tACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACg    <  2:2015272/148‑1 (MQ=255)
                                                                                                                                             cTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACgtt  >  1:1082225/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TTATTGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTT  >  CP009273/580221‑580509

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: