Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,496,379 T→C V790A (GTA→GCA)  nirB → nitrite reductase catalytic subunit NirB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,496,3790TC100.0% 1264.3 / NA 306V790A (GTA→GCA) nirBnitrite reductase catalytic subunit NirB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (172/134);  total (172/134)

Only 100 of 306 total aligned reads displayed.

CGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATG  >  NC_000913/3496281‑3496476
                                                                                                  |                                                                                                 
cGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCAc                                                                                                 >  2:5725930/1‑101 (MQ=255)
cGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCAc                                                                                                 >  1:5725927/1‑101 (MQ=255)
cGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCAc                                                                                                 <  1:5725082/101‑1 (MQ=255)
 gAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACt                                                                                                >  2:5725934/1‑101 (MQ=255)
 gAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACt                                                                                                <  2:5725477/101‑1 (MQ=255)
  aaGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACtg                                                                                               >  2:5725938/1‑101 (MQ=255)
  aaGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACtg                                                                                               >  1:5725937/1‑101 (MQ=255)
   aGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACtgt                                                                                              <  2:5725376/101‑1 (MQ=255)
    ggcggcATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGTGCGCCTGCGTGAAGCGGCACtgtg                                                                                             >  2:5725942/1‑101 (MQ=255)
    ggcggcATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACtgtg                                                                                             <  1:5725565/101‑1 (MQ=255)
      cggcATCTATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACtgtgtg                                                                                           >  2:5725947/1‑101 (MQ=255)
      cggcATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACtgtgtg                                                                                           >  1:5725944/1‑101 (MQ=255)
       ggcATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGa                                                                                          <  2:5725536/101‑1 (MQ=255)
          aTCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTg                                                                                       >  1:5725951/1‑101 (MQ=255)
          aTCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTg                                                                                       <  1:5725439/101‑1 (MQ=255)
           tCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTgg                                                                                      <  2:5725557/101‑1 (MQ=255)
            cGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAGGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGa                                                                                     <  1:5725764/101‑1 (MQ=255)
            cGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGa                                                                                     <  2:5725463/101‑1 (MQ=255)
            cGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGa                                                                                     <  1:5725529/101‑1 (MQ=255)
             gATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGAc                                                                                    >  1:5725957/1‑101 (MQ=255)
             gATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGAc                                                                                    >  2:5725954/1‑101 (MQ=255)
               ttACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTg                                                                                  >  2:5725959/1‑101 (MQ=255)
                 aCCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGaa                                                                                >  2:5725962/1‑101 (MQ=255)
                  ccTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGAGTGAAGCGGCACTGTGTGAGTGGACTGaaa                                                                               >  2:5725965/1‑101 (MQ=255)
                   cTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAAc                                                                              <  2:5725649/101‑1 (MQ=255)
                   cTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAAc                                                                              >  1:5725970/1‑101 (MQ=255)
                   cTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAAc                                                                              <  1:5725504/101‑1 (MQ=255)
                      aaaGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGt                                                                           >  1:5725973/1‑101 (MQ=255)
                      aaaGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGt                                                                           <  1:5725396/101‑1 (MQ=255)
                         gCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCaa                                                                        <  2:5725602/101‑1 (MQ=255)
                          cAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAAt                                                                       <  1:5725473/101‑1 (MQ=255)
                          cAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAAt                                                                       >  1:5725977/1‑101 (MQ=255)
                             tGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACg                                                                             >  2:5725983/1‑92 (MQ=255)
                             tGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACg                                                                    >  1:5725982/1‑101 (MQ=255)
                               aTCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGcc                                                                  >  2:5725987/1‑101 (MQ=255)
                               aTCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGcc                                                                  <  1:5725720/101‑1 (MQ=255)
                                tCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCg                                                                 >  1:5725990/1‑101 (MQ=255)
                                  aTTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTc                                                               >  2:5725993/1‑101 (MQ=255)
                                   ttGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCt                                                              <  1:5725561/101‑1 (MQ=255)
                                   ttGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCt                                                              >  1:5725996/1‑101 (MQ=255)
                                     gacgacAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTgc                                                            >  2:5726001/1‑101 (MQ=255)
                                     gacgacAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTgc                                                            >  1:5726002/1‑101 (MQ=255)
                                      acgacAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTgcg                                                           <  2:5725740/101‑1 (MQ=255)
                                       cgacAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTgcgc                                                          >  1:5726006/1‑101 (MQ=255)
                                         acAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAg                                                        >  1:5726010/1‑101 (MQ=255)
                                         acAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAg                                                        <  2:5725424/101‑1 (MQ=255)
                                           aaaCTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGAc                                                      >  1:5726013/1‑101 (MQ=255)
                                           aaaCTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGAc                                                      >  2:5726017/1‑101 (MQ=255)
                                           aaaCTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGAc                                                      >  1:5726021/1‑101 (MQ=255)
                                           aaaCTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGAc                                                      >  1:5726018/1‑101 (MQ=255)
                                              cTGGGGCTGAACGCACCTCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCg                                                   <  2:5725725/101‑1 (MQ=255)
                                              cTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCg                                                   >  2:5726026/1‑101 (MQ=255)
                                               tGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGc                                                  >  2:5726029/1‑101 (MQ=255)
                                               tGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGc                                                  >  1:5726030/1‑101 (MQ=255)
                                                   gCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCa                                              >  1:5726034/1‑101 (MQ=255)
                                                    cTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGt                                                                                        <  1:5725599/58‑1 (MQ=255)
                                                    cTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCaa                                             <  1:5725482/101‑1 (MQ=255)
                                                     tGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCaaa                                            >  2:5726037/1‑101 (MQ=255)
                                                     tGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCaaa                                            <  1:5725512/101‑1 (MQ=255)
                                                      gAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAac                                           <  1:5725360/101‑1 (MQ=255)
                                                       aaCGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAaca                                          >  1:5726040/1‑101 (MQ=255)
                                                        aCGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAacac                                         >  2:5726043/1‑101 (MQ=255)
                                                          gCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTg                                                                                       >  1:5726048/1‑53 (MQ=255)
                                                           cacaTCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtc                                      >  2:5726049/1‑101 (MQ=255)
                                                            acaTCTGGAAGAAGAGATGGCGCGCCTGGGTGAAGCGGCACTGTGTGAGTGGACTGAAACg                                                                             <  1:5725544/61‑1 (MQ=255)
                                                            acaTCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtca                                     <  2:5725575/101‑1 (MQ=255)
                                                            acaTCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtca                                     <  2:5725970/101‑1 (MQ=255)
                                                            acaTCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtca                                     >  2:5726054/1‑101 (MQ=255)
                                                             caTCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtcat                                    >  1:5726056/1‑101 (MQ=255)
                                                              aTCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtcatc                                   >  1:5726059/1‑101 (MQ=255)
                                                                cTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCaa                                 >  2:5726063/1‑101 (MQ=255)
                                                                 tGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtcatcagc                                <  2:5725507/101‑3 (MQ=255)
                                                                  ggAAGAAGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACa                               <  2:5725621/101‑1 (MQ=255)
                                                                   gaagaaGAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAg                              >  2:5726066/1‑101 (MQ=255)
                                                                    aagaagAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGc                             >  1:5726069/1‑101 (MQ=255)
                                                                    aagaagAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGc                             <  2:5725668/101‑1 (MQ=255)
                                                                    aagaagAGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGc                             <  2:5725317/101‑1 (MQ=255)
                                                                     agaagaGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAg                                                                                         >  2:5726070/1‑40 (MQ=255)
                                                                       aagaGATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGAc                          <  2:5725760/101‑1 (MQ=255)
                                                                        agagATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACa                         <  2:5725548/101‑1 (MQ=255)
                                                                        agagATGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACa                         <  1:5725425/101‑1 (MQ=255)
                                                                         gagaTGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACaa                        <  1:5725578/101‑1 (MQ=255)
                                                                          agaTGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAg                       >  1:5726075/1‑101 (MQ=255)
                                                                          agaTGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAg                       <  2:5725630/101‑1 (MQ=255)
                                                                           gaTGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGc                      >  2:5726076/1‑101 (MQ=255)
                                                                             tGGCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGt                    <  1:5725694/101‑1 (MQ=255)
                                                                              ggCGCGCCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTg                   <  1:5725460/101‑1 (MQ=255)
                                                                                  cgcCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGAccc               >  2:5726081/1‑101 (MQ=255)
                                                                                   gcCTGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCg              >  1:5726082/1‑101 (MQ=255)
                                                                                      tGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCa                                              >  2:5726088/1‑66 (MQ=255)
                                                                                      tGCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAAc           >  1:5726087/1‑101 (MQ=255)
                                                                                       gCGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGt                    <  2:5725585/91‑1 (MQ=255)
                                                                                        cGTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGt         <  1:5725521/101‑1 (MQ=255)
                                                                                         gTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGCACGTg        >  1:5726092/1‑101 (MQ=255)
                                                                                         gTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTg        <  2:5725707/101‑1 (MQ=255)
                                                                                         gTGAAGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTg        <  2:5725362/101‑1 (MQ=255)
                                                                                            aaGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAg     >  2:5726097/1‑101 (MQ=255)
                                                                                            aaGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAg     <  1:5725660/101‑1 (MQ=255)
                                                                                            aaGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAg     <  1:5725640/101‑1 (MQ=255)
                                                                                             aGCGGCACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGa    >  1:5726099/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATG  >  NC_000913/3496281‑3496476

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