Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,610,907:1 +A 9.7% coding (2820/3396 nt) mcoA ← Mn(II) copper oxidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,610,9071.A 9.7% 218.7 / 19.4 62coding (2820/3396 nt)mcoAMn(II) copper oxidase A
Reads supporting (aligned to +/- strand):  ref base . (24/32);  new base A (3/3);  total (27/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 = 7.56e-01

TGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTC  >  NC_002947/3610759‑3611050
                                                                                                                                                     |                                                                                                                                               
tgcCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCa                                                                                                                                                >  2:772935/1‑150 (MQ=255)
tgcCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAACACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCa                                                                                                                                                >  1:848922/1‑150 (MQ=255)
 gcCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCAc                                                                                                                                               <  1:658953/150‑1 (MQ=255)
  ccgtccgtGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACg                                                                                                                                              >  2:1097366/1‑150 (MQ=255)
      ccgtGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAg                                                                                                                                          <  1:772935/150‑1 (MQ=255)
       ccaacTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAgt                                                                                                                                         >  1:410159/5‑149 (MQ=255)
        gtGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAggt                                                                                                                                        >  1:1079305/1‑150 (MQ=255)
         agCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAgtgc                                                                                                                                       <  2:410159/149‑4 (MQ=255)
           cTCACGGTACCCTGGCCGGCGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGGTGCCGTCCGTGTCGCTGTCGTTGGCCAGCCCGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTGGTTGGCCACGACAggtggt                                                                                                                                     >  2:1266056/1‑150 (MQ=255)
           cTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACTACAggtggt                                                                                                                                     >  1:88666/1‑150 (MQ=255)
            tCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTGGCCAACGACAg                                                                                                                                          >  2:95065/1‑144 (MQ=255)
            tCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTGGCCAACGACAg                                                                                                                                          <  1:95065/144‑1 (MQ=255)
                   aCCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGc                                                                                                                             <  2:88666/150‑1 (MQ=255)
                    cccTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGcc                                                                                                                            <  2:1047915/150‑1 (MQ=255)
                                  cGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGTTCGGCTTGCCTGCAGTCGTTGCGGGACTGTCGTTTTCCACGACAGGTGGTTTGTTGGCCGGgtcgtcatcgtc                                                                                                              <  1:1097366/150‑1 (MQ=255)
                                   ggctggctGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGgtcgtcatcgtcg                                                                                                             <  1:1260175/150‑1 (MQ=255)
                                          tgAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTGGCCAACGACAgt                                                                                                                                         >  1:116111/1‑114 (MQ=255)
                                          tgAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTGGCCAACGACAgt                                                                                                                                         <  2:116111/115‑2 (MQ=39)
                                              gCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGtttt                                                                                                  >  2:893291/1‑150 (MQ=255)
                                               ccGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTg                                                                                                 <  1:625405/150‑1 (MQ=255)
                                                   tGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGc                                                                                             <  2:484849/150‑1 (MQ=255)
                                                      ccGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGa                                                                                          >  2:775111/1‑150 (MQ=255)
                                                       cGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGgtcgtcatcgtcg                                                                                                             >  1:753392/1‑130 (MQ=255)
                                                       cGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGgtcgtcatcgtcg                                                                                                             <  2:753392/130‑1 (MQ=255)
                                                                tACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAg                                                                                <  1:1191786/150‑1 (MQ=255)
                                                                    ttGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCt                                                                            >  2:294174/1‑150 (MQ=255)
                                                                          tCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAg                                                                      <  2:461836/150‑1 (MQ=255)
                                                                           ccGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAgg                                                                     <  1:211010/150‑1 (MQ=255)
                                                                           ccGGGTCGCTGTCGTTGGCCAACACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAgg                                                                     <  2:848922/150‑1 (MQ=255)
                                                                                  gCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTa                                                              >  2:1336450/1‑150 (MQ=255)
                                                                                    tGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAac                                                            >  1:263560/1‑150 (MQ=255)
                                                                                        gTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTTGTGCGGAACGTTGGCAGGc                                                                              <  1:842045/128‑1 (MQ=255)
                                                                                        gTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTTGTGCGGAACGTTGGCAGGc                                                                              >  2:842045/1‑128 (MQ=255)
                                                                                              cAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGccatc                                                  <  1:893291/150‑1 (MQ=255)
                                                                                                  aCGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGccatcccat                                              >  2:327169/1‑150 (MQ=255)
                                                                                                   cGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGccatcccatc                                             >  1:437792/1‑150 (MQ=255)
                                                                                                    gTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTg                                                                                    >  1:1300951/1‑110 (MQ=255)
                                                                                                    gTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTg                                                                                    <  2:1300951/110‑1 (MQ=255)
                                                                                                     ttCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCt                                           <  1:775111/150‑1 (MQ=255)
                                                                                                     ttCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCt                                           <  1:294174/150‑1 (MQ=255)
                                                                                                     ttCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCt                                           >  1:354249/1‑150 (MQ=255)
                                                                                                       cAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTGGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTgg                                         <  1:1266056/150‑1 (MQ=255)
                                                                                                          ggTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAg                                      >  2:124917/1‑150 (MQ=255)
                                                                                                              aTCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCggg                                  >  1:1250085/1‑150 (MQ=255)
                                                                                                              aTCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCggg                                  <  2:437792/150‑1 (MQ=255)
                                                                                                              aTCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGCGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCggg                                  <  1:46251/150‑1 (MQ=255)
                                                                                                                cGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTTATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGt                                >  1:885014/1‑150 (MQ=255)
                                                                                                                cGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCTGGGt                                <  2:354249/150‑1 (MQ=255)
                                                                                                                cGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGt                                <  1:1336450/150‑1 (MQ=255)
                                                                                                                     tGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCa                           <  1:347480/150‑1 (MQ=255)
                                                                                                                        cTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTa                        >  2:1334907/1‑150 (MQ=255)
                                                                                                                         tGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTAc                       >  1:118663/1‑150 (MQ=255)
                                                                                                                           cAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCt                     <  2:1250085/150‑1 (MQ=255)
                                                                                                                               cGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGaa                 <  2:1079305/150‑1 (MQ=255)
                                                                                                                                       cACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTgg                                                                                   <  1:1220833/76‑1 (MQ=255)
                                                                                                                                       cACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTgg                                                                                   >  2:1220833/1‑76 (MQ=255)
                                                                                                                                       accTGTCGTGGCCAACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCAt                                                                 >  2:427444/3‑94 (MQ=255)
                                                                                                                                       accTGTCGTGGCCAACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCAt                                                                 <  1:427444/92‑1 (MQ=255)
                                                                                                                                        aCTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGACGTCATCGTCGTCATCGGTTTTAGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCAGGGTGGGCATTACCTGGAACTGcccgag        >  2:488469/1‑147 (MQ=255)
                                                                                                                                          tGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACTGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGgt      <  1:488469/150‑1 (MQ=255)
                                                                                                                                             tcgtGGCCAACGACAGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGc                                                                                             <  2:428049/57‑1 (MQ=37)
                                                                                                                                             tcgtGGCCAACGACAGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGCTTTGGTGc                                                                                             >  1:428049/4‑60 (MQ=25)
                                                                                                                                             cGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGtt   <  2:534278/150‑1 (MQ=255)
                                                                                                                                              gTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTc  <  1:327169/150‑1 (MQ=255)
                                                                                                                                                    cACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTa                        >  1:1378347/1‑122 (MQ=255)
                                                                                                                                                    cACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTa                        <  2:1378347/122‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTC  >  NC_002947/3610759‑3611050

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