Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,759,556 C→T 17.9% P94S (CCT→TCT)  rnlA → CP4‑57 prophage; RNase LS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,759,5560CT17.9% 558.1 / 104.3 218P94S (CCT→TCT) rnlACP4‑57 prophage; RNase LS
Reads supporting (aligned to +/- strand):  ref base C (74/105);  new base T (20/19);  total (94/124)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 220 total aligned reads displayed.

TTCGGGTCGTATTTACCAGACCAGGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATT  >  CP009273/2759407‑2759702
                                                                                                                                                     |                                                                                                                                                  
   gggTCGTATTTACCAGACCAGGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTctg                                                                                                                                                 <  1:4445279/150‑1 (MQ=255)
     gTCGTATTTACCAGACCAGGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCctgct                                                                                                                                               >  2:3489301/1‑150 (MQ=255)
           tttACCAGACCAGGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAAttt                                                                                                                                         <  1:4292541/150‑1 (MQ=255)
            ttACCAGACCAGGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTg                                                                                                                                        >  1:1115171/1‑150 (MQ=255)
             tACCAGACCAGGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGa                                                                                                                                       >  1:2567156/1‑150 (MQ=255)
              accagaccagGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAg                                                                                                                                      <  1:1232989/150‑1 (MQ=255)
                 agaccagGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAg                                                                                                                                   <  1:1949652/150‑1 (MQ=255)
                    ccagGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTa                                                                                                                                <  2:411419/150‑1 (MQ=255)
                     cagGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTaa                                                                                                                               <  2:4318450/150‑1 (MQ=255)
                      agGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTaaa                                                                                                                              >  1:2673370/1‑150 (MQ=255)
                         ggTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAatat                                                                                                                           <  2:2574970/150‑1 (MQ=255)
                           ttAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATgg                                                                                                                         <  1:2370652/150‑1 (MQ=255)
                                ttGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTg                                                                                                                    >  2:4424433/1‑150 (MQ=255)
                                  ggCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCa                                                                                                                  <  2:4402172/150‑1 (MQ=255)
                                   gCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCaa                                                                                                                 <  1:2739624/150‑1 (MQ=255)
                                    cAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAg                                                                                                                >  1:2104214/1‑150 (MQ=255)
                                    cAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAg                                                                                                                >  1:756187/1‑150 (MQ=255)
                                    cAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAg                                                                                                                >  1:422583/1‑150 (MQ=255)
                                         gtgGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGAttt                                                                                                           <  1:3444455/150‑1 (MQ=255)
                                          tgGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTg                                                                                                          >  1:1424290/1‑150 (MQ=255)
                                             aCATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTag                                                                                                       >  2:3757976/1‑150 (MQ=255)
                                              cATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTaga                                                                                                      <  1:457112/150‑1 (MQ=255)
                                                 tttttATAATGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGAc                                                                                                   <  1:3645027/150‑1 (MQ=255)
                                                  ttttATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACa                                                                                                  >  2:1302666/1‑150 (MQ=255)
                                                  ttttATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACa                                                                                                  <  1:388255/150‑1 (MQ=255)
                                                  ttttATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACa                                                                                                  <  1:387970/150‑1 (MQ=255)
                                                   tttATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACaa                                                                                                 <  2:4268065/150‑1 (MQ=255)
                                                    ttATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAg                                                                                                >  2:859394/1‑150 (MQ=255)
                                                    ttATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAg                                                                                                <  1:38307/150‑1 (MQ=255)
                                                     tataACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGc                                                                                               >  1:2962871/1‑150 (MQ=255)
                                                       taACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGt                                                                                             >  2:3172287/1‑150 (MQ=255)
                                                           ggggACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTa                                                                                         >  2:1717538/1‑150 (MQ=255)
                                                           ggggACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTa                                                                                         <  1:3702899/150‑1 (MQ=255)
                                                           ggggACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTa                                                                                         >  2:3073209/1‑150 (MQ=255)
                                                            gggACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTAc                                                                                        <  1:1733709/150‑1 (MQ=255)
                                                            gggACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTAc                                                                                        >  1:1394358/1‑150 (MQ=255)
                                                              gACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCt                                                                                      <  2:1846296/150‑1 (MQ=255)
                                                              gACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCt                                                                                      <  2:2708053/150‑1 (MQ=255)
                                                                cGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGt                                                                                    <  2:1023153/150‑1 (MQ=255)
                                                                 ggTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTa                                                                                   >  1:4553874/1‑150 (MQ=255)
                                                                 ggTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTa                                                                                   >  1:3362452/1‑150 (MQ=255)
                                                                   ttCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACt                                                                                 <  2:815330/150‑1 (MQ=255)
                                                                     cGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTg                                                                               <  2:647787/150‑1 (MQ=255)
                                                                      gACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGa                                                                              <  1:2545303/150‑1 (MQ=255)
                                                                         tACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCt                                                                           <  1:2797781/150‑1 (MQ=255)
                                                                          aCAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTa                                                                          <  1:4103436/150‑1 (MQ=255)
                                                                          aCAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTa                                                                          >  1:527237/1‑150 (MQ=255)
                                                                           cAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTAt                                                                         <  2:1441662/150‑1 (MQ=255)
                                                                             aTTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATca                                                                       <  2:2993849/150‑1 (MQ=255)
                                                                             aTTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATca                                                                       <  2:1530266/150‑1 (MQ=255)
                                                                              ttCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATcag                                                                      <  1:462918/150‑1 (MQ=255)
                                                                              ttCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATcag                                                                      <  1:463062/150‑1 (MQ=255)
                                                                                   tatCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGAt                                                                 <  1:2176711/150‑1 (MQ=255)
                                                                                    atCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATg                                                                >  2:759866/1‑150 (MQ=255)
                                                                                    atCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATg                                                                >  2:1137532/1‑150 (MQ=255)
                                                                                    atCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATg                                                                >  2:3597396/1‑150 (MQ=255)
                                                                                     tctcACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGa                                                               <  1:1611063/150‑1 (MQ=255)
                                                                                          cTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGc                                                          <  2:778679/150‑1 (MQ=255)
                                                                                          cTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGc                                                          >  1:2949958/1‑150 (MQ=255)
                                                                                           tGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCa                                                         >  2:3923494/1‑150 (MQ=255)
                                                                                            ggAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCat                                                        <  2:3887201/150‑1 (MQ=255)
                                                                                              aGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCatat                                                      <  1:180335/150‑1 (MQ=255)
                                                                                               gCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCAtata                                                     <  1:3565895/150‑1 (MQ=255)
                                                                                                ccAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATag                                                    <  2:1091725/150‑1 (MQ=255)
                                                                                                   aTCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGt                                                 <  1:3979757/150‑1 (MQ=255)
                                                                                                   aTCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGt                                                 <  2:2172907/150‑1 (MQ=255)
                                                                                                   aTCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGt                                                 >  2:3576693/1‑150 (MQ=255)
                                                                                                    tCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGtt                                                <  1:493182/150‑1 (MQ=255)
                                                                                                    tCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACAATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGtt                                                <  1:492158/150‑1 (MQ=255)
                                                                                                     cGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTc                                               >  2:1244833/1‑150 (MQ=255)
                                                                                                      gTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCa                                              <  2:4515912/150‑1 (MQ=255)
                                                                                                       ttCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCag                                             <  1:3139702/150‑1 (MQ=255)
                                                                                                        tctcTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCaga                                            >  1:3579617/1‑150 (MQ=255)
                                                                                                        tctcTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCaga                                            >  2:2105022/1‑150 (MQ=255)
                                                                                                        tctcTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCaga                                            <  2:1178028/150‑1 (MQ=255)
                                                                                                         ctctGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCagag                                           >  1:2077765/1‑150 (MQ=255)
                                                                                                             gggCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAaca                                       <  1:608304/150‑1 (MQ=255)
                                                                                                              ggCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAacac                                      >  2:2836171/1‑150 (MQ=255)
                                                                                                              ggCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAacac                                      >  2:468929/1‑150 (MQ=255)
                                                                                                                 cAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACActct                                   >  2:3680609/1‑150 (MQ=255)
                                                                                                                   agagTTAGCGGATCATCTTTTTGAAACCATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCg                                 >  1:3375554/1‑150 (MQ=255)
                                                                                                                   agagTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCg                                 >  2:2343938/1‑150 (MQ=255)
                                                                                                                     agTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTa                               <  1:3489301/150‑1 (MQ=255)
                                                                                                                      gTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTaa                              <  1:3217891/150‑1 (MQ=255)
                                                                                                                      gTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTaa                              >  1:3028102/1‑150 (MQ=255)
                                                                                                                       ttAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAAc                             >  1:4136815/1‑150 (MQ=255)
                                                                                                                       ttAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAAc                             >  1:1680279/1‑150 (MQ=255)
                                                                                                                          gctgATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCt                          <  1:1235725/147‑1 (MQ=255)
                                                                                                                              atcatcTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATa                      <  2:1891399/150‑1 (MQ=255)
                                                                                                                                catcTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATAcc                    <  1:32613/150‑1 (MQ=255)
                                                                                                                                 atcTTTTTGAAACCATCACTCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCg                   <  1:3458862/150‑1 (MQ=255)
                                                                                                                                 atcTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCg                   <  1:3515483/150‑1 (MQ=255)
                                                                                                                                      tttGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAgtg              <  1:1592537/150‑1 (MQ=255)
                                                                                                                                         gAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTgg           >  2:3543810/1‑150 (MQ=255)
                                                                                                                                           aaCCATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGaa         <  2:4773395/150‑1 (MQ=255)
                                                                                                                                             ccATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGaaaa       >  1:4691876/1‑150 (MQ=255)
                                                                                                                                              cATCAATTCTGCTGAATTTGAGCAGGTAAATATGGTACTGCCAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATCTAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAat      >  2:3269412/1‑150 (MQ=255)
                                                                                                                                              cATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAat      >  2:4555253/1‑150 (MQ=255)
                                                                                                                                               aTCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAatt     <  1:4088591/150‑1 (MQ=255)
                                                                                                                                                  aaTCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAattatt  <  2:4407822/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTCGGGTCGTATTTACCAGACCAGGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGGCCAAGAGTTAGCGGATCATCTTTTTGAAACCATCAATCCTGCTGAATTTGAGCAGGTAAATATGGTACTGCAAGGATTTGTAGAGACAAGCGTTCTACCTGTACTTGAGCTATCAGCAGATGAATCGCATATAGAGTTCAGAGAACACTCTCGTAACGCTCATACCGTAGTGTGGAAAATTATT  >  CP009273/2759407‑2759702

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: