Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 449,933 G→A 11.3% M2I (ATG→ATA bolA → stationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273449,9330GA11.3% 116.8 / 6.1 49M2I (ATG→ATAbolAstationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base A (4/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CCCGCCATCATAACATTGCCAACGGCGAGGGGAAGTGGGTAAGGCATGTAAATTCATCATGTTGACGAAATAATCGCCCCTGGTAAAAGAAACACTGATGCGAGGCCTGTGTTTCAATCTTTAAATCAGTAAACTTCATACGCTTGACGGAAAAACCAGGACGAAACCTAAATATTTGTTGTTAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGAACTCTCTACTACCGTTCATGCGCTGGCTCTGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTT  >  CP009273/449666‑450205
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                
cccGCCATCATAACATTGCCAACGGCGAGGGGAAGGGGGTAAAGCATGTAAGTTCATCATGTTGACGAAATAATAGCACCGGGTAAAAGAAACACTGATGCGAGGCCTGTGTTTAAATCTTTAAGTCAGTAAACTTAATACGCTTGAACGACACACCAGGCCGAATCCTAACTATTAGTTGTAACGATACAATGGAAACGATCAACGCGAATAGACTTTAAAGGGATGGCTGACATCTCAGCGAGGTCGGAGGAGCTATTTCAGGAGTCTACGTg                                                                                                                                                                                                                                                                           >  2:330829/1‑275 (MQ=255)
                 gCCAACGGCGAGGGGAAGTGGGTAAGGCATGTAAATTCATCATGTTGACGAAATAATCGCCCCTGGTAAAAGAAACACTGATGCGAGGCCTGTGTTTCAATCTTTAAATCAGTAAACTTCATACGCTTGACGGAAAACCCAGGACAAAACCTAAATATTTGTTGTTAATCTGCAATGGAAACGGTAAAAGCGGCTAGTATTCAAAGGGATGAATAACATCTCAGAGTTGTCGGAGGAAATATTTCATGATAATACGATAGCGCATAGAaaaaaa                                                                                                                                                                                                                                                           >  1:763015/1‑274 (MQ=255)
                            agggAAATGGGTGNAGGCATTTATATTCATTATTGTTACGAAATAATCGCTCCTCGTAAAAGAAACACATATTCGATGCCTGTGTTTCAATCTTTAAATCAGTAAACTTCATACGCTTGACGGAAAAACCAAGACAAAAGCTAAATGTTTGTTGTTAAGCTGGAATGGAAACGGGAAAAGCGGCGAATATTTAAGGTGATGGAAGACGTCTCAGCGTTGTCGGAGGGGATATTGCATGATGACACGTGAGCGGATAGAAGAAAAAATAAGggcg                                                                                                                                                                                                                                                 <  2:595184/273‑1 (MQ=255)
                                     ggTAAGACATGTAACTTCATCATGTTGATGAAATAATCGCCCCTTGTAAAAGAAACACTGATGCGAGGCCTGTGTTTCAATCTTTAAATCAGTAAACTTCATACGCTTGACGGAAAAACCAGGACGAAACCTAAATATTTGTTGTTAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTANAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGCTAGAAGAAAAATTAAGGGCGACGTTACAAcc                                                                                                                                                                                                                                      >  1:595189/1‑275 (MQ=255)
                                                                                         gAAACAATGATGAGAGGCCTGTGTTTATATCTTTATATCAGTAAGCTTCATTCGCTTGACGGAAATACCAGGACGAAAGCTAATTATTTGTTGTTAAGCTGCAATGGAAACGGTAAAAGCGGCTTGTATTTACAGGGATGGATGTCATCTCAGCGTTGTCGGAGGAGATCTTTAATTATGAAACAGGAGCAGATAGTGCATAAATTAAGGGCGTCGTCGCAACCAGTATTCCTCGAAGTAGTGGATGAAAGCCACAGTCACAATGTCACAGCCgt                                                                                                                                                                                   <  2:571636/275‑2 (MQ=255)
                                                                                             cacTGAGGCGAGGCCTGTGTTTCAATCTTTAAATCAGTAAACTTCATACGCTGGACGGAAAAACCAGGACGAAACCTAAATATTTGTTGTTAAGCTGCAATGGAAACGGTAAAAGAGGCTAGTCTTAAAAGGGATGGATGACATATCAGCGTTGTTGGAGGAGATATTTCATGATAATACGTGAGAGAAGAGACGAAAAATGAAGGGCGACGATAGAAAACGTATAACGCGAAGGAGTGGATGAAAGCTATAGCAACAATGtacagggcgnactn                                                                                                                                                                               >  2:217/1‑262 (MQ=255)
                                                                                                  aTGCGAGGCCTGTGTTTCAATCTTTAAATCAGTAAACTTCATACGCTTGACGGAAAAACCAGGACGAAACCTAAATATTTGTTGTTAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAg                                                                                                                                                                          >  1:137099/1‑275 (MQ=255)
                                                                                                   tgtgAGTCCTGTGTTTCAATCTTTAAATCTGTAAACTTCATACGCTTGACGGAAAAACCAGGACGAAACCTAAATATTTGTTGTTAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGc                                                                                                                                                                         <  1:352507/272‑1 (MQ=255)
                                                                                                       aGGCCTGTGTTTCAATCTTTAAATCAGTAAACTTCATACGCTTGACGGAAAAACTAGGACGAAACCTAAATATTTGTTGTTAAGCTGCAATGGAAACGGTAAAAGCGGCCAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGAGAGAAGAAAAATTAAGGGAGGCGTTACAACCAGTATTCCTCGAAGTAGTGGATGAAGGATATGATCAAAATGTCGCAGCGGGCTCTGAAAGctttt                                                                                                                                                                     >  2:257513/1‑271 (MQ=255)
                                                                                                                          tAAACCAGTAAACTTCATACGCTTTACGGAAAAACCAGGACGAAACCTACATACTTGTTGTTAAGCTGCAATGGAAACGGTAAACACGCCTAGCATTTAAAGGGATGGAGGACATCTCAGTGGTGTCGGATGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCATTCTAACCCGTATTCATCGAAGAAGTGGATGAAAGATAGAGTCACAAAGTCCCGGCCGGCCCTGAAAGCCATTTAAAAgattttgttttctg                                                                                                                                                  >  2:626465/1‑262 (MQ=255)
                                                                                                                                                 tGACGGAAAAACCAGGACGAAACCTAGATATCGATTGCTCCGCTGCATTGGAAACGCATATAGCGGTTAGAATTTAAAGGCATGGTTGACATCTCAGGGCTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGACTCAGAAAGCCATTTTAAAGTTGTGCTGGTCAGAGATCGTTTTANCGGTGaaact                                                                                                                            >  2:66980/1‑271 (MQ=255)
                                                                                                                                                    cGGAAAAACCAGGACGAAACCTAAATATTTGTTGTTAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATAc                                                                                                                                                                                                                                                                              >  1:184736/1‑125 (MQ=255)
                                                                                                                                                        aaaaCCAGGACGAAACCTAAATATTTGTTGTTAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGa                                                                                                                                                                                                                                                                          <  2:184733/125‑1 (MQ=255)
                                                                                                                                                                   gAAACCTAAATATTTGTTGTTATGCTGCAATGGAAACGGTAAAAGCGGCTAGTAATTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGaa                                                                                                         <  1:66980/275‑1 (MQ=255)
                                                                                                                                                                           catatTTGTTGTTAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGAGGATAGAAGAAAAATCAAGGGCGGCGTTCCAACCCGTATTCGTCGAAGTAGTGGATGAACGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTAAAAGTTGTGCTGGTCAGTGATCGTTTTACGGGAGAACGTTTTTTGCATCTTAATTGAATGCTTTAc                                                                                                 >  2:368259/2‑275 (MQ=255)
                                                                                                                                                                              atTTGGTGTTTAGCGGTAATTGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGAGTACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAACGGATAGAAGAAAAATTAAGGGCGTCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAg                                                                                               <  1:217/274‑1 (MQ=255)
                                                                                                                                                                                 tgttgtTAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAACAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAACGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTAc                                                                                            >  1:43925/1‑274 (MQ=255)
                                                                                                                                                                                  gttgttAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTAt                                                                                                                                                                                                                                 >  1:772656/1‑140 (MQ=255)
                                                                                                                                                                                  gttgttAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATAATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGAt                                                                                                                                                                                                                >  2:44217/1‑157 (MQ=255)
                                                                                                                                                                                      ttAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCt                                                                                                                                                                                                                             <  2:772647/140‑1 (MQ=255)
                                                                                                                                                                                      ttAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATAATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGaaa                                                                                                                                                                                                            <  1:44217/157‑1 (MQ=255)
                                                                                                                                                                                             ggaaTGGAAAAGGTAAAAGCGGATATTATTTAAGGGGATGGATGACGTCTCAGCTTTGTCGGAGGATATATTTCATGATGCTACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGGTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGaa                                                                               <  1:30536/273‑1 (MQ=255)
                                                                                                                                                                                             gCAATGGAAACTGTAAATGCGGCTAGTATTTATAGGGATGGATGACATCTCAGCGTTGTCGGAGGGGATATTTCATAATGATACGTGGGCGGATAGAAGAAAAAGTAAGGGCGGCGTTCCAATCCGTACTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGaa                                                                               <  1:722802/275‑1 (MQ=255)
                                                                                                                                                                                             gCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAAGTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGaa                                                                               <  1:257517/275‑1 (MQ=255)
                                                                                                                                                                                              agaTGGAAACGGTAAAAGCGGCTATTATTTAAAGGGAGGGATGACATCTCAGTGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTACCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTTAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGaa                                                                               <  1:330833/272‑1 (MQ=255)
                                                                                                                                                                                                                  gCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGaaa                                                                                                                                                                                                                                                            >  1:760545/1‑81 (MQ=255)
                                                                                                                                                                                                                   ttAGTAATAAATGGGGTGGAGGACAGCTCAGCGTTGTCGGAGGAGATATTTCATGATTATACGTTAGCTGAGAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGGTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGAACTCTCTACTACCGTTCATgcgc                                                         <  2:43925/274‑1 (MQ=255)
                                                                                                                                                                                                                   gtcctctatgaaGGGGAGGATGACAGCGCAGCGTTGTCGGAGGAGATATTTCATGAGGATACGTGGGCGGATAGAAGAAAAATTAATGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCGCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGAACTCTCTACTACCGTTCATgcgc                                                         <  2:137099/265‑1 (MQ=255)
                                                                                                                                                                                                                      gTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAAtt                                                                                                                                                                                                                                                        <  2:760536/81‑1 (MQ=255)
                                                                                                                                                                                                                         tttAAAGGGATGGTTGTCATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGATCGGAGAGAAGAAAAATTAAGGGCGGCGTTNCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGGACTCTCTACTACCGTTCATGCGCTGGct                                                    <  1:626472/274‑1 (MQ=255)
                                                                                                                                                                                                                              aGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGaaa                                                                                                                                                                                                            >  1:48473/1‑117 (MQ=255)
                                                                                                                                                                                                                                 gatggatgACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCACAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGACTGATTTACAGTACTTTAGCGGAGGAACTCTCTACTACCGTTCATGCCATGGCTATGCATACt                                           >  1:16101/1‑275 (MQ=255)
                                                                                                                                                                                                                                 gatggatgACATATCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTACTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGATAAGCGATCGTTTAACGGGTGAACGTTTTCTGAATAGTCATCGAATGATATACAGTACTTTATCGGAGGAACTCACTAATACCGTTCATGCGCAGGCGATGCAAACt                                           >  1:61914/1‑275 (MQ=255)
                                                                                                                                                                                                                                  atggatgACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTa                                                                                                                                                                                                        <  2:48473/117‑1 (MQ=255)
                                                                                                                                                                                                                                   tggatgACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCg                                                                                                                                                >  2:657716/1‑172 (MQ=255)
                                                                                                                                                                                                                                   tggatgACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCAGGTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGAActctct                                                                         >  1:661261/1‑243 (MQ=255)
                                                                                                                                                                                                                                   tggatgACATCTCAGCGTTGTCGGAGGAGATATTTCATGATAATACGTGCGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGt                                                                                                                                                                                                    >  2:25312/1‑120 (MQ=255)
                                                                                                                                                                                                                                       tgACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCg                                                                                                                                            <  1:657723/172‑1 (MQ=255)
                                                                                                                                                                                                                                       tgACATCTCAGCGTTGTCGGAGGAGATATTTCATGATAATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTcaca                                                                                                                                                                                                <  1:25312/120‑1 (MQ=255)
                                                                                                                                                                                                                                          cATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAg                                                                                                                                                                                                                                                     >  1:20954/1‑64 (MQ=255)
                                                                                                                                                                                                                                            tctcAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGtt                                                                                                                                                             >  2:468067/1‑150 (MQ=255)
                                                                                                                                                                                                                                              tcAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGggcg                                                                                                                                                                                                                                                 <  2:20954/64‑1 (MQ=255)
                                                                                                                                                                                                                                                aGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGc                                                                                                                                                         <  1:468072/150‑1 (MQ=255)
                                                                                                                                                                                                                                                    ttGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAAcc                                                                                                                                                                                                                                      >  2:233110/1‑69 (MQ=255)
                                                                                                                                                                                                                                                     agacGGAGGAGATATTTCCTGATGATACGTGAGTGGATAGAAGAAAAATTAACGGCGGCGTTCCAACCCGTATTCCTCGa                                                                                                                                                                                                                          >  2:773709/4‑80 (MQ=255)
                                                                                                                                                                                                                                                        cGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTa                                                                                                                                                                                                                                  <  1:233114/69‑1 (MQ=255)
                                                                                                                                                                                                                                                         ggaggaGATATTTCATGATGATACGTGAGCGGCTAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAagta                                                                                                                                                                                                                      <  1:773718/80‑1 (MQ=255)
                                                                                                                                                                                                                                                           atgagaTATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGAACTCTCTACTACCGTTCATGCGCTGGCTCTGCATACTTaca                                       <  1:463604/251‑1 (MQ=255)
                                                                                                                                                                                                                                                                atatTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGAACTCTCTACTACCGTTCATGCGCTGGCTCTGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAg             >  1:271546/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                       atgatgATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGAACTCTCTACTACCGTTCATGCGCTGGCTCTGCATACt                                           >  1:288448/1‑237 (MQ=255)
                                                                                                                                                                                                                                                                           tgatACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGAACTCTCTACTACCGTTCATGCGCTGGCTCTGCATACTTaca                                       <  2:288444/237‑1 (MQ=255)
                                                                                                                                                                                                                                                                           tgatACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGAACTCTCTACTACCGTTCATGCGCTGGCTCTGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCtt  <  2:271542/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                
CCCGCCATCATAACATTGCCAACGGCGAGGGGAAGTGGGTAAGGCATGTAAATTCATCATGTTGACGAAATAATCGCCCCTGGTAAAAGAAACACTGATGCGAGGCCTGTGTTTCAATCTTTAAATCAGTAAACTTCATACGCTTGACGGAAAAACCAGGACGAAACCTAAATATTTGTTGTTAAGCTGCAATGGAAACGGTAAAAGCGGCTAGTATTTAAAGGGATGGATGACATCTCAGCGTTGTCGGAGGAGATATTTCATGATGATACGTGAGCGGATAGAAGAAAAATTAAGGGCGGCGTTCCAACCCGTATTCCTCGAAGTAGTGGATGAAAGCTATCGTCACAATGTCCCAGCCGGCTCTGAAAGCCATTTTAAAGTTGTGCTGGTCAGCGATCGTTTTACGGGTGAACGTTTTCTGAATCGTCATCGAATGATTTACAGTACTTTAGCGGAGGAACTCTCTACTACCGTTCATGCGCTGGCTCTGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTT  >  CP009273/449666‑450205

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