Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,833,272 T→G 24.1% S28R (AGC→CGC)  nlpA ← cytoplasmic membrane lipoprotein‑28

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,833,2720TG24.1% 385.2 / 172.3 216S28R (AGC→CGC) nlpAcytoplasmic membrane lipoprotein‑28
Reads supporting (aligned to +/- strand):  ref base T (68/96);  new base G (24/28);  total (92/124)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 221 total aligned reads displayed.

TGTCGCATCGTTCGGTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAAC  >  CP009273/3833123‑3833421
                                                                                                                                                     |                                                                                                                                                     
tGTCGCATCGTTCGGTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCgctgct                                                                                                                                                       >  1:3403858/1‑150 (MQ=255)
      aTCGTTCGGTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTg                                                                                                                                                 <  2:4244962/150‑1 (MQ=255)
      aTCGTTCGGTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTg                                                                                                                                                 <  2:2306617/150‑1 (MQ=255)
       tCGTTCGGTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTgg                                                                                                                                                >  2:2228163/1‑150 (MQ=255)
        cGTTCGGTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGt                                                                                                                                               >  2:4531628/1‑150 (MQ=255)
        cGTTCGGTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGt                                                                                                                                               >  1:3467786/1‑150 (MQ=255)
         gTTCGGTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTc                                                                                                                                              >  2:2121890/1‑150 (MQ=255)
         gTTCGGTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTc                                                                                                                                              >  2:4031922/1‑150 (MQ=255)
           tCGGTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGc                                                                                                                                            >  2:4567080/1‑150 (MQ=255)
           aacgTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATATCACTTTTGCCACCTATATATCGACTTCCGCGACATCTTGTTCTGCTCAAATTAAAACGCCAACTTTAATGTGTTTTGAATCGCTGCGACAATGTTCGc                                                                                                                                            >  2:2472334/4‑150 (MQ=255)
                       aaCCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCcagca                                                                                                                                >  1:3862300/1‑150 (MQ=255)
                           gCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAAt                                                                                                                            >  1:1169009/1‑150 (MQ=255)
                           gCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAAt                                                                                                                            <  2:555146/150‑1 (MQ=255)
                             tAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTc                                                                                                                          >  1:1521620/1‑150 (MQ=255)
                             tAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTc                                                                                                                          >  1:704997/1‑150 (MQ=255)
                             tAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTc                                                                                                                          >  1:704717/1‑150 (MQ=255)
                                   cTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCa                                                                                                                    <  1:1334736/150‑1 (MQ=255)
                                        aGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCaat                                                                                                               <  1:3833062/150‑1 (MQ=255)
                                          ttCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCaataa                                                                                                             <  1:2683005/150‑1 (MQ=255)
                                          ttCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCaataa                                                                                                             <  1:2683067/150‑1 (MQ=255)
                                             aaCATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGc                                                                                                          >  1:2986693/1‑150 (MQ=255)
                                                aTCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGc                                                                                                       >  1:46213/1‑150 (MQ=255)
                                                   gAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGcccc                                                                                                    >  2:3172452/1‑150 (MQ=255)
                                                     ggCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTg                                                                                                  <  2:4736821/150‑1 (MQ=255)
                                                       ccATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTc                                                                                                >  1:3492733/1‑150 (MQ=255)
                                                        cATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTcc                                                                                               <  2:3712670/150‑1 (MQ=255)
                                                         atatTTCTCTTTTGCCACCTTTTTAGAGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCAAACTATAATTTGTTTTGCATCGCTGCGACTCTGTTCGCAACCTGCCAGCAGAATTCCGCCCAGCAATAATGCGGCCCCTGTCCg                                                                                              >  1:2881094/1‑150 (MQ=255)
                                                          tatTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGt                                                                                             <  1:3340678/150‑1 (MQ=255)
                                                           atTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTa                                                                                            >  2:1844030/1‑150 (MQ=255)
                                                           atTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTa                                                                                            >  2:769347/1‑150 (MQ=255)
                                                             ttCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAga                                                                                          <  2:1067661/150‑1 (MQ=255)
                                                              tctcTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAgat                                                                                         <  1:2356207/150‑1 (MQ=255)
                                                              tctcTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAgat                                                                                         <  1:2376656/150‑1 (MQ=255)
                                                              tctcTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAgat                                                                                         >  1:3713222/1‑150 (MQ=255)
                                                               ctctTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAgatg                                                                                        >  1:1473961/1‑150 (MQ=255)
                                                                  ttttGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGatgatg                                                                                     <  1:1091926/150‑1 (MQ=255)
                                                                     tGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATgtgt                                                                                  <  1:4818439/150‑1 (MQ=255)
                                                                          ccTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCa                                                                             <  1:2941783/150‑1 (MQ=255)
                                                                           cTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAg                                                                            <  2:4867682/150‑1 (MQ=255)
                                                                            tttttAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGt                                                                           <  1:2179354/150‑1 (MQ=255)
                                                                            tttttAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGt                                                                           <  1:2826911/150‑1 (MQ=255)
                                                                            tttttAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGt                                                                           <  1:670000/150‑1 (MQ=255)
                                                                             ttttAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGtt                                                                          <  1:3111203/150‑1 (MQ=255)
                                                                             ttttAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGtt                                                                          >  2:1061530/1‑150 (MQ=255)
                                                                              tttAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGttt                                                                         <  2:2845811/150‑1 (MQ=255)
                                                                               ttAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTc                                                                        <  1:2659264/150‑1 (MQ=255)
                                                                                tAGCGACTTCCGCGACATCTTATTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCa                                                                       <  1:1416612/150‑1 (MQ=255)
                                                                                   cGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCAttt                                                                    <  1:442947/150‑1 (MQ=255)
                                                                                    gACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCAtttt                                                                   <  1:4419848/150‑1 (MQ=255)
                                                                                    gACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCAtttt                                                                   <  1:1968424/150‑1 (MQ=255)
                                                                                    gACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCAtttt                                                                   >  1:607473/1‑150 (MQ=255)
                                                                                     aCTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCAttttt                                                                  <  1:4268878/150‑1 (MQ=255)
                                                                                     aCTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCAttttt                                                                  <  1:1815868/150‑1 (MQ=255)
                                                                                       ttCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCAttttttt                                                                <  2:4695983/150‑1 (MQ=255)
                                                                                        tCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTa                                                               <  2:141048/150‑1 (MQ=255)
                                                                                          cgcgACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATc                                                             >  1:2589069/1‑150 (MQ=255)
                                                                                          cgcgACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATc                                                             <  1:539095/150‑1 (MQ=255)
                                                                                            cgACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCt                                                           >  1:952022/1‑150 (MQ=255)
                                                                                            cgACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCt                                                           <  1:3592669/150‑1 (MQ=255)
                                                                                             gACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCtt                                                          >  2:2795064/1‑150 (MQ=255)
                                                                                              aCATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCttt                                                         <  1:3762656/150‑1 (MQ=255)
                                                                                                aTCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTtta                                                       >  2:3830433/1‑150 (MQ=255)
                                                                                                aTCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTGTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTtta                                                       <  2:4009999/150‑1 (MQ=255)
                                                                                                 tCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTttat                                                      <  1:3661300/150‑1 (MQ=255)
                                                                                                 tCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTttat                                                      >  1:121088/1‑150 (MQ=255)
                                                                                                  cTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTttatt                                                     <  2:3080776/150‑1 (MQ=255)
                                                                                                     gTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAAt                                                  <  1:3275131/150‑1 (MQ=255)
                                                                                                     gTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAAt                                                  >  1:3951899/1‑150 (MQ=255)
                                                                                                     gTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAAt                                                  >  1:2480313/1‑150 (MQ=255)
                                                                                                      ttCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCATCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATg                                                 <  2:3260625/150‑1 (MQ=255)
                                                                                                       tCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGa                                                <  2:3403858/150‑1 (MQ=255)
                                                                                                         tGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGAtt                                              <  2:2707403/150‑1 (MQ=255)
                                                                                                         tGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGAtt                                              <  2:3115869/150‑1 (MQ=255)
                                                                                                              cATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTg                                         <  1:827586/150‑1 (MQ=255)
                                                                                                                 ttATAACGCCAACTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACt                                      <  2:3552155/150‑1 (MQ=255)
                                                                                                                   ataACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACtat                                    <  1:4333246/150‑1 (MQ=255)
                                                                                                                    taACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACtata                                   <  2:1599819/150‑1 (MQ=255)
                                                                                                                       cGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCg                                >  1:4104804/1‑150 (MQ=255)
                                                                                                                          cAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTa                             >  2:3881598/1‑150 (MQ=255)
                                                                                                                           aaCTTTAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTACCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTaa                            <  1:1702883/150‑1 (MQ=255)
                                                                                                                             cTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATTGCGGTAAgg                          >  2:2980433/1‑150 (MQ=255)
                                                                                                                             cTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAgg                          <  2:1700013/150‑1 (MQ=255)
                                                                                                                               ttAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTc                        <  1:2241284/150‑1 (MQ=255)
                                                                                                                                tAATGTGTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCg                       <  1:506404/150‑1 (MQ=255)
                                                                                                                                 aaTGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGc                      >  1:1717656/1‑150 (MQ=255)
                                                                                                                                  aTGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCt                     >  2:4443859/1‑150 (MQ=255)
                                                                                                                                   tgtgTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTg                    >  2:1974060/1‑150 (MQ=255)
                                                                                                                                   tgtgTTTTGCATCGCTGCGACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTg                    >  2:1880759/1‑150 (MQ=255)
                                                                                                                                       ttttGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGt                >  1:4189019/1‑150 (MQ=255)
                                                                                                                                        tttGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGtt               >  2:708611/1‑150 (MQ=255)
                                                                                                                                          tGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTa             <  2:1245645/150‑1 (MQ=255)
                                                                                                                                          tGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTa             <  2:2986693/150‑1 (MQ=255)
                                                                                                                                            cATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTaaa           <  1:2540942/150‑1 (MQ=255)
                                                                                                                                            cATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTaaa           <  2:2748960/150‑1 (MQ=255)
                                                                                                                                            cATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTaaa           <  2:2339598/150‑1 (MQ=255)
                                                                                                                                             aTCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTaaaa          >  1:2949898/1‑150 (MQ=255)
                                                                                                                                             aTCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTaaaa          >  1:1351942/1‑150 (MQ=255)
                                                                                                                                             aTCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTaaaa          >  1:583510/1‑150 (MQ=255)
                                                                                                                                             aTCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTaaaa          >  1:1969055/1‑150 (MQ=255)
                                                                                                                                                  tgctACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTAAAaaacga     <  2:1642930/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTCGCATCGTTCGGTAGCAACGAACCGCTAAATCCTACCAGTTCAACATCGAGGCCATATTTCTCTTTTGCCACCTTTTTAGCGACTTCCGCGACATCTTGTTCTGCGCCATTTATAACGCCAACTTTAATGTGTTTTGCATCGCTGCTACTCTGGTCGCAACCTGCCAGCAGAATTCCGGCCAGCAATAATGCGGCCCCTGTCCGTAGATGATGTGTTGTCAGTTTCATTTTTTTATCCTTTTATTAATGATTACCTGACTATAGCGGTAAGGTCGCTGCGGTTTAAAAAACGAAAC  >  CP009273/3833123‑3833421

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: