Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,597,517 T→G 15.7% L81* (TTA→TGA)  lsrC → autoinducer 2 import system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,597,5170TG15.7% 691.5 / 157.9 306L81* (TTA→TGA) lsrCautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base T (124/134);  new base G (20/28);  total (144/162)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 309 total aligned reads displayed.

TTTAAGTGTGCAAACGCTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAGGGTTGTACAGAGGCA  >  CP009273/1597368‑1597666
                                                                                                                                                     |                                                                                                                                                     
     gtgtGCAAACGCTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGa                                                                                                                                                  <  1:1771504/150‑1 (MQ=255)
        tgCAAACGCTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACg                                                                                                                                               <  2:1029905/150‑1 (MQ=255)
        tgCAAACGCTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACg                                                                                                                                               <  1:2328845/150‑1 (MQ=255)
        tgCAAACGCTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTGACTGAACg                                                                                                                                               <  2:3545019/150‑1 (MQ=255)
          cAAACGCTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCa                                                                                                                                             >  1:3741354/1‑150 (MQ=255)
           aaaCGCTGACCATGGTTGATAGCAGCGCGCAAATCCTGATCCTGCTGGCAAAGGGCGCGACGCTGGTAATGCTTACGCGCAAGATTGATGTGTCAGTGGGTTCGATTACCGGAATGTGCGCGGAGCTGTTGGGGATGATACTGAACGCAg                                                                                                                                            >  2:1002439/1‑150 (MQ=255)
             aCGCTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGa                                                                                                                                          <  2:1362266/150‑1 (MQ=255)
              cGCTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGat                                                                                                                                         >  1:3328447/1‑150 (MQ=255)
               gCTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGata                                                                                                                                        <  2:2310590/150‑1 (MQ=255)
               gCTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGata                                                                                                                                        <  1:4092701/150‑1 (MQ=255)
                cTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGatat                                                                                                                                       <  1:1052705/150‑1 (MQ=255)
                cTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGatat                                                                                                                                       >  1:1428854/1‑150 (MQ=255)
                cTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGatat                                                                                                                                       >  2:745028/1‑150 (MQ=255)
                 tGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGGACGCAGGATAtt                                                                                                                                      >  2:1810053/1‑150 (MQ=255)
                 tGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATAtt                                                                                                                                      >  2:2643733/1‑150 (MQ=255)
                 tGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATAtt                                                                                                                                      >  2:3315406/1‑150 (MQ=255)
                  gACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTc                                                                                                                                     >  1:2390174/1‑150 (MQ=255)
                    ccATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCAc                                                                                                                                   >  2:1069459/1‑150 (MQ=255)
                     cATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACt                                                                                                                                  >  2:2729448/1‑150 (MQ=255)
                      aTGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTa                                                                                                                                 >  2:20795/1‑150 (MQ=255)
                      aTGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTa                                                                                                                                 <  1:2747642/150‑1 (MQ=255)
                       tGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTAc                                                                                                                                >  1:3223380/1‑150 (MQ=255)
                          tttATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTg                                                                                                                             <  2:4236029/150‑1 (MQ=255)
                            tataGCAGCGCGCAACTCCTGATCCTGCTGGCAATGGGCGCGAGGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGtt                                                                                                                           <  2:4064600/150‑1 (MQ=255)
                            tataGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGtt                                                                                                                           <  1:858910/150‑1 (MQ=255)
                            tataGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGtt                                                                                                                           <  2:1970233/150‑1 (MQ=255)
                             ataGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTGACTGAACGCAGGATATTCACTACCTGTTg                                                                                                                          <  1:504052/150‑1 (MQ=255)
                                  agcGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATCTTCACTACCTGTTGCTtgt                                                                                                                     >  1:4343036/1‑150 (MQ=255)
                                  agcGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTtgt                                                                                                                     <  2:4372152/150‑1 (MQ=255)
                                    cgcgcAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTgtgt                                                                                                                   >  2:1701624/1‑150 (MQ=255)
                                     gcgcAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTc                                                                                                                  <  2:3937225/150‑1 (MQ=255)
                                      cgcAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTcg                                                                                                                 <  1:28371/150‑1 (MQ=255)
                                         aaaTCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTGACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGa                                                                                                              >  1:2670130/1‑150 (MQ=255)
                                           atcctgatcctgCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTGACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACt                                                                                                            >  1:727923/1‑150 (MQ=255)
                                              ctgatcctgCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTa                                                                                                         <  1:3480453/150‑1 (MQ=255)
                                                    ctgctgGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTGACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCtt                                                                                                   <  2:2749337/150‑1 (MQ=255)
                                                      gctgGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTgg                                                                                                 <  2:3059848/150‑1 (MQ=255)
                                                      gctgGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTgg                                                                                                 <  2:2504224/150‑1 (MQ=255)
                                                       ctgGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGt                                                                                                <  2:2914028/150‑1 (MQ=255)
                                                       ctgGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTGACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGt                                                                                                >  1:3750216/1‑150 (MQ=255)
                                                         ggCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGttt                                                                                              <  2:1340004/150‑1 (MQ=255)
                                                              tGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGTGCGGTGCTGTTGGGGATGTGACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTcg                                                                                         <  2:2901624/150‑1 (MQ=255)
                                                              tGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTcg                                                                                         >  1:2244482/1‑150 (MQ=255)
                                                               gggCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTcgc                                                                                        >  1:418413/1‑150 (MQ=255)
                                                                ggCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGACTGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTcgcg                                                                                       >  1:3539875/1‑150 (MQ=255)
                                                                    cgACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGAt                                                                                   >  2:2377175/1‑150 (MQ=255)
                                                                    cgACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGAt                                                                                   >  1:1358045/1‑150 (MQ=255)
                                                                    cgACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGAt                                                                                   >  1:2978390/1‑150 (MQ=255)
                                                                     gACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGAtt                                                                                  <  2:1298286/150‑1 (MQ=255)
                                                                     gACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCGTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGAtt                                                                                  >  2:3649684/1‑150 (MQ=255)
                                                                      aCGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGAttt                                                                                 <  1:1143411/150‑1 (MQ=255)
                                                                       cGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGAtttt                                                                                >  2:1953095/1‑150 (MQ=255)
                                                                        gCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGACTGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGAttttt                                                                               >  1:3597439/1‑150 (MQ=255)
                                                                          tGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCa                                                                             <  1:4067383/150‑1 (MQ=255)
                                                                                 gCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGgtgt                                                                      <  1:4037567/150‑1 (MQ=255)
                                                                                  cTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTc                                                                     >  2:1411106/1‑150 (MQ=255)
                                                                                   ttACGCGCAATATTGATGTTTCTGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACGTGTTGCTTGTGTCGCGACTTTTCTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTcc                                                                    <  1:2186920/150‑1 (MQ=255)
                                                                                   ttACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTcc                                                                    <  1:4375063/150‑1 (MQ=255)
                                                                                    tACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCt                                                                   <  1:1980044/150‑1 (MQ=255)
                                                                                       gcgcAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGt                                                                >  1:1826663/1‑150 (MQ=255)
                                                                                        cgcAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTc                                                               <  2:200158/150‑1 (MQ=255)
                                                                                        cgcAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTc                                                               >  2:511497/1‑150 (MQ=255)
                                                                                         gcAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTcg                                                              <  2:4353586/150‑1 (MQ=255)
                                                                                                tGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCt                                                       >  2:3237480/1‑150 (MQ=255)
                                                                                                tGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTGACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCt                                                       <  1:4595984/150‑1 (MQ=255)
                                                                                                       tCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATc                                                <  2:4135173/150‑1 (MQ=255)
                                                                                                         aGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATccc                                              >  1:465563/1‑150 (MQ=255)
                                                                                                          gTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCt                                             <  1:4612340/150‑1 (MQ=255)
                                                                                                          gTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCt                                             >  1:2702693/1‑150 (MQ=255)
                                                                                                             ggTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGcc                                          <  1:2895954/150‑1 (MQ=255)
                                                                                                               ttCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTGACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCAt                                        >  1:2557386/1‑150 (MQ=255)
                                                                                                                tCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCAtt                                       >  1:4166050/1‑150 (MQ=255)
                                                                                                                tCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTGACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCAtt                                       <  2:2788761/150‑1 (MQ=255)
                                                                                                                 cGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCAttg                                      >  2:3152901/1‑150 (MQ=255)
                                                                                                                  gATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCAttgt                                     >  2:4494758/1‑150 (MQ=255)
                                                                                                                   aTTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCAttgtt                                    <  2:3223490/150‑1 (MQ=255)
                                                                                                                      aCCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGcc                                 >  1:151083/1‑150 (MQ=255)
                                                                                                                       ccGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCa                                <  2:2871367/150‑1 (MQ=255)
                                                                                                                       ccGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTGTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCa                                <  2:2280675/150‑1 (MQ=255)
                                                                                                                         ggAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTGAAGATCCCTGCCATTGTTGCCAcc                              >  1:4523432/1‑150 (MQ=255)
                                                                                                                          gAATGTGCGCTGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCAccc                             <  1:3380057/150‑1 (MQ=255)
                                                                                                                           aaTGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCt                            <  1:3462131/150‑1 (MQ=255)
                                                                                                                             tgtgCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTg                          >  2:2838513/1‑150 (MQ=255)
                                                                                                                               tgCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGc                        <  2:4153722/150‑1 (MQ=255)
                                                                                                                                gcgcGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCa                       >  2:3179734/1‑150 (MQ=255)
                                                                                                                                gcgcGGTGCTGTTGGGGATGTGACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCa                       <  2:3750210/150‑1 (MQ=255)
                                                                                                                                    ggTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGtt                   <  1:2859018/150‑1 (MQ=255)
                                                                                                                                    ggTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGATAAAGATCCCTGCCATTGTTTCCACCCTTGGCACGtt                   >  2:4131417/1‑150 (MQ=255)
                                                                                                                                     gTGCTGTTGGGGATGTGACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTa                  >  2:2014700/1‑150 (MQ=255)
                                                                                                                                      tGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAg                 >  2:1120103/1‑150 (MQ=255)
                                                                                                                                       gCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAgg                <  2:1805699/150‑1 (MQ=255)
                                                                                                                                       gCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAgg                <  1:4145766/150‑1 (MQ=255)
                                                                                                                                       gCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAgg                <  2:626556/150‑1 (MQ=255)
                                                                                                                                        cTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAggg               >  1:78479/1‑150 (MQ=255)
                                                                                                                                          gTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAGGGtt             <  1:362107/150‑1 (MQ=255)
                                                                                                                                             ggggATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAGGGTTGTa          >  2:3797753/1‑150 (MQ=255)
                                                                                                                                              gggATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAGGGTTGTAc         <  1:20795/150‑1 (MQ=255)
                                                                                                                                               ggATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAGGGTTGTACa        <  2:3398936/150‑1 (MQ=255)
                                                                                                                                               ggATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAGGGTTGTACa        <  2:4518351/150‑1 (MQ=255)
                                                                                                                                                  tGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAGGGTTGTACagag     <  1:3837856/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTAAGTGTGCAAACGCTGACCATGGTTTATAGCAGCGCGCAAATCCTGATCCTGCTGGCAATGGGCGCGACGCTGGTAATGCTTACGCGCAATATTGATGTTTCAGTGGGTTCGATTACCGGAATGTGCGCGGTGCTGTTGGGGATGTTACTGAACGCAGGATATTCACTACCTGTTGCTTGTGTCGCGACTTTACTGCTTGGTTTGCTCGCGGGATTTTTCAACGGTGTCCTGGTCGCGTGGCTAAAGATCCCTGCCATTGTTGCCACCCTTGGCACGTTAGGGTTGTACAGAGGCA  >  CP009273/1597368‑1597666

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: