Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,297,994 T→G D89A (GAT→GCT)  napC ← periplasmic nitrate reductase cytochrome c protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,297,9940TG98.8% 300.5 / ‑4.5 81D89A (GAT→GCT) napCperiplasmic nitrate reductase cytochrome c protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (34/46);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCTGTGCCATCGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAATCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTA  >  NC_000913/2297860‑2298128
                                                                                                                                      |                                                                                                                                      
tCTGTGCCATCGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGa                                                                                                                                   >  1:1162270/1‑140 (MQ=255)
tCTGTGCCATCGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGa                                                                                                                                   >  2:615981/1‑140 (MQ=255)
    tgCCATCGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAgg                                                                                                                               >  1:268070/1‑140 (MQ=255)
      ccATCGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTc                                                                                                                             >  2:44172/1‑140 (MQ=255)
      ccATCGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTc                                                                                                                             >  1:1331697/1‑140 (MQ=255)
        aTCGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCg                                                                                                                            <  2:933188/139‑1 (MQ=255)
        aTCGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCg                                                                                                                            >  1:933188/1‑139 (MQ=255)
         tCGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCa                                                                                                                          >  1:5709/1‑140 (MQ=255)
         tCGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCa                                                                                                                          <  2:1362262/140‑1 (MQ=255)
          cGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCAc                                                                                                                         >  2:666288/1‑140 (MQ=255)
            tCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACgg                                                                                                                       >  2:53921/1‑140 (MQ=255)
             cAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGa                                                                                                                      <  1:615981/140‑1 (MQ=255)
             cAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGa                                                                                                                      <  1:666288/140‑1 (MQ=255)
              aGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGAc                                                                                                                     >  2:785456/1‑140 (MQ=255)
              aGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGAc                                                                                                                     >  1:1076810/1‑140 (MQ=255)
                    tGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCt                                                                                                               <  2:1003655/140‑1 (MQ=255)
                    tGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCt                                                                                                               <  1:688170/140‑1 (MQ=255)
                    tGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCt                                                                                                               <  2:268070/140‑1 (MQ=255)
                    tGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCt                                                                                                               <  2:1331697/140‑1 (MQ=255)
                    tGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCt                                                                                                               <  2:893942/140‑1 (MQ=255)
                    tGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCt                                                                                                               <  2:5709/140‑1 (MQ=255)
                      aGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTAc                                                                                                             >  1:181671/1‑140 (MQ=255)
                      aGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTAc                                                                                                             <  2:452199/140‑1 (MQ=255)
                           cAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGgttg                                                                                                        >  1:938029/1‑140 (MQ=255)
                                 tCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAg                                                                                                  <  2:32624/140‑1 (MQ=255)
                                  cTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGt                                                                                                 >  1:342123/1‑140 (MQ=255)
                                   tGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTg                                                                                                >  2:157185/1‑140 (MQ=255)
                                       gcgTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACg                                                                                            <  2:1074387/140‑1 (MQ=255)
                                       gcgTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACg                                                                                            <  2:901306/140‑1 (MQ=255)
                                          tgtCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGaa                                                                                         >  2:284719/1‑140 (MQ=255)
                                          tgtCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCGTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGaa                                                                                         <  1:478602/140‑1 (MQ=255)
                                           gtCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAAt                                                                                        <  1:620957/140‑1 (MQ=255)
                                           gtCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAAt                                                                                        <  1:1207119/140‑1 (MQ=255)
                                               aTAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCAt                                                                                    <  2:344496/140‑1 (MQ=255)
                                                tAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATg                                                                                   <  1:497020/140‑1 (MQ=255)
                                                tAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATg                                                                                   <  2:342123/140‑1 (MQ=255)
                                                 aaCGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGt                                                                                  <  1:546720/140‑1 (MQ=255)
                                                 aaCGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGt                                                                                  <  1:801201/140‑1 (MQ=255)
                                                 aaCGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGt                                                                                  <  2:455997/140‑1 (MQ=255)
                                                   cGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTAt                                                                                <  2:258651/140‑1 (MQ=255)
                                                   cGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTAt                                                                                <  2:1162270/140‑1 (MQ=255)
                                                   cGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTAt                                                                                <  2:986555/140‑1 (MQ=255)
                                                   cGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTAt                                                                                <  1:1269639/140‑1 (MQ=255)
                                                   cGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTAt                                                                                <  1:553017/140‑1 (MQ=255)
                                                    gCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTAtt                                                                               >  1:1328319/1‑140 (MQ=255)
                                                      cAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGt                                                                                  <  2:738427/135‑1 (MQ=255)
                                                      cAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGt                                                                                  >  1:738427/1‑135 (MQ=255)
                                                      cAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTcc                                                                             >  1:499103/1‑140 (MQ=255)
                                                               tACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACAcc                                                                    <  1:53921/140‑1 (MQ=255)
                                                                       aGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTgcg                                                            <  2:352751/140‑1 (MQ=255)
                                                                                cTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTg                                                   <  1:1226332/140‑1 (MQ=255)
                                                                                cTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTg                                                   <  1:157185/140‑1 (MQ=255)
                                                                                   gCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCa                                                >  1:1091715/1‑140 (MQ=255)
                                                                                   gCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCa                                                <  2:1076810/140‑1 (MQ=255)
                                                                                   gCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCa                                                >  1:1252223/1‑140 (MQ=255)
                                                                                   gCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCa                                                >  2:492252/1‑140 (MQ=255)
                                                                                   gCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCa                                                >  1:688229/1‑140 (MQ=255)
                                                                                     tttGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGc                                              <  2:944405/140‑1 (MQ=255)
                                                                                     tttGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGc                                              <  1:44172/140‑1 (MQ=255)
                                                                                     tttGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGc                                              <  2:688229/140‑1 (MQ=255)
                                                                                      ttGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCt                                             <  2:938029/140‑1 (MQ=255)
                                                                                             tGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAg                                      <  1:1130475/140‑1 (MQ=255)
                                                                                             tGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAg                                      <  1:1287918/140‑1 (MQ=255)
                                                                                              gCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGa                                     >  2:626338/1‑140 (MQ=255)
                                                                                               cGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGaa                                    <  2:181671/140‑1 (MQ=255)
                                                                                               cGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGaa                                    <  1:586231/140‑1 (MQ=255)
                                                                                               cGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGaa                                    >  1:87054/1‑140 (MQ=255)
                                                                                                   tcatcTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAActct                                >  1:1142789/1‑140 (MQ=255)
                                                                                                   tcatcTTCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAActct                                >  2:113493/1‑140 (MQ=255)
                                                                                                        ttCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGt                           >  2:111426/1‑140 (MQ=255)
                                                                                                         tCGGCACAAACTCGTGCGGAACGTGACAATCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTa                          >  2:738310/1‑140 (MQ=255)
                                                                                                         tCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTa                          >  2:663950/1‑140 (MQ=255)
                                                                                                         tCGGCACAAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTa                          >  2:148421/1‑140 (MQ=255)
                                                                                                              acaAACTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTATTGGc                     <  1:492252/140‑1 (MQ=255)
                                                                                                                aaaCTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTa                          >  1:1196654/1‑133 (MQ=255)
                                                                                                                aaaCTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTa                          <  2:1196654/133‑1 (MQ=255)
                                                                                                                 aaCTCGTGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTATTGGCttt                  >  1:1134578/1‑140 (MQ=255)
                                                                                                                       tGCGGAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCAt            >  2:1180168/1‑140 (MQ=255)
                                                                                                                          ggAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATccc         >  2:1361482/1‑140 (MQ=255)
                                                                                                                           gAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCg        <  1:113493/140‑1 (MQ=255)
                                                                                                                           gAACGTGACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCg        <  2:1173165/140‑1 (MQ=255)
                                                                                                                                 gACAAGCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTa  >  2:1014164/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                      
TCTGTGCCATCGTCAGACGATGAGCTTCAAATTTCTGCGGCGTGTCAATAACGCCAAAAATTTTACCATACAGCTCTTTACTTGCTTTGAGCTTGCGTATCATCTTCGGCACAAACTCGTGCGGAACGTGACAATCCGGACAGGTCGCACGGACGCCGCTACGGTTGTTGTAGTGCACGGAATCCATGTATTCCTGATACACCGTGTTGCGCATTTCGTGGCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTA  >  NC_000913/2297860‑2298128

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: