Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,383,967 G→C 7.6% V406V (GTC→GTG nasA ← nitrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,383,9670GC 7.6% 202.0 / 6.8 66V406V (GTC→GTGnasAnitrate transporter
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base C (2/3);  total (31/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 = 8.91e-01

CGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGCCCGGCCGCCAGCAG  >  NC_002947/2383819‑2384116
                                                                                                                                                    |                                                                                                                                                     
cGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGa                                                                                                                                                      <  1:90679/150‑1 (MQ=255)
  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:867798/1‑150 (MQ=255)
                                                                                                             ttGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGcc                                         >  1:1420255/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                                                                                                                     >  1:899259/1‑54 (MQ=255)
                                                                                                                                 cTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCt                                                                                                                     <  2:899259/54‑1 (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                                                                                                    <  2:335242/64‑1 (MQ=255)
                                                                                                                                        ccGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAgg                                                                                                    >  1:335242/1‑64 (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)
                                                                                                                                                    |                                                                                                                                                     
CGAACAGGTAGCCTTTGCTGGCGGCCAGTTCTGGTGCCGGCGTGACCAGGCGCAGGGCGTCCTGGTTTTCCGCGCGCGGCCCGGTGGCGCTGGCTTGGGCAAAGCTCAGTTGCAGGCTCATCGGCGTGCCTCAGACCCGTGCTGCGGTGACTGCAGCCGAGCCCCAGGTGGTGCGCCAGCGCTGTTTCACCCCGTGCAGGCCGAACCAGGCCAGCAGGCCCAGGCTGGCGAACAGCCACAGCCCCAGTTGGTAATCGCCAGTGTGCTGCTTGATGGTGCCAAGCCCGGCCGCCAGCAG  >  NC_002947/2383819‑2384116

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: