Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,019,631 G→T 31.1% M169I (ATG→ATT)  guaD → guanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,019,6310GT31.1% 85.9 / 46.7 72M169I (ATG→ATT) guaDguanine deaminase
Reads supporting (aligned to +/- strand):  ref base G (27/21);  new base T (8/14);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACA  >  CP009273/3019487‑3019777
                                                                                                                                                |                                                                                                                                                  
taTGGCGTCGTTCTTCATCATGCAGCTTGTACGTACCGGAATCACCACGGCTCTGTAGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGGGATTATgg                                                                                                                                                <  1:738740/148‑1 (MQ=255)
 atgccggtattcttcATCAAGTGGCATTAACGTAACTGAACCACCCCGTCGCCGGTGTTTGGCACTGTTCCTCCCCAATCTGTTGATGCCCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCCGTAAGGTGATGATGGa                                                                                                                                               <  1:828228/140‑1 (MQ=255)
 aggtccgggTTCTTCATCAAGCATCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCCTCCGCAATCTGTTGATGCGCTCTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGa                                                                                                                                               <  1:633038/141‑1 (MQ=255)
 aTGTCGGCGTTCTTCATCAGGCAGCTTTTACGTCACGGACCCACCACGGCGCGGGGGTTGGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGATGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGa                                                                                                                                               <  2:201818/149‑1 (MQ=255)
 aTGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTCTGAAGGAGCCATTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGa                                                                                                                                               <  1:791343/149‑1 (MQ=255)
    tCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCg                                                                                                                                            >  1:521264/1‑149 (MQ=255)
      ggCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCa                                                                                                                                          <  1:320662/149‑1 (MQ=255)
       gCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCaa                                                                                                                                         <  1:784733/149‑1 (MQ=255)
          ttcttcATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGc                                                                                                                                      <  1:721883/149‑1 (MQ=255)
          ttcttcATCAAGCAGCTTTTACGTAACGGAACACCCCCCCCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGc                                                                                                                                      <  2:162153/149‑1 (MQ=255)
            cttcATCAAGGCGCTTGTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCCTATCAATATGCGTAGGATTGCCGGTAAGGTGATGATGGACCGCTACGCAc                                                                                                                                    <  1:405787/149‑1 (MQ=255)
            cttcATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGAATGATTGCCGGAAAGGTGATGATGGACCgacacgaac                                                                                                                                    >  2:448115/1‑141 (MQ=255)
              tcatcaAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGAGTGCCGGTAAGGTGATGATGGACCGCAACGCACCg                                                                                                                                  >  2:693495/1‑149 (MQ=255)
                    agcagcTTTTACGGAACGGAACCAGCACGGCGCTTGGGGTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCCTATCCATATGCGTATATTTGCCGGTAAGGTGAGTATGGACCGCAACGCACCGGATTAt                                                                                                                            <  1:209361/149‑1 (MQ=255)
                     gcagcTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATc                                                                                                                           >  2:2784/1‑149 (MQ=255)
                         cTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCt                                                                                                                       >  2:666827/1‑149 (MQ=255)
                         cTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCGGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCCACGCACCGGATTATCTGcg                                                                                                                       >  1:758071/1‑148 (MQ=255)
                            ttACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGa                                                                                                                    >  2:818106/1‑149 (MQ=255)
                                gTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGTATTAAATGACCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACt                                                                                                                >  2:417323/1‑149 (MQ=255)
                                gTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACt                                                                                                                >  1:336981/1‑149 (MQ=255)
                                gTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTCTCTGCTCGacacc                                                                                                                >  2:191101/1‑148 (MQ=255)
                                gTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGGAAGGTGATCATGGACCGCCACGCACCGGATTATCTGCTAGACACt                                                                                                                >  1:272122/1‑149 (MQ=255)
                                 tAACGGACCCACCACGGCGCTGGTGTTTGGCTCTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTg                                                                                                               <  1:339415/149‑1 (MQ=255)
                                   aCGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGcc                                                                                                             >  2:616277/1‑149 (MQ=255)
                                   aCGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGcc                                                                                                             <  2:6545/149‑1 (MQ=255)
                                     ggAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGGTTATCTGCCCGACACTGCCGa                                                                                                           >  1:591080/1‑149 (MQ=255)
                                      ggcccgcctcGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTGGAGGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGaa                                                                                                          <  1:278472/140‑1 (MQ=255)
                                      gAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGaa                                                                                                          >  2:731767/1‑149 (MQ=255)
                                           accacGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTAGACACTGCCGAAagcag                                                                                                     >  2:676233/1‑149 (MQ=255)
                                                 gcgcTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTGTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGACGGACCGCAACGCACCGGATTATCTCCCCGCCACTGCCGGAAACAGCTCTCa                                                                                               >  1:42261/1‑149 (MQ=255)
                                                           ttGGCACTGTTCATCCGCACTCTGTTGCTGCGCTGTTGGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAACGTGCTTATGCAACGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACcaaagcaaa                                                                                     <  1:185953/149‑1 (MQ=255)
                                                           ttGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACcaaagcaaa                                                                                     <  2:376794/149‑1 (MQ=255)
                                                           ttGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACcaaagcaaa                                                                                     <  2:617674/149‑1 (MQ=255)
                                                             ggCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGa                                                                                   >  2:706293/1‑149 (MQ=255)
                                                                aCTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACt                                                                                >  1:226083/1‑149 (MQ=255)
                                                                aCTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACt                                                                                >  2:708811/1‑149 (MQ=255)
                                                                  tGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGAAAAGAACTGa                                                                              >  1:697442/1‑149 (MQ=255)
                                                                   gTTCATCAGCAGTCTGTTGATGAGCTGTTTGCAGCCGCCAGTCGTATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGAt                                                                             <  1:449893/149‑1 (MQ=255)
                                                                    ttCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATc                                                                            <  1:317498/149‑1 (MQ=255)
                                                                    ttCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATc                                                                            <  1:133848/149‑1 (MQ=255)
                                                                          cGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATCTCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACCCTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTgaacgaaccc                                                                      >  1:371550/1‑147 (MQ=255)
                                                                            cAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTg                                                                    >  1:83901/1‑149 (MQ=255)
                                                                            cAATCTGGTGTTGCGCTGTTTGAAGCCACCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTg                                                                    <  2:66238/149‑1 (MQ=255)
                                                                             aaTCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTgg                                                                   <  2:765734/149‑1 (MQ=255)
                                                                                      aTGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAAt                                                          >  1:10710/1‑149 (MQ=255)
                                                                                           cTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgtcg                                                     >  2:269908/1‑149 (MQ=255)
                                                                                               ttGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTg                                                 <  2:336981/149‑1 (MQ=255)
                                                                                                tGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGc                                                >  1:599535/1‑149 (MQ=255)
                                                                                                 gAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCt                                               <  1:2784/149‑1 (MQ=255)
                                                                                                   aGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCtat                                             >  2:732215/1‑149 (MQ=255)
                                                                                                     ccgccAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTatat                                           >  2:110163/1‑149 (MQ=255)
                                                                                                         cAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGGTAAGTTTTAGCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGa                                       <  2:417055/149‑1 (MQ=255)
                                                                                                         cAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGa                                       <  2:468415/149‑1 (MQ=255)
                                                                                                         cAGGCATATCAATATGCGTATGATTTCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGa                                       <  1:26387/149‑1 (MQ=255)
                                                                                                            tCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTa                                    >  2:328369/1‑149 (MQ=255)
                                                                                                               tatCAATATGCGTATGATTTCCGGTAAGGTAATGATGGACCGCAACGCACCGGATTATCGGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGGCTGCTATATGCGATTACGc                                 >  2:651618/1‑149 (MQ=255)
                                                                                                               tatCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGc                                 >  2:20920/1‑149 (MQ=255)
                                                                                                                 tccatttGCGTATGATTGCCGGTAAGGCGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCCCCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCa                               <  1:680824/143‑1 (MQ=255)
                                                                                                                 tCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCAGCCATCTGCGATTAccccc                               >  1:155109/1‑148 (MQ=255)
                                                                                                                  cAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCAc                              >  1:733167/1‑149 (MQ=255)
                                                                                                                    atatGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATAGGCGATGACGCCACGc                            >  2:199970/1‑149 (MQ=255)
                                                                                                                     tatGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCGGCTATATGCGGTTACGCCACGCt                           >  1:808758/1‑149 (MQ=255)
                                                                                                                      atGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCtt                          <  2:277987/149‑1 (MQ=255)
                                                                                                                       tGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCCGCTATCACCAAAGCAGAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCCATATGCGGTTCCGCCCCGCTTc                         >  1:340751/1‑149 (MQ=255)
                                                                                                                       ggCGTATGATTGCCTGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTc                         <  2:805916/148‑1 (MQ=255)
                                                                                                                            cttaTTGCCGGTAAGGGGTTCATGGACCGCAACCAACCGGTTTATCTGCTCGACACTGCCGAAAGCAGCTATCACCACAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGcccc                    <  2:57923/146‑1 (MQ=255)
                                                                                                                             tGATTGCCGGTAAGGTGATTCTGGACCGCAACGCTCCCGATTATCTGCCCGACACTGCCGAATGCATCTATCACCAAAGCAAAGTTCTGATCGATCGCTTGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCg                   <  1:818106/149‑1 (MQ=255)
                                                                                                                                 tGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCt               >  2:698927/1‑149 (MQ=255)
                                                                                                                                   ccGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCa             <  2:212253/149‑1 (MQ=255)
                                                                                                                                       tAAGGTGCTTATGGACCGCAACGCACCGGATTATCTGCTCGACGCTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCAtctc         <  1:180356/149‑1 (MQ=255)
                                                                                                                                         aGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCt       <  2:653961/149‑1 (MQ=255)
                                                                                                                                           gTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGa     <  2:797828/149‑1 (MQ=255)
                                                                                                                                              atgatgGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACa  <  2:697442/149‑1 (MQ=255)
                                                                                                                                              atgatgGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCAATTACGCCACCCTTCGCCCCGACCTCATCTCCTGAACa  >  2:453160/1‑149 (MQ=255)
                                                                                                                                              atgaggGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATTAAATATTACGCCACGATTCGCCCCGACCTCATCTCCTGAACa  <  1:148925/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAACA  >  CP009273/3019487‑3019777

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑