Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,383,970 T→C 7.7% A405A (GCA→GCG nasA ← nitrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,383,9700TC 7.7% 180.6 / 4.7 65A405A (GCA→GCGnasAnitrate transporter
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base C (2/3);  total (30/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.33e-01

AACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGCCCGGCCGCCAGCAGGAA  >  NC_002947/2383821‑2384119
                                                                                                                                                     |                                                                                                                                                     
aaCAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACt                                                                                                                                                       >  2:575813/1‑150 (MQ=255)
    ggTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAg                                                                                                                                                   <  1:18006/150‑1 (MQ=255)
     gTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGc                                                                                                                                                  >  1:109869/1‑150 (MQ=255)
     gTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGc                                                                                                                                                  >  2:328770/1‑150 (MQ=255)
      tAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCAGTCACCGCAGca                                                                                                                                                 >  2:767458/1‑149 (MQ=255)
        gCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGa                                                                                                                                               >  1:52365/1‑150 (MQ=255)
        gCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGa                                                                                                                                               <  2:688834/150‑1 (MQ=255)
         ccTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCAGTCACCGCAGcacgg                                                                                                                                              <  1:1347314/150‑5 (MQ=255)
         ccTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCAGTCACCGCAGcacgg                                                                                                                                              >  2:1347314/1‑146 (MQ=255)
         ccTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGACAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCAGTCACCGCAGcacgg                                                                                                                                              <  2:1146524/150‑5 (MQ=255)
           tttGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCAGTCACCGCAGcactggt                                                                                                                                            <  1:767458/150‑7 (MQ=255)
              gCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCa                                                                                                                                         <  1:941011/150‑1 (MQ=255)
                          ttCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCa                                                                                                                             >  2:490544/1‑150 (MQ=255)
                             tGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAgcg                                                                                                                          >  1:492757/1‑150 (MQ=255)
                               gTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCt                                                                                                                        >  2:10102/1‑150 (MQ=255)
                                  ccGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGAGCGGCCCGCATCCGCCCTCTTGCGCAAAGCTCAGTTGCAGGATCATCGGCGTGCCTCAGCCCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGtt                                                                                                                     <  1:1386564/150‑1 (MQ=255)
                                            cAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTg                                                                                                           >  1:1349518/1‑150 (MQ=255)
                                                      gCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAAc                                                                                                 <  1:575813/150‑1 (MQ=255)
                                                        gTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCa                                                                                               <  2:729435/150‑1 (MQ=255)
                                                          ccTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAgg                                                                                             <  1:598921/150‑1 (MQ=255)
                                                           cTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGc                                                                                            >  1:888105/1‑150 (MQ=255)
                                                           cTGGTTATCCACGCGCGGCACGACAGACCTGGATTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTAAGACCCGTGCTGCGGTGACTGCAACCGAGCACCCGGAGGTGCGCCAACGCTGTTTCACCCCGTGCAGGACGAAAAAGGc                                                                                            >  2:1361776/1‑150 (MQ=255)
                                                                           ggCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCANGGCNGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTg                                                                            <  1:1344550/150‑1 (MQ=255)
                                                                            gCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTgg                                                                           <  1:990227/150‑1 (MQ=255)
                                                                                 gTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAac                                                                      >  2:706800/1‑150 (MQ=255)
                                                                                 gTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAac                                                                      <  2:900155/150‑1 (MQ=255)
                                                                                     cgcTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAacagcc                                                                  <  1:10102/150‑1 (MQ=255)
                                                                                      gcTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAacagcca                                                                 >  2:806456/1‑150 (MQ=255)
                                                                                          gCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAacagccacagc                                                             >  1:1337111/1‑150 (MQ=255)
                                                                                              gggCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCa                                                         >  1:1234487/1‑150 (MQ=255)
                                                                                                  aaaGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGt                                                                                                                      >  2:806843/1‑85 (MQ=255)
                                                                                                  aaaGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGt                                                                                                                      <  1:806843/85‑1 (MQ=255)
                                                                                                  aaaGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTg                                                     <  1:328770/150‑1 (MQ=255)
                                                                                                  aaaGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTg                                                     <  1:46248/150‑1 (MQ=255)
                                                                                                        cAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATc                                               <  2:109869/150‑1 (MQ=255)
                                                                                                           ttGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGcc                                            >  1:1420255/1‑150 (MQ=255)
                                                                                                           ttGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGcc                                            >  1:867798/1‑150 (MQ=255)
                                                                                                           ttGCAGGCTAATCGGCGTGCCTCAGACCAATGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGGGCTGTTTCACGCACTGCAGCCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACCGCCCCAGTTTGTAATCGcc                                            >  1:678861/1‑150 (MQ=255)
                                                                                                            tGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCa                                           <  1:638315/150‑1 (MQ=255)
                                                                                                              cAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAgt                                         <  1:806456/150‑1 (MQ=255)
                                                                                                                 gCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGtgtg                                      <  1:871685/150‑1 (MQ=255)
                                                                                                                     aTCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCTGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTg                                                     >  1:808345/1‑131 (MQ=255)
                                                                                                                     aTCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTg                                                     <  2:808345/131‑1 (MQ=255)
                                                                                                                      tCGGCGTGCCTCAGCCCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGGCAGCGAGGGTGCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTgctgct                                 <  1:1236895/150‑1 (MQ=255)
                                                                                                                          cGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTgg                                                                           <  1:1384413/104‑1 (MQ=255)
                                                                                                                          cGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTgg                                                                           >  2:1384413/1‑104 (MQ=255)
                                                                                                                           gTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATg                            <  2:1337111/150‑1 (MQ=255)
                                                                                                                             gCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGt                          <  1:50237/150‑1 (MQ=255)
                                                                                                                              ccTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTg                         <  1:490544/150‑1 (MQ=255)
                                                                                                                               cTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCt                                                                                                                        <  2:899259/54‑1 (MQ=255)
                                                                                                                               cTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCt                                                                                                                        >  1:899259/1‑54 (MQ=255)
                                                                                                                               cTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGc                        >  2:833420/1‑150 (MQ=255)
                                                                                                                               cTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGc                        >  1:626185/1‑150 (MQ=255)
                                                                                                                                tCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGt                                                   <  2:1274382/122‑1 (MQ=255)
                                                                                                                                tCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGt                                                   >  1:1274382/1‑122 (MQ=255)
                                                                                                                                tCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCCGTGTGCTGCTTGATGGTGcc                       >  2:112172/1‑150 (MQ=255)
                                                                                                                                   gACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAg                    >  1:1392627/1‑150 (MQ=255)
                                                                                                                                    aCCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGc                   <  2:1392627/150‑1 (MQ=255)
                                                                                                                                      ccGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAgg                                                                                                       >  1:335242/1‑64 (MQ=255)
                                                                                                                                      ccGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAgg                                                                                                       <  2:335242/64‑1 (MQ=255)
                                                                                                                                      ccGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGccc                 >  1:905241/1‑150 (MQ=255)
                                                                                                                                       cGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGCCCg                <  1:1107508/150‑1 (MQ=255)
                                                                                                                                            tgcGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGCCCGgccgc           <  2:1064447/150‑1 (MQ=255)
                                                                                                                                              cGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGCCCGGCCGCca         <  2:1234487/150‑1 (MQ=255)
                                                                                                                                               ggTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCACCAGTGTGCTGCTTGATGGTGCCAAGCCCGGCCGCcag        >  2:197810/1‑150 (MQ=255)
                                                                                                                                                  gACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGCCCGGCCGCcagcag     <  1:706800/150‑1 (MQ=255)
                                                                                                                                                     tGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGCCCGgcc             <  2:711495/139‑1 (MQ=255)
                                                                                                                                                     tGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGCCCGgcc             >  1:711495/1‑139 (MQ=255)
                                                                                                                                                     tGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGCCCGGCCGCCCGCAGGaa  >  1:17781/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGCCCGGCCGCCAGCAGGAA  >  NC_002947/2383821‑2384119

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: