Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,610,902 Δ1 bp 6.5% coding (2825/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,9020T. 6.5% 188.5 / 7.2 62coding (2825/3396 nt)mcoAMn(II) copper oxidase A
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base . (2/2);  total (28/34)
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.48e-01

TGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTC  >  NC_002947/3610754‑3611050
                                                                                                                                                    |                                                                                                                                                     
tgtgCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTTGCCAGCACGTTCAGGGTGATCGGCTTGCGTGCAGTCGTTGCGGCACTGTCGtt                                                                                                                                                      >  2:46251/1‑150 (MQ=255)
     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                                                                                    <  2:1300951/110‑1 (MQ=255)
                                                                                                         gTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTg                                                                                    >  1:1300951/1‑110 (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)
                                                                                                                                                  cGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGtt   <  2:534278/150‑1 (MQ=255)
                                                                                                                                                   gTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTc  <  1:327169/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTGCTGCCGTCCGTGCTCACGGTACCCTGGCCGGAGTCCGGCTGGCTGAGGCCGGTGACCGTGAGCGGTACGTTGCCGTCCGGGTCGCTGTCGTTGGCCAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTC  >  NC_002947/3610754‑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: