Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,599,194 C→A 5.3% A189A (GCG→GCT dapA ← dihydrodipicolinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,599,1940CA 5.3% 233.8 / 3.5 75A189A (GCG→GCTdapAdihydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  ref base C (35/36);  new base A (2/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

GGCATCAGACGCTGATTAATAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAG  >  NC_000913/2599048‑2599341
                                                                                                                                                  |                                                                                                                                                   
cgCATCAGACGCTGATTAATAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCg                                                                                                                                                    >  4:88755/2‑148 (MQ=255)
 gCATCAGACGCTGATTAACAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCGTATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGc                                                                                                                                                   <  3:192887/148‑1 (MQ=255)
    tCAGACGCTGATTAATAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGcatca                                                                                                                                               <  3:88755/149‑1 (MQ=255)
       gACGCTGATTAATAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACACCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATcg                                                                                                                                            >  7:72353/1‑149 (MQ=255)
        aCGCTGATTAATAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATcgc                                                                                                                                           <  4:84806/149‑1 (MQ=255)
           cTGATTAATAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATcgccgc                                                                                                                                        >  3:204379/1‑149 (MQ=255)
            tGATTAATAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCt                                                                                                                                       >  1:81275/1‑149 (MQ=255)
             gATTAATAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTc                                                                                                                                      >  6:252150/1‑149 (MQ=255)
              aTTAATAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCt                                                                                                                                       <  7:28927/147‑1 (MQ=255)
               ttAATAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTcag                                                                                                                                    <  4:67578/149‑1 (MQ=255)
                  ataaCGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTcagcag                                                                                                                                 >  6:353399/1‑149 (MQ=255)
                     aCGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAAc                                                                                                                              <  4:19679/149‑1 (MQ=255)
                         gTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACaaaa                                                                                                                          <  7:105425/149‑1 (MQ=255)
                            ccTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAatca                                                                                                                       >  5:306063/1‑149 (MQ=255)
                             cTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAatcat                                                                                                                      <  4:204379/149‑1 (MQ=255)
                              tCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTAGCATCATCGCCGCTCAGCAGAACAAAatcatc                                                                                                                     >  3:233749/1‑149 (MQ=255)
                              tCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTAGCATCATCGCCGCTCAGCAGAACAAAatcatc                                                                                                                     >  3:238700/1‑149 (MQ=255)
                                    aaaTGCCCTTCTGCTGCCAGTTTGCACCTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCACAATTTAATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAAc                                                                                                               <  1:257317/149‑1 (MQ=255)
                                       tGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAg                                                                                                            >  7:142175/1‑149 (MQ=255)
                                              ctgctgCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCtt                                                                                                       <  2:59310/147‑1 (MQ=255)
                                                gctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGAt                                                                                                   >  7:167531/1‑149 (MQ=255)
                                                 ctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATCCCCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATc                                                                                                  <  2:260177/149‑1 (MQ=255)
                                                  tgcCAGTTTGCCCATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCt                                                                                                 >  2:121989/1‑149 (MQ=255)
                                                     cAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTAGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGt                                                                                              <  4:233749/149‑1 (MQ=255)
                                                     cAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTAGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGt                                                                                              <  4:238700/149‑1 (MQ=255)
                                                     cAGTTTGCACATCTGGGCAATAGCACGCGCAGCGACGTTAGCCGTAACGCAAATTACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGt                                                                                              <  6:31606/149‑1 (MQ=255)
                                                     cAGTGTGCACATCTGGGCCATATCACGCGCGGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGGCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGt                                                                                              <  6:310021/149‑1 (MQ=255)
                                                      aGTTTGCACATATGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGtt                                                                                             >  7:284673/1‑149 (MQ=255)
                                                         ttGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTa                                                                                           >  6:24908/1‑148 (MQ=255)
                                                             acctCTGCGCCCTAGCCCGCGCTGCGACGTACGCCCTAACGGAAATAACCCCATGACCGCCCAATAGCATGAAGCCCACCGCGCTCGCATCATCGCCGCTCAGCAGAACAACATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgc                                                                                      <  2:194696/146‑1 (MQ=255)
                                                               aTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGt                                                                                    <  5:24908/149‑1 (MQ=255)
                                                               aTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCAAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGt                                                                                    >  5:72097/1‑149 (MQ=255)
                                                                    ggCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTaa                                                                                 <  5:197322/147‑1 (MQ=255)
                                                                    ggCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTaa                                                                                 <  8:254262/147‑1 (MQ=255)
                                                                             aCGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTg                                                                      <  5:353399/149‑1 (MQ=255)
                                                                                     cGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCttt                                                               <  6:50283/148‑1 (MQ=255)
                                                                                     cGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCtt                                                                <  8:284673/147‑1 (MQ=255)
                                                                                      gACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGa                                                             <  1:121989/149‑1 (MQ=255)
                                                                                      gACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGa                                                             >  6:287723/1‑149 (MQ=255)
                                                                                               cGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGataa                                                    >  6:319375/1‑149 (MQ=255)
                                                                                                gTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCACGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACCCGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGataa                                                    >  7:373361/1‑148 (MQ=255)
                                                                                                 tAACCGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATTTGAAACCACCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGataata                                                  <  5:275922/149‑1 (MQ=255)
                                                                                                           aaCCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACtt                                        >  6:38070/1‑149 (MQ=255)
                                                                                                             ccccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTc                                      <  8:167531/149‑1 (MQ=255)
                                                                                                             ccccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTc                                      <  8:142175/149‑1 (MQ=255)
                                                                                                              cccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCg                                     >  7:238557/1‑149 (MQ=255)
                                                                                                               ccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAACTCATCTGAAACCAGCTCTTTGATGTGGGTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGc                                    <  2:244440/149‑1 (MQ=255)
                                                                                                               ccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGc                                    <  5:319375/149‑1 (MQ=255)
                                                                                                               ccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGc                                    >  3:356352/1‑149 (MQ=255)
                                                                                                               ccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGGTGCCTCTTTGATTCCGTTATTATTTTTTACTTTCGc                                    >  3:188689/1‑149 (MQ=255)
                                                                                                                    cccGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGAc                               >  7:226697/2‑149 (MQ=255)
                                                                                                                    aCCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGAc                               >  7:232296/1‑149 (MQ=255)
                                                                                                                    aCCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGAc                               >  7:226701/1‑149 (MQ=255)
                                                                                                                       gCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGc                            >  6:127522/1‑149 (MQ=255)
                                                                                                                           aaTTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCAc                        >  6:32470/1‑149 (MQ=255)
                                                                                                                           aaTTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCAc                        >  4:365012/1‑149 (MQ=255)
                                                                                                                           aaTTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCAc                        >  5:304729/1‑149 (MQ=255)
                                                                                                                           aaTTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCAc                        >  5:380017/1‑149 (MQ=255)
                                                                                                                               gCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGtt                    >  1:65303/1‑149 (MQ=255)
                                                                                                                                cATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCCGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGttt                   <  5:287723/149‑1 (MQ=255)
                                                                                                                                 aTGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTc                  >  4:56968/1‑149 (MQ=255)
                                                                                                                                 aTGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTc                  <  4:356352/149‑1 (MQ=255)
                                                                                                                                    aaGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCgg               <  8:232296/149‑1 (MQ=255)
                                                                                                                                    aaGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCgg               <  8:226701/149‑1 (MQ=255)
                                                                                                                                    aaGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCgg               <  8:226697/149‑1 (MQ=255)
                                                                                                                                      gTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAACCCAGCTCTTTGATCTGGTTTTCACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCATACGGCCCACCGTTTCCGGGa             <  6:306063/149‑1 (MQ=255)
                                                                                                                                      gTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg              >  4:250735/1‑148 (MQ=255)
                                                                                                                                      gTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg              <  3:250735/148‑1 (MQ=255)
                                                                                                                                      gTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGa             <  8:308719/149‑1 (MQ=255)
                                                                                                                                       tCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAg            >  6:363510/1‑149 (MQ=255)
                                                                                                                                         cAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCa          <  8:72353/149‑1 (MQ=255)
                                                                                                                                               gcTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCg     >  2:92465/1‑148 (MQ=255)
                                                                                                                                                cTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCa   <  2:65303/149‑1 (MQ=255)
                                                                                                                                                 tCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCCCCGTTTCCGGGAGCAGATCGCAg  >  8:82609/1‑149 (MQ=255)
                                                                                                                                                 tCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAg  >  2:374779/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGCATCAGACGCTGATTAATAACGCGTGCCTCGGCAAAATGCCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAG  >  NC_000913/2599048‑2599341

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: