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.8% 998.6 / ‑6.3 305R242S (CGC→AGC) ompTDLP12 prophage; outer membrane protease VII (outer membrane protein 3b)
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (171/133);  total (171/134)
Fisher's exact test for biased strand distribution p-value = 4.39e-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 310 total aligned reads displayed.

TATTATTGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTT  >  CP009273/580219‑580512
                                                                                                                                                   |                                                                                                                                                  
  ttattGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATaa                                                                                                                                                <  1:5098609/150‑1 (MQ=255)
  ttattGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATaa                                                                                                                                                >  1:2813250/1‑150 (MQ=255)
  ttattGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATaa                                                                                                                                                <  2:3846218/150‑1 (MQ=255)
   tattGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAg                                                                                                                                               >  2:3422674/1‑150 (MQ=255)
   tattGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAg                                                                                                                                               <  1:3687443/150‑1 (MQ=255)
    attGTGATCATAAAGTGAAGTATTACCTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTg                                                                                                                                             >  1:2403011/1‑150 (MQ=255)
         gATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTc                                                                                                                                         <  1:5147135/150‑1 (MQ=255)
         gATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTc                                                                                                                                         >  1:3861702/1‑150 (MQ=255)
          aTCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCt                                                                                                                                        >  1:345968/1‑150 (MQ=255)
           tCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCtt                                                                                                                                       >  1:3224033/1‑150 (MQ=255)
            cATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCttt                                                                                                                                      >  1:25381/1‑150 (MQ=255)
            cATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCttt                                                                                                                                      >  2:3493448/1‑150 (MQ=255)
            cATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCttt                                                                                                                                      >  2:3796990/1‑150 (MQ=255)
            cATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCttt                                                                                                                                      <  1:5219728/150‑1 (MQ=255)
             aTAAAGTGAAGTATTACCTTTTTTATTAGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCtttt                                                                                                                                     >  2:103781/1‑150 (MQ=255)
               aaaGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTc                                                                                                                                   <  1:4572016/150‑1 (MQ=255)
               aaaGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTc                                                                                                                                   <  1:101981/150‑1 (MQ=255)
               aaaGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGAGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTc                                                                                                                                   >  2:1448394/1‑150 (MQ=255)
                  gTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCg                                                                                                                                >  2:2505718/1‑150 (MQ=255)
                   tGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCgg                                                                                                                               >  2:2231202/1‑150 (MQ=255)
                    gAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCggg                                                                                                                              >  1:4399488/1‑150 (MQ=255)
                     aaGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGt                                                                                                                             >  1:2681742/1‑150 (MQ=255)
                           tACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCAGGGTCATAgt                                                                                                                       >  1:3017462/1‑150 (MQ=255)
                            aCCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgtg                                                                                                                      >  1:3407208/1‑150 (MQ=255)
                             ccTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAgtgt                                                                                                                     >  1:1303968/1‑150 (MQ=255)
                                 ttttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCAt                                                                                                                 <  1:86913/150‑1 (MQ=255)
                                   ttATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCg                                                                                                               <  1:4964816/150‑1 (MQ=255)
                                    tATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGt                                                                                                              >  1:170985/1‑150 (MQ=255)
                                     aTTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGtt                                                                                                             <  2:1164121/150‑1 (MQ=255)
                                           aaCCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAga                                                                                                       >  2:4210939/1‑150 (MQ=255)
                                             cccGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAgatg                                                                                                     >  2:1297423/1‑150 (MQ=255)
                                                gATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                                  <  1:2925218/150‑1 (MQ=255)
                                                gATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGAtt                                                                                                  <  1:4550281/150‑1 (MQ=255)
                                                 aTTCCATGCGCCTTCAACTTAATCTTTTGCGTTAGGTGTGACGTAATTACCTGCTTTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTc                                                                                                 <  1:1770064/150‑1 (MQ=255)
                                                  ttCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTcc                                                                                                >  1:1402172/1‑150 (MQ=255)
                                                  gtCCATGCGCCGTCAACATAAACTTGTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTcc                                                                                                <  2:2776056/149‑1 (MQ=255)
                                                    ccATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCAc                                                                                              >  2:1409373/1‑150 (MQ=255)
                                                    ccATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGGCATAGTGTTCATCGTTATCAGATGATTCCAc                                                                                              <  1:824911/150‑1 (MQ=255)
                                                     cATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCAcc                                                                                             >  2:4284302/1‑150 (MQ=255)
                                                      aTGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCAccc                                                                                            >  2:1119416/1‑150 (MQ=255)
                                                      aTGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCAccc                                                                                            >  2:4248346/1‑150 (MQ=255)
                                                       tGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCa                                                                                           <  1:3343199/150‑1 (MQ=255)
                                                       tGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCa                                                                                           <  1:2514568/150‑1 (MQ=255)
                                                       tGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCa                                                                                           <  1:549991/150‑1 (MQ=255)
                                                         cgcCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGc                                                                                         >  2:4085173/1‑150 (MQ=255)
                                                            cTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGc                                                                                      <  1:1211600/150‑1 (MQ=255)
                                                             ttCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCt                                                                                     >  2:983482/1‑150 (MQ=255)
                                                             ttCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCt                                                                                     <  2:3089885/150‑1 (MQ=255)
                                                                  cATATACTTTTGCGTTAGGTGTGACGTTATAAGCTGCATTGACTGCAACAGAATAGTTATTTTGGTCTTTGACCTTACTGCTATTAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAtt                                                                                <  1:1494191/150‑1 (MQ=255)
                                                                   aTAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAttt                                                                               >  2:5249800/1‑150 (MQ=255)
                                                                   aTAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTAttt                                                                               <  1:1945939/150‑1 (MQ=255)
                                                                     aaaCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTaa                                                                             <  2:286344/150‑1 (MQ=255)
                                                                       aCTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAt                                                                           <  2:4464093/150‑1 (MQ=255)
                                                                        cTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAtg                                                                          >  2:548849/1‑150 (MQ=255)
                                                                         ttttGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAAtgt                                                                         <  1:2579147/150‑1 (MQ=255)
                                                                             gCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCa                                                                     <  2:917542/150‑1 (MQ=255)
                                                                                  aGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGat                                                                >  1:927403/1‑149 (MQ=255)
                                                                                  aGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAg                                                                >  2:2360130/1‑150 (MQ=255)
                                                                                   ggTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGt                                                               >  1:3423416/1‑150 (MQ=255)
                                                                                     tgtgACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTc                                                             >  2:334506/1‑150 (MQ=255)
                                                                                       tgtcGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCAGCTGTATTTAAATGTGCCACCGAGTTCaa                                                           <  2:2943594/147‑1 (MQ=255)
                                                                                       tgACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCaa                                                           >  2:176062/1‑150 (MQ=255)
                                                                                         aCGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCaaaa                                                         >  2:4057646/1‑150 (MQ=255)
                                                                                           gTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAAta                                                       >  1:4723326/1‑149 (MQ=255)
                                                                                            taataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCt                                                      <  1:3562542/150‑1 (MQ=255)
                                                                                             aataaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCtt                                                     <  1:1172735/150‑1 (MQ=255)
                                                                                                aaCCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCAt                                                  <  2:5038320/150‑1 (MQ=255)
                                                                                                      cATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGAt                                            <  2:4778852/150‑1 (MQ=255)
                                                                                                      cATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGAt                                            <  2:2708264/150‑1 (MQ=255)
                                                                                                          gACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACt                                        <  2:1526873/150‑1 (MQ=255)
                                                                                                          gACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACt                                        <  1:995365/150‑1 (MQ=255)
                                                                                                          gACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACt                                        <  2:5174355/150‑1 (MQ=255)
                                                                                                          gACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGGCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACt                                        <  2:2178979/150‑1 (MQ=255)
                                                                                                           aCTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACtt                                       >  1:2219383/1‑150 (MQ=255)
                                                                                                            cTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTc                                      <  2:5226907/150‑1 (MQ=255)
                                                                                                            cTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTc                                      >  1:4229987/1‑150 (MQ=255)
                                                                                                            cTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTc                                      >  2:522335/1‑150 (MQ=255)
                                                                                                             tGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTcc                                     <  2:2668558/150‑1 (MQ=255)
                                                                                                             tGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTcc                                     <  1:480712/150‑1 (MQ=255)
                                                                                                               cAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAg                                   >  2:1640470/1‑150 (MQ=255)
                                                                                                                aaCAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGt                                  <  1:4996337/150‑1 (MQ=255)
                                                                                                                 aCAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATTTGCCACCGAGTTCAAAATATTCATAACGATAACATCCAGTc                                 >  1:4014582/1‑150 (MQ=255)
                                                                                                                  cAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCa                                <  1:2048051/150‑1 (MQ=255)
                                                                                                                  cAGAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCa                                >  1:1169307/1‑150 (MQ=255)
                                                                                                                    gAATAGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGCGGTCAAAATCTTCATAACGATAAATTACAGTCAAg                              >  1:4900282/1‑150 (MQ=255)
                                                                                                                        aGTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCGGTATTTAAATTTGCCACCGAGTTCAAAATCTTCATAACGATAAATTACAGTCAATCCaa                          >  1:4465856/1‑150 (MQ=255)
                                                                                                                         gTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAAt                         <  2:3224033/150‑1 (MQ=255)
                                                                                                                         gTAATTTTGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAAt                         <  2:1703225/150‑1 (MQ=255)
                                                                                                                               ttGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAggg                   <  1:216509/150‑1 (MQ=255)
                                                                                                                                tGGTCTTTGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGc                  <  1:1676481/150‑1 (MQ=255)
                                                                                                                                 ggTCTTTGACCTTACTGCTATAAGTTATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATTTGACACCGAGTTCAAAATCTTCATAACGATAAATTCCAGTCAAGCAAATGTAGGGCa                 >  2:3059611/1‑150 (MQ=255)
                                                                                                                                     tttGACCTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCAtttt             >  1:4183243/1‑150 (MQ=255)
                                                                                                                                          ccTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAAc        <  1:4076094/150‑1 (MQ=255)
                                                                                                                                           cTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACg       >  1:2195580/1‑150 (MQ=255)
                                                                                                                                           cTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACg       >  2:1941392/1‑150 (MQ=255)
                                                                                                                                           cTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACg       >  2:74486/1‑150 (MQ=255)
                                                                                                                                           cTTACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACg       >  1:662348/1‑150 (MQ=255)
                                                                                                                                             tACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACgtt     <  2:667595/150‑1 (MQ=255)
                                                                                                                                             tACTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACgtt     >  2:792272/1‑150 (MQ=255)
                                                                                                                                               cTGCTATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACgttgt   >  1:4323047/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TATTATTGTGATCATAAAGTGAAGTATTACCTTTTTTATTCGTAACCCGATTCCATGCGCCTTCAACATAAACTTTTGCGTTAGGTGTGACGTAATAACCTGCATTGACTGCAACAGAATAGTAATTTTGGTCTTTGACCTTACTGCGATAAGTGATTCTTTTTCCCGGGTCATAGTGTTCATCGTTATCAGATGATTCCACCCAGCCGCTGTATTTAAATGTGCCACCGAGTTCAAAATCTTCATAACGATAACTTCCAGTCAAGCCAATGTAGGGCATTTTAAAACGTTGTT  >  CP009273/580219‑580512

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