Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,842,587 C→T Q73* (CAG→TAG)  gdhA → glutamate dehydrogenase, NADP‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,842,5870CT96.9% 222.6 / ‑1.5 65Q73* (CAG→TAG) gdhAglutamate dehydrogenase, NADP‑specific
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (29/34);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AGTTCGCGCAAGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTG  >  NC_000913/1842438‑1842727
                                                                                                                                                     |                                                                                                                                            
aGTTCGCGCAAGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACCa                                                                                                                                             <  2:962822/151‑1 (MQ=255)
     gcgcAAGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATAc                                                                                                                                        >  2:649941/1‑151 (MQ=255)
      cgcAAGCCGTGCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACa                                                                                                                                       <  1:759140/151‑1 (MQ=255)
        cAAGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAgg                                                                                                                                     <  2:526047/151‑1 (MQ=255)
         aaGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGt                                                                                                                                    >  2:365893/1‑151 (MQ=255)
            ccGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCa                                                                                                                                  >  2:870971/1‑150 (MQ=255)
             cGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAAc                                                                                                                                >  1:494771/1‑151 (MQ=255)
              gctCGGGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGGATGGGTTGATGATCGCAACTAGATACAGGTCAAc                                                                                                                                >  1:520613/3‑150 (MQ=255)
              gTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAAcc                                                                                                                               <  1:299558/151‑1 (MQ=255)
                 cGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTg                                                                                                                            <  1:1022205/151‑1 (MQ=255)
                    gAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCAt                                                                                                                         <  2:646758/151‑1 (MQ=255)
                          aTGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCgtg                                                                                                                   <  1:12398/151‑1 (MQ=255)
                            gACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCgtgt                                                                                                                  <  1:582290/150‑1 (MQ=255)
                                   ctctGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTc                                                                                                           >  1:315072/1‑150 (MQ=255)
                                     ctGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGc                                                                                                        >  1:206969/1‑151 (MQ=255)
                                         ccTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGTTCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGctct                                                                                                     >  1:1174729/1‑150 (MQ=255)
                                             ttCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCa                                                                                                 <  2:481410/150‑1 (MQ=255)
                                                ttGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCgg                                                                                             >  2:702911/1‑151 (MQ=255)
                                                  gAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGc                                                                                            <  1:730685/150‑1 (MQ=255)
                                                     cAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCg                                                                                         >  2:602974/1‑150 (MQ=255)
                                                       aaaTCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTAc                                                                                      <  2:1165310/151‑1 (MQ=255)
                                                        aaTCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACa                                                                                     <  2:183451/151‑1 (MQ=255)
                                                             aaaaTATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAgg                                                                                 >  2:346073/1‑150 (MQ=255)
                                                               aaTATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCgg                                                                              >  1:928312/1‑151 (MQ=255)
                                                                atatCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGt                                                                             <  1:599353/151‑1 (MQ=255)
                                                                     gCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATgcg                                                                        >  1:537893/1‑151 (MQ=255)
                                                                       cAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCCGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGACCCAGGTCAACCGTGCATGGCGAGTGCCGTTCAGCTGTGCCATCGGCCCGTACAACGGCGGTATGCGCt                                                                      <  1:135015/151‑1 (MQ=255)
                                                                       cAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATgcgc                                                                       >  2:473031/1‑150 (MQ=255)
                                                                        aGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCtt                                                                     <  2:625421/151‑1 (MQ=255)
                                                                         gATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGTGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTc                                                                    <  1:1220019/151‑1 (MQ=255)
                                                                          aTGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTcc                                                                   >  1:859789/1‑151 (MQ=255)
                                                                           tGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTcc                                                                   >  2:129725/1‑150 (MQ=255)
                                                                            gTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCa                                                                  <  1:1190516/150‑1 (MQ=255)
                                                                             tCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATc                                                                >  2:116572/1‑151 (MQ=255)
                                                                              cATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATcc                                                               <  2:631049/151‑1 (MQ=255)
                                                                              cATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATcc                                                               >  2:940321/1‑151 (MQ=255)
                                                                               aTTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATcc                                                               <  1:915691/150‑1 (MQ=255)
                                                                               aTTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCg                                                              >  2:249747/1‑151 (MQ=255)
                                                                                 tACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGt                                                             <  1:458342/150‑1 (MQ=255)
                                                                                 tACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTc                                                            >  2:169138/1‑151 (MQ=255)
                                                                                  aCTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTc                                                            <  1:817185/150‑1 (MQ=255)
                                                                                   cTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCa                                                           <  1:350243/150‑1 (MQ=255)
                                                                                     ggggCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCGCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAg                                                          <  1:361248/146‑1 (MQ=255)
                                                                                     ggAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAg                                                          <  1:1185586/149‑1 (MQ=255)
                                                                                          gTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCt                                                   <  1:197972/151‑1 (MQ=255)
                                                                                          gTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCt                                                   <  1:467107/151‑1 (MQ=255)
                                                                                               gTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCttt                                                 <  1:785388/148‑1 (MQ=255)
                                                                                                  gAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCAtt                                            >  2:920148/1‑150 (MQ=255)
                                                                                                  gAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGCATGCGCTTCCATCCGTCAGTTAACCTTTCCAtt                                            >  2:195195/1‑150 (MQ=255)
                                                                                                   aaCCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATtc                                           >  1:972776/1‑150 (MQ=255)
                                                                                                   aaCCGGAGCGCGTGATCCAGTTTCGCGTGGGATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATtct                                          <  2:684856/151‑1 (MQ=255)
                                                                                                     ccGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCa                                        <  2:7178/151‑1 (MQ=255)
                                                                                                     ccGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCa                                        <  1:324475/151‑1 (MQ=255)
                                                                                                      cGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATtctc                                         >  1:108323/1‑149 (MQ=255)
                                                                                                        gAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAAt                                     >  2:244661/1‑151 (MQ=255)
                                                                                                        gAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAAt                                     <  1:387321/151‑1 (MQ=255)
                                                                                                            gcgTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTcc                                  >  1:1162498/1‑150 (MQ=255)
                                                                                                              gTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCg                               >  2:64219/1‑151 (MQ=255)
                                                                                                                gATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGc                             <  2:595811/151‑1 (MQ=255)
                                                                                                                    cAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTg                         <  2:503435/151‑1 (MQ=255)
                                                                                                                    cAGTTTCGCGTGGTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTg                         <  1:362420/151‑1 (MQ=255)
                                                                                                                           gcgTGGTATGGGTTGATGATCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAAc                  >  1:727260/1‑151 (MQ=255)
                                                                                                                               ggTATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTc              <  2:612890/151‑1 (MQ=255)
                                                                                                                                 tATGGGTTGATGATCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCa             >  2:722069/1‑150 (MQ=255)
                                                                                                                                           tgatCGCAACTAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTg  >  2:906847/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
AGTTCGCGCAAGCCGTTCGTGAAGTAATGACCACACTCTGGCCTTTTCTTGAACAAAATCCAAAATATCGCCAGATGTCATTACTGGAGCGTCTGGTTGAACCGGAGCGCGTGATCCAGTTTCGCGTGGTATGGGTTGATGATCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTG  >  NC_000913/1842438‑1842727

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑