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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,3221.G99.2% 503.6 / ‑2.9 132noncoding (80/87 nt)leuQtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (101/30);  total (102/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01
*NC_0009134,606,3222.G80.3% 389.4 / 109.3 132noncoding (80/87 nt)leuQtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (20/6);  new base G (82/24);  total (102/30)
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.51e-01

Only 100 of 161 total aligned reads displayed.

ACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATT  >  NC_000913/4606222‑4606413
                                                                                                     ||                                                                                           
  ccccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gg                                                                                           >  2:7187965/1‑101 (MQ=255)
  ccccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gg                                                                                           <  2:7187109/101‑1 (MQ=255)
  ccccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gg                                                                                           >  1:7187964/1‑101 (MQ=255)
   cccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑ggg                                                                                          <  1:7187413/101‑1 (MQ=255)
   cccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑ggg                                                                                          <  1:7187328/101‑1 (MQ=255)
   cccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑ggg                                                                                          <  2:7187398/101‑1 (MQ=255)
     cACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑ggggg                                                                                        <  1:7187420/101‑1 (MQ=255)
      aCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gggggg                                                                                       >  2:7187967/1‑101 (MQ=255)
      aCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gggggg                                                                                       <  2:7187414/101‑1 (MQ=255)
       cGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑ggggggg                                                                                      <  2:7187517/101‑1 (MQ=255)
        gTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gggggggg                                                                                     <  2:7187436/101‑1 (MQ=255)
         tCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑ggggggggg                                                                                    >  2:7187970/1‑100 (MQ=255)
          ccGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gggggggggg                                                                                   >  1:7187971/1‑99 (MQ=255)
            ttAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑ggggggggggac                                                                                 >  2:7187974/2‑97 (MQ=255)
                  acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gggggg                                                                                       >  2:7187977/1‑89 (MQ=255)
                   cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAAcc                                                                            >  2:7187979/1‑101 (MQ=255)
                    acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAAccc                                                                           >  2:7187982/1‑101 (MQ=255)
                    acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAAcccc                                                                          >  1:7187981/1‑101 (MQ=255)
                      tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGGCTTGAAccccc                                                                         >  1:7187984/1‑101 (MQ=25)
                       aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCAc                                                                       <  2:7187772/101‑1 (MQ=255)
                       aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCAc                                                                       <  1:7187543/101‑1 (MQ=255)
                        acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCAc                                                                       >  1:7187986/1‑101 (MQ=255)
                        acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCAc                                                                       <  1:7187504/101‑1 (MQ=255)
                         cacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGCCTTGAAcccc                                                                          >  1:7187989/1‑97 (MQ=25)
                            cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTcc                                                                   >  2:7187994/1‑101 (MQ=255)
                              gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGt                                                                 >  1:7187995/1‑101 (MQ=255)
                               aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGt                                                                     >  2:7187999/1‑96 (MQ=38)
                               aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTa                                                                >  2:7187996/1‑101 (MQ=255)
                                aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTaa                                                               <  1:7187789/101‑1 (MQ=255)
                                aGCTAGCGCGTCTACCAATTCCGCCACCTTCGACCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTaa                                                               >  2:7188000/1‑101 (MQ=38)
                                 gCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                              >  1:7188001/1‑101 (MQ=255)
                                   tAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                                                                                               >  2:7188005/1‑66 (MQ=255)
                                   tAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGa                                                            >  2:7188003/1‑101 (MQ=255)
                                    aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gggggggggg                                                                                   >  1:7188006/1‑73 (MQ=255)
                                       gcgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTa                                                       >  2:7188008/1‑101 (MQ=255)
                                        cgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTa                                                       >  1:7188012/1‑101 (MQ=255)
                                        cgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTa                                                       >  2:7188011/1‑101 (MQ=255)
                                        cgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                      >  2:7188010/1‑101 (MQ=255)
                                          tCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAac                                                     >  2:7188014/1‑101 (MQ=255)
                                           cTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                      >  2:7188018/1‑99 (MQ=37)
                                           cTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca                                                    >  1:7188016/1‑101 (MQ=255)
                                              ccAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAacac                                                   >  2:7188020/1‑99 (MQ=37)
                                              ccAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCt                                                 >  1:7188019/1‑101 (MQ=255)
                                                aaTTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa                                               >  2:7188026/1‑101 (MQ=25)
                                                 aTTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGaa                                              >  2:7188027/1‑101 (MQ=255)
                                                 aTTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGaa                                              <  1:7187515/101‑1 (MQ=255)
                                                 aTTCCGCCACCTTCACCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGaa                                              <  1:7187463/101‑1 (MQ=25)
                                                     cGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGac                                                           >  2:7188029/1‑84 (MQ=37)
                                                      gCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                         >  2:7188031/1‑101 (MQ=255)
                                                       ccACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc                                        >  2:7188032/1‑101 (MQ=255)
                                                        cACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                       >  1:7188034/1‑101 (MQ=255)
                                                         aCCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                     >  2:7188035/1‑101 (MQ=255)
                                                          ccTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAAc                                                                             >  1:7188040/1‑61 (MQ=38)
                                                          ccTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGCCCGTAAGGacac                                                         >  2:7188038/1‑81 (MQ=25)
                                                          ccTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGt                                    >  2:7188037/1‑101 (MQ=255)
                                                           cTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                   >  2:7188041/1‑101 (MQ=255)
                                                            ttCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc                                        >  2:7188044/1‑96 (MQ=21)
                                                            ttCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                   >  1:7188043/1‑101 (MQ=255)
                                                             tCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCt                                  >  1:7188046/1‑101 (MQ=255)
                                                               gCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAAcccc                                                                          >  2:7188050/1‑59 (MQ=37)
                                                               gCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc                                        >  1:7188049/1‑92 (MQ=255)
                                                                cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAAcccc                                                                          <  1:7187444/58‑1 (MQ=37)
                                                                 ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCa                              >  1:7188053/1‑101 (MQ=255)
                                                                  cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAAccccc                                                                         >  1:7188057/1‑57 (MQ=37)
                                                                  cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACTTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                             >  1:7188055/1‑101 (MQ=18)
                                                                  cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                             >  2:7188054/1‑101 (MQ=255)
                                                                   aGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt                            >  2:7188059/1‑101 (MQ=255)
                                                                     tGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACg                                                                      >  2:7188062/1‑57 (MQ=25)
                                                                     tGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGc                                            >  2:7188061/1‑83 (MQ=18)
                                                                     tGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                         <  1:7187525/101‑1 (MQ=255)
                                                                      gCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                         >  2:7188064/1‑101 (MQ=255)
                                                                      gCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                         >  2:7188066/1‑101 (MQ=255)
                                                                         aGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGc                                            <  1:7187643/79‑1 (MQ=18)
                                                                         aGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                      >  1:7188069/1‑101 (MQ=255)
                                                                         aGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                      >  1:7188068/1‑101 (MQ=255)
                                                                         aGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                     >  2:7188067/1‑101 (MQ=255)
                                                                            aaTATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                         <  1:7187541/95‑1 (MQ=17)
                                                                             atatCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAAc                                                                             <  1:7187458/42‑1 (MQ=25)
                                                                             atatCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                       >  2:7188072/1‑80 (MQ=17)
                                                                              tatCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                 >  2:7188073/1‑101 (MQ=255)
                                                                              tatCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                 >  1:7188074/1‑101 (MQ=255)
                                                                               atCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                  >  1:7188077/1‑99 (MQ=14)
                                                                               atCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                >  1:7188076/1‑101 (MQ=255)
                                                                                  aaCGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAAccccc                                                                         <  2:7187526/41‑1 (MQ=18)
                                                                                  aaCGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt                            <  1:7187836/86‑1 (MQ=14)
                                                                                  aaCGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                <  1:7187475/98‑1 (MQ=14)
                                                                                   aCGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACg                                                                      >  1:7188082/1‑43 (MQ=17)
                                                                                   aCGTGGTTTTTGGTGCGAGGGGGGGGGCCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATa           >  2:7188080/1‑101 (MQ=255)
                                                                                    cGTGGTTTTTGGTGCGAGGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac                                                         >  1:7188084/1‑55 (MQ=17)
                                                                                    cGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATa           >  2:7188083/1‑101 (MQ=255)
                                                                                     gTGGTTTTTGGTGGGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa                                               <  1:7187381/64‑1 (MQ=14)
                                                                                     gTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATAc          >  1:7188087/1‑101 (MQ=255)
                                                                                     gTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATAcc         >  2:7188085/1‑101 (MQ=255)
                                                                                      tGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTaa                                                               <  1:7187600/47‑1 (MQ=14)
                                                                                       ggTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACg                                                                      >  2:7188090/1‑39 (MQ=14)
                                                                                       ggTTTTTGGGGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGaa                                              <  1:7187478/63‑1 (MQ=11)
                                                                                         tttttGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                      <  2:7187466/53‑1 (MQ=11)
                                                                                         tttttGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                           <  2:7187858/79‑1 (MQ=255)
                                                                                         tttttGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCAt       >  1:7188093/1‑99 (MQ=255)
                                                                                         tttttGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCa     >  1:7188091/1‑101 (MQ=255)
                                                                                                     ||                                                                                           
ACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATT  >  NC_000913/4606222‑4606413

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: