Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,534,038 T→C T202A (ACC→GCC)  narV ← cryptic nitrate reductase 2 gamma subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,534,0380TC87.2% 623.2 / 67.2 251T202A (ACC→GCC) narVcryptic nitrate reductase 2 gamma subunit
Reads supporting (aligned to +/- strand):  ref base T (13/19);  new base C (101/118);  total (114/137)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

GGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGATCCACCA  >  NC_000913/1533948‑1534129
                                                                                          |                                                                                           
ggCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGa                                                                                           <  2:6007963/93‑1 (MQ=255)
ggCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGa                                                                                           >  1:1757251/1‑93 (MQ=255)
 gCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGaa                                                                                          <  2:7291313/93‑1 (MQ=255)
   cAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACg                                                                                        >  1:949458/1‑93 (MQ=255)
   cAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACg                                                                                        >  2:5194531/1‑93 (MQ=255)
    aGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACgg                                                                                       <  1:3691051/93‑1 (MQ=255)
    aGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACgg                                                                                       >  1:5207274/1‑93 (MQ=255)
    aGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACgg                                                                                       >  2:4911965/1‑93 (MQ=255)
       gTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGaa                                                                                    >  2:2909137/1‑93 (MQ=255)
        tGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAag                                                                                   >  1:1534349/1‑93 (MQ=255)
        tGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAag                                                                                   <  1:2324755/93‑1 (MQ=255)
        tGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAag                                                                                   >  2:2171101/1‑93 (MQ=255)
        tGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAag                                                                                   <  2:7418783/93‑1 (MQ=255)
          aaaTTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTTGGGCGAACGGGAagag                                                                                 >  1:1163140/1‑93 (MQ=255)
            aTTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCa                                                                               <  2:5743491/93‑1 (MQ=255)
             ttATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAg                                                                              <  1:6893996/93‑1 (MQ=255)
               aTCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGa                                                                            <  2:6433059/93‑1 (MQ=255)
                tCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGaa                                                                           >  1:1007964/1‑93 (MQ=255)
                 cgccgcGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAg                                                                          <  2:6279139/93‑1 (MQ=255)
                  gccgcGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGa                                                                         >  2:571998/1‑93 (MQ=255)
                   ccgcGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGAt                                                                        >  1:1121546/1‑93 (MQ=255)
                    cgcgAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATg                                                                       <  1:4282611/93‑1 (MQ=255)
                    cgcgAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATg                                                                       >  1:5043071/1‑93 (MQ=255)
                        aaCGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAATGGGAAGAGCAGGAAGATGGTCa                                                                   <  2:4446372/93‑1 (MQ=255)
                        aaCGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCa                                                                   >  2:1418909/1‑93 (MQ=255)
                         aCGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCAt                                                                  <  2:4766324/93‑1 (MQ=255)
                         aCGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCAt                                                                  >  1:249710/1‑93 (MQ=255)
                           gCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTc                                                                >  2:5883213/1‑93 (MQ=255)
                           gCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTc                                                                >  2:3473726/1‑93 (MQ=255)
                           gCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTc                                                                <  1:1727806/93‑1 (MQ=255)
                           gCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTc                                                                >  2:2544465/1‑93 (MQ=255)
                              caATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCa                                                             <  1:2568265/93‑1 (MQ=255)
                               aaTTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCaa                                                            >  2:2926723/1‑93 (MQ=255)
                                aTTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAac                                                           <  2:5859565/93‑1 (MQ=255)
                                aTTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAac                                                           >  2:2649356/1‑93 (MQ=255)
                                   tGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACAcc                                                        <  1:2850617/93‑1 (MQ=255)
                                    gATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCa                                                       <  1:2027449/93‑1 (MQ=255)
                                    gATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCa                                                       >  2:4026740/1‑93 (MQ=255)
                                     atatCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAg                                                      >  1:7138082/1‑93 (MQ=255)
                                      tatCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGa                                                     <  2:380267/93‑1 (MQ=255)
                                       atCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGAt                                                    <  1:2111293/93‑1 (MQ=255)
                                       atCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGAt                                                    <  1:299632/93‑1 (MQ=255)
                                       atCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGAt                                                    <  1:4888221/93‑1 (MQ=255)
                                          gacgaGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCa                                                 <  1:569135/93‑1 (MQ=255)
                                          gacgaGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCa                                                 <  2:1508091/93‑1 (MQ=255)
                                          gacgaGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCa                                                 >  1:3973466/1‑93 (MQ=255)
                                           acgaGTAAAGTATTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAg                                                <  2:1283060/93‑1 (MQ=255)
                                           acgaGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAg                                                <  1:4187132/93‑1 (MQ=255)
                                           acgaGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAg                                                <  1:4079833/93‑1 (MQ=255)
                                           acgaGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAg                                                <  2:5367503/93‑1 (MQ=255)
                                            cgaGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAgg                                               >  1:4728397/1‑93 (MQ=255)
                                             gaGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGc                                              <  2:2775488/93‑1 (MQ=255)
                                             gaGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGc                                              <  2:2223036/93‑1 (MQ=255)
                                              aGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCg                                             >  2:3329089/1‑93 (MQ=255)
                                                tAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGa                                           >  1:5832655/1‑93 (MQ=255)
                                                 aaaGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGaa                                          <  1:3416011/93‑1 (MQ=255)
                                                  aaGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAg                                         <  1:1339989/93‑1 (MQ=255)
                                                  aaGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAg                                         <  2:5442644/93‑1 (MQ=255)
                                                   aGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGa                                        >  1:1469753/1‑93 (MQ=255)
                                                    gTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGAc                                       >  1:5481494/1‑93 (MQ=255)
                                                    gTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGAc                                       >  2:1877703/1‑93 (MQ=255)
                                                     tACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACg                                      >  1:658461/1‑93 (MQ=255)
                                                     tACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACg                                      <  2:5174918/93‑1 (MQ=255)
                                                       cTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGATGaa                                    <  2:6506577/93‑1 (MQ=255)
                                                       cTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGaa                                    <  1:6274680/93‑1 (MQ=255)
                                                       cTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGaa                                    <  2:4698947/93‑1 (MQ=255)
                                                       cTCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGaa                                    >  1:4546900/1‑93 (MQ=255)
                                                        tCAAACGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAAc                                   <  2:4330217/93‑1 (MQ=255)
                                                        tCAAACGGCACGCTCCACACATGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAAc                                   <  2:6893363/93‑1 (MQ=255)
                                                          aaaCGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGc                                 <  2:7201390/93‑1 (MQ=255)
                                                          aaaCGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGc                                 <  1:4722302/93‑1 (MQ=255)
                                                            aCGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTa                               <  2:804152/93‑1 (MQ=255)
                                                            aCGGCGCGCTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTa                               >  2:6777607/1‑93 (MQ=255)
                                                               gcgcgcTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGc                            <  2:2290712/93‑1 (MQ=255)
                                                                 gcgcTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCg                          >  1:4396804/1‑93 (MQ=255)
                                                                  cgcTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGt                         <  1:7471205/93‑1 (MQ=255)
                                                                  cgcTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGt                         <  2:4193483/93‑1 (MQ=255)
                                                                   gcTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGtt                        <  1:3319962/93‑1 (MQ=255)
                                                                    cTCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTg                       >  1:363899/1‑93 (MQ=255)
                                                                     tCCACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGa                      >  1:7494130/1‑93 (MQ=255)
                                                                      ccACACATGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAg                     >  1:7347289/1‑93 (MQ=255)
                                                                       cacgcaTGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGc                    >  1:6975566/1‑93 (MQ=255)
                                                                       cacacaTGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGc                    <  2:5235878/93‑1 (MQ=255)
                                                                        acacaTGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCa                   >  2:3869556/1‑93 (MQ=255)
                                                                         cacaTGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCAt                  >  1:3205743/1‑93 (MQ=255)
                                                                         cacaTGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCAt                  >  2:1987376/1‑93 (MQ=255)
                                                                          acaTGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCAtc                 >  2:2385787/1‑93 (MQ=255)
                                                                          acaTGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCAtc                 >  1:2160249/1‑93 (MQ=255)
                                                                           caTGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCAtct                >  1:6107143/1‑93 (MQ=255)
                                                                           caTGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCAtct                >  2:3166875/1‑93 (MQ=255)
                                                                            aTGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCAtctc               >  2:1184097/1‑93 (MQ=255)
                                                                            aTGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCAtctc               <  2:5673091/93‑1 (MQ=255)
                                                                             tGCACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCATctct              >  2:1760142/1‑93 (MQ=255)
                                                                               cACCAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGa            <  1:6778784/93‑1 (MQ=255)
                                                                                  cAGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGa         >  2:3386283/1‑93 (MQ=255)
                                                                                   aGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGAt        <  1:2652162/93‑1 (MQ=255)
                                                                                   aGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGAt        <  1:3082904/93‑1 (MQ=255)
                                                                                   aGTCGGGCGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGAt        >  1:1872243/1‑93 (MQ=255)
                                                                                    gTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGATc       <  2:5026911/93‑1 (MQ=255)
                                                                                      cGGGTGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGATcca     <  2:744952/93‑1 (MQ=255)
                                                                                          |                                                                                           
GGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGTGAACGGGAAGAGCAGGAAGATGGTCATTCCCAACACCAGATGCAGGCGGAAGACGAACGCTACGCCGTTGAGCATCTCTGACGATCCACCA  >  NC_000913/1533948‑1534129

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