Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 374,745 A→C 9.1% E12A (GAG→GCG)  PP_0312 → electron transfer flavoprotein subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947374,7450AC 9.1% 214.0 / 7.8 73E12A (GAG→GCG) PP_0312electron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base C (5/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCACAGATCAAGGACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGA  >  NC_002947/374600‑374892
                                                                                                                                                  |                                                                                                                                                   
cccACAGATCAAGGACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTg                                                                                                                                                  <  2:282459/149‑1 (MQ=255)
  cacaGATCAAGGACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGa                                                                                                                                                <  1:88954/149‑1 (MQ=255)
    caGATCAAGGACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATc                                                                                                                                              <  4:219373/149‑1 (MQ=255)
        tCAAGGACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACTCGGCGCGTGGATCGcgg                                                                                                                                          <  1:167571/149‑3 (MQ=255)
           aGGACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGt                                                                                                                                        >  3:17412/1‑148 (MQ=255)
                   gCCGAACTTGTGGCCGACGTGCTGATCGAGGAGGACGCGGCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTg                                                                                                                               <  8:108779/149‑1 (MQ=255)
                   gCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTg                                                                                                                               >  3:185729/1‑149 (MQ=255)
                     cGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACTCGGCGCGTGGATCGcgg                                                                                                                                          <  5:206106/136‑3 (MQ=255)
                     cGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACTCGGCGCGTGGATCGcgg                                                                                                                                          >  6:206106/1‑134 (MQ=255)
                     cGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCa                                                                                                                             >  4:242533/1‑149 (MQ=255)
                          tggtggCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCAcc                                                                                                                           <  3:3000/146‑1 (MQ=255)
                          tggtggCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCAcc                                                                                                                           >  4:3000/1‑146 (MQ=255)
                               gCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCa                                                                                                                    <  4:185729/148‑1 (MQ=255)
                                        cTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATg                                                                                                          >  6:204896/1‑149 (MQ=255)
                                           aTCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCAATGATGAGAGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACACGCTGCATGCCGCGATg                                                                                                          >  2:9080/1‑146 (MQ=255)
                                            tCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGc                                                                                                         >  3:178856/1‑146 (MQ=255)
                                            tCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGc                                                                                                         <  4:178856/146‑1 (MQ=255)
                                                    gaCGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCaa                                                                                                 >  5:2453/1‑146 (MQ=255)
                                                       gCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATc                                                                                                                                              >  7:133637/1‑98 (MQ=255)
                                                       gCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATc                                                                                                                                              <  8:133637/98‑1 (MQ=255)
                                                       gCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGAc                                                                                            >  3:294049/1‑148 (MQ=255)
                                                         aTCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAg                                                                                         >  3:164279/1‑149 (MQ=255)
                                                         aTCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAg                                                                                         >  1:216193/1‑149 (MQ=255)
                                                               tGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCTCCCGCGATCCACGCGCCGCGTGCAACGCCCGTAACCGTATGAACCCCCTCCATGCCGCGATGCAGCCGCAAAAGACCAGATGGa                                                                                    >  1:246506/1‑149 (MQ=255)
                                                                 tCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATg                                                                                  <  4:17412/148‑1 (MQ=255)
                                                                 tCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGGCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATg                                                                                  <  1:111887/148‑1 (MQ=255)
                                                                   ccAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATgggg                                                                               >  1:100375/1‑149 (MQ=255)
                                                                    gaaaGCCGCAACCGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGc                                                                              <  8:278774/148‑1 (MQ=255)
                                                                    cAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGc                                                                              >  8:220612/1‑149 (MQ=255)
                                                                     aaaGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGcc                                                                                                                                           >  7:131568/1‑87 (MQ=255)
                                                                     aaaGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGcc                                                                                                                                           <  8:131568/87‑1 (MQ=255)
                                                                     aaaGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGcc                                                                             >  7:218239/1‑149 (MQ=255)
                                                                        gCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCaa                                                                          >  4:230554/1‑149 (MQ=255)
                                                                         ccGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAAc                                                                         <  3:242533/149‑1 (MQ=255)
                                                                                  gCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCt                                                                                       <  2:277854/126‑1 (MQ=255)
                                                                                  gCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCt                                                                                       >  1:277854/1‑126 (MQ=255)
                                                                                  gCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCatcatc                                                                <  5:204896/149‑1 (MQ=255)
                                                                                   cTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATcc                                                               <  8:218239/149‑1 (MQ=255)
                                                                                    tAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGAGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCg                                                              >  4:145892/1‑149 (MQ=255)
                                                                                       aCGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGc                                                             >  8:115815/1‑147 (MQ=255)
                                                                                         gTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCaa                                                                          <  8:134559/132‑1 (MQ=255)
                                                                                         gTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCaa                                                                          >  7:134559/1‑132 (MQ=255)
                                                                                           gAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAAc                                                       >  7:114368/1‑149 (MQ=255)
                                                                                              ccTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCca                                                    <  3:230554/149‑1 (MQ=255)
                                                                                                 gCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACAt                                                 >  4:7325/1‑149 (MQ=255)
                                                                                                  cGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACAt                                                 <  7:220612/148‑1 (MQ=255)
                                                                                                  cGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACAt                                                 <  2:216193/148‑1 (MQ=255)
                                                                                                   ggAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGc                                               <  6:2453/149‑1 (MQ=255)
                                                                                                    gAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGc                                               >  5:125545/1‑148 (MQ=255)
                                                                                                     aGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGa                                             <  4:294049/149‑1 (MQ=255)
                                                                                                       gTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGAtt                                           >  8:242224/1‑149 (MQ=255)
                                                                                                             tgatgaGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCaa                                                                          >  8:173003/1‑112 (MQ=255)
                                                                                                             tgatgaGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCaa                                                                          <  7:173003/112‑1 (MQ=255)
                                                                                                                   gCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCCCCCTCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGGACCCAAATGCGATTGCCGCAGGCttt                               <  1:9080/149‑1 (MQ=255)
                                                                                                                   gCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCttt                               <  3:184111/149‑1 (MQ=255)
                                                                                                                   gCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCttt                               >  4:184111/1‑149 (MQ=255)
                                                                                                                       cATTATCCGCCGCGCTCCACGCGCCGTGGGGATCGCCCGTACCCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTcg                           <  4:185763/149‑1 (MQ=255)
                                                                                                                           aTCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTcggcc                        <  8:114368/148‑1 (MQ=255)
                                                                                                                                gccgcGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATgggg                                                                               >  6:210007/1‑89 (MQ=255)
                                                                                                                                gccgcGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATgggg                                                                               <  5:210007/89‑1 (MQ=255)
                                                                                                                                gccgcGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTcggcccggcc                   >  7:73405/1‑149 (MQ=255)
                                                                                                                                gccgcGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTcggcccggcc                   <  2:100375/149‑1 (MQ=255)
                                                                                                                                   gcgATCCACTCGGCGCGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGATGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCccggccgg                 >  2:232762/1‑148 (MQ=255)
                                                                                                                                   gcgATCCACTCGGCGCGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCccggccgg                 >  4:20634/1‑148 (MQ=255)
                                                                                                                                   gcgATCCACTCGGCGCGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCccggccgg                 >  8:14987/1‑148 (MQ=255)
                                                                                                                                     gATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCAt              <  7:115815/149‑1 (MQ=255)
                                                                                                                                      aTCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCAt              >  1:155453/1‑148 (MQ=255)
                                                                                                                                         cACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAg          <  2:246506/149‑1 (MQ=255)
                                                                                                                                          aCGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCaacc         >  1:109774/1‑147 (MQ=255)
                                                                                                                                          aCGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAgc         <  3:7325/149‑1 (MQ=255)
                                                                                                                                            gcgcCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAgcgc       <  8:73405/149‑1 (MQ=255)
                                                                                                                                                cGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATc    >  3:27281/1‑148 (MQ=255)
                                                                                                                                                 gAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGACCCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGa  >  3:234848/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCCACAGATCAAGGACTTGGCCGAACTGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGA  >  NC_002947/374600‑374892

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: