Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,285,852 A→T 5.6% I19N (ATT→AAT)  PP_2887 ← catalase‑like peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,285,8520AT 5.6% 216.0 / 4.0 71I19N (ATT→AAT) PP_2887catalase‑like peroxidase
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base T (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

GCCCTTGGCGTGGTTCTTGCGGTAGCCTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCGGGAATAGAATTTGTTCACAGGCGTCTGC  >  NC_002947/3285704‑3285994
                                                                                                                                                     |                                                                                                                                              
gCCCTTGGCGTGGTTCTTGCGGTAGCCTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAAt                                                                                                                                               <  2:1001166/150‑1 (MQ=255)
gCCCTTGGCGTGGTTCTTGCGGTAGCCTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAAt                                                                                                                                               >  1:345685/1‑150 (MQ=255)
  ccTTGGCGTGGTTCTTGCGGTAGCCTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATgg                                                                                                                                             >  2:410213/1‑150 (MQ=255)
   cTTGGCGTGGTTCTTGCGGTAGCCTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATggc                                                                                                                                            >  2:1329190/1‑150 (MQ=255)
   cTTGGCGTGGTTCTTGCGGTAGCCTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATggc                                                                                                                                            >  1:187476/1‑150 (MQ=255)
         gTGGTTCTTGCGGTAGCCTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGCCACCTGCGGCCAGCAGCCCGCCGCCAATCGCGcccgg                                                                                                                                      >  2:870209/1‑148 (MQ=255)
               cTTGCGGTAGCCTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGgccgga                                                                                                                                <  1:410213/150‑5 (MQ=255)
                   cGGTAGCCTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTC‑GGGGGCGCAGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCtt                                                                                                                           <  2:1169665/150‑1 (MQ=255)
                     gTAGCCTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCttt                                                                                                                          >  2:1230286/1‑150 (MQ=255)
                          cTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTggccg                                                                                                                     >  2:103318/1‑150 (MQ=255)
                             ggTAATGCCCGGCCTGGGCTTCGAAGAAGACGATGATGCGCATTTGCGGGAGTCGTGGGGGGCCGAGCCAGCCAGC‑CCAGCGATCACCTGGGGCTTTGTCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTggccgggc                                                                                                                  >  1:325185/1‑150 (MQ=255)
                                    cccGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCaa                                                                                                           >  2:504841/1‑150 (MQ=255)
                                      cGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAgg                                                                                                         <  1:1329190/150‑1 (MQ=255)
                                      cGGCCTGGGCTTCGAAGGTGTCGATGATGCGCGGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAgg                                                                                                         <  2:123111/150‑1 (MQ=255)
                                        gCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGgt                                                                                                       <  1:395594/150‑1 (MQ=255)
                                             ggCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGtt                                                                                                  >  1:1337413/1‑150 (MQ=255)
                                               cTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCGATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCa                                                                                                >  1:224000/1‑150 (MQ=255)
                                                ttCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAgg                                                                                                         <  1:1417185/140‑1 (MQ=255)
                                                ttCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAgg                                                                                                         >  2:1417185/1‑140 (MQ=255)
                                                ttCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCAt                                                                                               <  2:325185/150‑1 (MQ=255)
                                                ttCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCAt                                                                                               >  2:599086/1‑150 (MQ=255)
                                                   gAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATggcg                                                                                                                                           <  1:1216127/103‑1 (MQ=255)
                                                   gAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATggcg                                                                                                                                           >  2:1216127/1‑103 (MQ=255)
                                                   gAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGt                                                                                            >  2:985195/1‑150 (MQ=255)
                                                       gtgtCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATc                                                                                        <  1:932565/150‑1 (MQ=255)
                                                          tCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTggccggg                                                                                                                   <  2:1114360/120‑1 (MQ=255)
                                                          tCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTggccggg                                                                                                                   >  1:1114360/1‑120 (MQ=255)
                                                           cGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTcc                                                                                    <  1:985195/150‑1 (MQ=255)
                                                            gatgatGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGGGGCGTAG‑TCAAAGTCGGGACCTGCGGCCAGCAGCGCGCCGCCAATGGCGGCCAGGAGCAGGTCTTTGGCCGGGTCGTGCAAGGGTGTGTTCATGGTCATCGTCCg                                                                                   >  1:519606/1‑150 (MQ=255)
                                                            gatgatGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTcc                                                                                    <  1:1329069/149‑1 (MQ=255)
                                                            gatgatGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTcc                                                                                    >  2:1329069/1‑149 (MQ=255)
                                                            gatgatGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCg                                                                                   >  1:1322615/1‑150 (MQ=255)
                                                            gatgatGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGAAGACCCCTGGCAATTGTGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCg                                                                                   >  1:324543/1‑150 (MQ=255)
                                                                 tgCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCa                                                                              <  1:1230286/150‑1 (MQ=255)
                                                                  gcgcTGTGGAGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAg                                                                             <  1:833013/150‑1 (MQ=255)
                                                                             tCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGAGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCa                                                                  <  1:870209/150‑1 (MQ=255)
                                                                               aGTCGGGGTGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGtt                                                                                <  1:72271/134‑1 (MQ=255)
                                                                               aGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGtt                                                                                >  2:72271/1‑134 (MQ=255)
                                                                                   ggggggCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTatga                                                            <  1:991208/150‑1 (MQ=255)
                                                                                      gggCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACgt                                                         <  2:187476/150‑1 (MQ=255)
                                                                                      gggCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAACAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGAGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACgt                                                         >  1:328451/1‑150 (MQ=255)
                                                                                        gcgcCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGgt                                                                                                       >  1:765552/1‑102 (MQ=255)
                                                                                        gcgcCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGgt                                                                                                       <  2:765552/102‑1 (MQ=255)
                                                                                            cGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCaga                                                   <  1:599086/150‑1 (MQ=255)
                                                                                            cGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCaga                                                   <  2:716209/150‑1 (MQ=255)
                                                                                             gAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCagag                                                  <  2:345685/150‑1 (MQ=255)
                                                                                                          gCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGGCAGGACGACGTGTCAGAGAAAAAGTTATTct                                     >  1:1276885/1‑149 (MQ=255)
                                                                                                            aTAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTggccgg                                                                                                                    <  1:1007504/69‑1 (MQ=255)
                                                                                                            aTAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTggccgg                                                                                                                    >  2:1007504/1‑69 (MQ=255)
                                                                                                              aG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGa                                 >  1:757457/1‑150 (MQ=255)
                                                                                                               g‑gCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGaaa                                               <  1:152784/135‑1 (MQ=255)
                                                                                                               g‑gCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGaaa                                               >  2:152784/1‑135 (MQ=255)
                                                                                                                 gCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATc                               >  2:1364373/1‑150 (MQ=255)
                                                                                                                  cAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCg                              >  1:860149/1‑150 (MQ=255)
                                                                                                                    aaGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCggg                            <  2:1337413/150‑1 (MQ=255)
                                                                                                                    aaGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCggg                            <  2:860149/150‑1 (MQ=255)
                                                                                                                        cGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCGGGAATa                        <  1:504841/150‑1 (MQ=255)
                                                                                                                             ccTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCGGGAATAGAAtt                   <  2:1322615/150‑1 (MQ=255)
                                                                                                                             ccTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCGGGAATAGAAtt                   >  2:646331/1‑150 (MQ=255)
                                                                                                                                    ccAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAg                                                                 >  2:1074591/1‑97 (MQ=255)
                                                                                                                                    ccAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAg                                                                 <  1:1074591/97‑1 (MQ=255)
                                                                                                                                     cagcagCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCGGGAATAGAATTTGTTcaca           <  2:324543/150‑1 (MQ=255)
                                                                                                                                      agctgcACGCCGCCAATGGCGGCCAGCCGCAAGGCTTAGGCCGGGCAGAGCAACAGTCTCTTAATGGTAAACCTCCCACGCTGCGCAATTGATATTTGGACTTGTCTTAGAAATATTTAGTCCTTGATCGTGAATAGAATGTGTTCACAg          <  2:519606/150‑1 (MQ=255)
                                                                                                                                        cagcaCGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCGGGAATAGAATTTGTTCACAGGc        >  1:511283/1‑150 (MQ=255)
                                                                                                                                           ctggccgccATTGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCt                                   <  1:877408/117‑1 (MQ=255)
                                                                                                                                           ctggccgccATTGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCt                                   >  2:877408/4‑120 (MQ=255)
                                                                                                                                            tggccgccATTGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCGGGAATAGAATTTg                 <  2:39292/135‑1 (MQ=255)
                                                                                                                                            tggccgccATTGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCGGGAATAGAATTTg                 >  1:39292/3‑137 (MQ=255)
                                                                                                                                              gccgccAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCGGGAATAGAATTTGTTCACAGGCGTCTGc  >  1:1323023/1‑150 (MQ=255)
                                                                                                                                                 gccAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCggg                            >  2:671304/1‑121 (MQ=255)
                                                                                                                                                 gccAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCggg                            <  1:671304/121‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GCCCTTGGCGTGGTTCTTGCGGTAGCCTGGGTAATGCCCGGCCTGGGCTTCGAAGGTGTCGATGATGCGCTGTGGCGTCAGTCGGGGGGCGCCGAGCCAGCCAGCGGCATAG‑GCAAAGCCGGCACCTGCGGCCAGCAGCACGCCGCCAATGGCGGCCAGGCGCAAGGCTTTGGCCGGGCCGTGCAAGGGTGTGTTCATGGTCATCGTCCGGTTCAGTGAAGGTGCCAGTACGACGTGTCAGAGAAAAAGTTATTCCCTGATCGGGAATAGAATTTGTTCACAGGCGTCTGC  >  NC_002947/3285704‑3285994

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