Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,606,215 (G)8→9 noncoding (72/87 nt) leuP ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,2071.G85.4% 135.1 / 17.6 48noncoding (80/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (3/4);  new base G (20/21);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

AGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCA  >  NC_000913/4606072‑4606343
                                                                                                                                            |                                                                                                                                          
aGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑gggg                                                                                                                                       <  2:875113/140‑1 (MQ=255)
  aTAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑ggggg‑g                                                                                                                                    <  1:66799/140‑1 (MQ=255)
    aaTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCAACCAAATCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑ggggg‑ggg                                                                                                                                  <  2:543267/140‑1 (MQ=255)
      ttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑ggggg‑gggga                                                                                                                                >  1:1103729/1‑138 (MQ=255)
         gTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑ggggg‑ggggactt                                                                                                                             <  2:1103729/140‑6 (MQ=255)
           gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGACTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGa                                                                                                                            >  1:971551/1‑140 (MQ=255)
            cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACCGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGCGG‑G‑‑GGGGG‑GGGACTTGaa                                                                                                                           >  2:1090950/1‑140 (MQ=255)
            cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGaa                                                                                                                           >  2:1213801/1‑140 (MQ=255)
            cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGa                                                                                                                            <  2:701562/139‑1 (MQ=255)
            cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGa                                                                                                                            >  1:701562/1‑139 (MQ=255)
              aGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑ggggg‑g                                                                                                                                    >  1:460949/1‑128 (MQ=255)
              aGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCCGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑ggggg‑g                                                                                                                                    <  2:460949/128‑1 (MQ=255)
                     ggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AGGGGGGG‑GGGACTTGAACCCCCACg                                                                                                                   <  2:1141653/140‑1 (MQ=255)
                        cTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCg                                                                                                               >  2:128908/1‑140 (MQ=255)
                        cTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCg                                                                                                               >  1:843098/1‑140 (MQ=255)
                         ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACg                                                                                                                   <  1:260599/135‑1 (MQ=255)
                         ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACg                                                                                                                   >  2:260599/1‑135 (MQ=255)
                            aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAg                                                                                                           >  1:1149599/1‑140 (MQ=255)
                                cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAaca                                                                                                       <  2:1149599/140‑1 (MQ=255)
                                 ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGCCTTGAACCCCCACGTCCGTAAGaaatc                                                                                                      >  2:331149/1‑137 (MQ=255)
                                   aCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTa                                                                                                    >  2:890066/1‑140 (MQ=255)
                                               acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCt                                                                                        >  2:1002908/1‑140 (MQ=255)
                                                 acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCCTCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGTCTTGAACCCCCACGACCGTAAGAACACTAACACCTGAAGCTAg                                                                                      >  2:389286/1‑140 (MQ=37)
                                                 acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑GAA‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                      >  2:412449/1‑140 (MQ=255)
                                                 acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                      >  2:57817/1‑140 (MQ=255)
                                                  cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                                                                     >  1:1152833/1‑140 (MQ=255)
                                                  cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCCTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                    <  1:610100/140‑1 (MQ=255)
                                                           gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                            >  2:239736/1‑140 (MQ=255)
                                                                   cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                                                                  <  2:209943/126‑1 (MQ=255)
                                                                   cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                                                                  >  1:209943/1‑126 (MQ=255)
                                                                   cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                    <  2:843098/140‑1 (MQ=255)
                                                                   cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                    <  1:1090950/140‑1 (MQ=255)
                                                                   cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                    <  1:128908/140‑1 (MQ=255)
                                                                   cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                    <  1:1002908/140‑1 (MQ=255)
                                                                   cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTCC‑G‑A‑‑GGGGGCGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                    <  2:1152833/140‑1 (MQ=255)
                                                                          aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                             <  1:412449/140‑1 (MQ=255)
                                                                             aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                          <  1:60305/140‑1 (MQ=255)
                                                                             aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGACTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                          <  2:971551/140‑1 (MQ=255)
                                                                             aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATAGCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                          <  1:1194059/140‑1 (MQ=255)
                                                                                  cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AGGGGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                      >  1:342350/1‑140 (MQ=255)
                                                                                  cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                     >  2:721912/1‑140 (MQ=255)
                                                                                    ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                    >  1:504915/1‑123 (MQ=255)
                                                                                    ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                    <  2:504915/123‑1 (MQ=255)
                                                                                    ccACCTACGAACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                   <  2:91001/140‑1 (MQ=255)
                                                                                     cgccTTCGCACAGTCATCTTACTTTTTTTGATATCGACTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGa                                                  <  2:1164585/138‑1 (MQ=255)
                                                                                      aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                 <  1:57817/140‑1 (MQ=255)
                                                                                       ccTTCGCACAGTCATTTTACTTTTTTTGGAATCCCCCCGTTTTGGGCGG‑G‑G‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                 <  2:1213648/140‑1 (MQ=17)
                                                                                        cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa                                              >  2:454722/1‑140 (MQ=255)
                                                                                                  tCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                       <  1:497996/138‑1 (MQ=255)
                                                                                                  tCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                       >  2:497996/1‑138 (MQ=255)
                                                                                                        tACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttttt                              <  1:454722/140‑1 (MQ=255)
                                                                                                            ttttttGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                           <  1:239736/140‑1 (MQ=255)
                                                                                                              ttttGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTGGCGCGTCTACCAATTCCGCCACCTTCGCCCAGAGCGAGCAATATCAACGTGGTTTTTGGAGc                         >  2:778218/1‑140 (MQ=21)
                                                                                                               tttGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                        <  1:389286/140‑1 (MQ=255)
                                                                                                               tttGATATCGCCTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                        >  1:1192160/1‑140 (MQ=255)
                                                                                                                         ccTCGTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggggac              >  1:629581/1‑138 (MQ=255)
                                                                                                                         catCGTTT‑‑GGTGCGG‑G‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTCCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggggac              <  2:846634/138‑3 (MQ=21)
                                                                                                                             gTTT‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGa         <  1:680310/140‑1 (MQ=255)
                                                                                                                              ttt‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCAGGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGaa        <  1:721912/140‑1 (MQ=21)
                                                                                                                              ttt‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGaa        >  1:530947/1‑140 (MQ=255)
                                                                                                                              ttt‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACACCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggg                 <  2:638153/132‑1 (MQ=255)
                                                                                                                              ttt‑‑GGTGC‑G‑AG‑GGGGG‑GGGACTTGAACACCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggg                 >  1:638153/1‑132 (MQ=255)
                                                                                                                                    gTGC‑G‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAAccccc   >  2:938177/1‑140 (MQ=255)
                                                                                                                                     tGC‑G‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGACTTGAACCCCCa  <  1:938177/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
AGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTT‑‑GGTGC‑G‑A‑‑GGGGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCA  >  NC_000913/4606072‑4606343

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