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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,0870G.98.4% 233.6 / ‑2.3 61noncoding (79/87 nt)leuVtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (26/34);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

CAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG  >  NC_000913/4605950‑4606220
                                                                                                                                         |                                                                                                                                     
caAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgg                                                                                                                                      >  2:493390/1‑139 (MQ=255)
  aaCAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggg                                                                                                                                    >  2:1030101/1‑139 (MQ=255)
  aaCAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggg                                                                                                                                    >  2:1028045/1‑139 (MQ=255)
     aGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGAAAAGATAATTGGTGCGAggggggg                                                                                                                                 >  2:919383/1‑139 (MQ=255)
       ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggac                                                                                                                               >  2:929869/1‑137 (MQ=255)
       ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggac                                                                                                                               >  2:1202213/1‑137 (MQ=255)
       ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggac                                                                                                                               >  2:1152260/1‑137 (MQ=255)
       ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggac                                                                                                                               >  1:958303/1‑137 (MQ=255)
        tGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggact                                                                                                                              >  1:537211/1‑136 (MQ=255)
           ggCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGa                                                                                                                          <  1:1028045/139‑1 (MQ=255)
             cATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAAc                                                                                                                        >  1:1198009/1‑139 (MQ=255)
                tACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAAcccc                                                                                                                     <  2:1564728/139‑1 (MQ=255)
                 aCATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAAccccc                                                                                                                    <  1:1202213/139‑1 (MQ=255)
                       tGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTcc                                                                                                              >  2:427505/1‑139 (MQ=255)
                          tACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                           <  2:270536/139‑1 (MQ=255)
                             cATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAgg                                                                                                        <  2:279731/139‑1 (MQ=255)
                                       ggCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAacac                                                                                              >  2:1642321/1‑139 (MQ=255)
                                                tGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                                                                                     <  2:1355791/139‑1 (MQ=255)
                                                tGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATACACGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCCA‑GGGGGGGACTTGAGCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                                                                                     <  1:884408/139‑1 (MQ=37)
                                                 gCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                                    <  2:1645633/139‑1 (MQ=255)
                                                  cTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAgg                                                                                                        >  2:13359/1‑118 (MQ=255)
                                                  cTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAgg                                                                                                        <  1:13359/118‑1 (MQ=255)
                                                   tGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                                                                  <  1:314893/139‑1 (MQ=255)
                                                     gCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                                <  1:995759/139‑1 (MQ=255)
                                                          atGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggactt                                                                                                                             >  1:356840/1‑86 (MQ=255)
                                                          atGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggactt                                                                                                                             <  2:356840/90‑5 (MQ=255)
                                                            gCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCa                                                                         <  1:929869/139‑1 (MQ=255)
                                                            gCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCa                                                                         <  2:416344/139‑1 (MQ=255)
                                                            gCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCa                                                                         <  1:919383/139‑1 (MQ=255)
                                                                     aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                                >  1:156027/1‑139 (MQ=255)
                                                                      aTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                   >  1:942208/1‑103 (MQ=255)
                                                                      aTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                   <  2:942208/103‑1 (MQ=255)
                                                                      aTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                               >  2:743334/1‑139 (MQ=255)
                                                                       tGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                              >  2:1033180/1‑139 (MQ=255)
                                                                          aTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCa                                                                                                                   <  1:575325/83‑1 (MQ=255)
                                                                          aTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCa                                                                                                                   >  2:575325/1‑83 (MQ=255)
                                                                          aTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                           >  1:1579441/1‑139 (MQ=255)
                                                                                   aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGAGGCTCATACCTGAAGCTAGCGCGTCTACCGATTCCGCCACCTTCGcacactgc                                                  >  1:136715/1‑135 (MQ=21)
                                                                                   aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCa                                                  >  1:1129069/1‑139 (MQ=255)
                                                                                     ggATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATc                                                >  1:520177/1‑139 (MQ=255)
                                                                                      gATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCt                                               <  2:15207/139‑1 (MQ=255)
                                                                                      gATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCt                                               <  1:1152260/139‑1 (MQ=255)
                                                                                       aTGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCtt                                              <  2:1616977/139‑1 (MQ=255)
                                                                                         gAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTa                                                                                                  >  1:997356/1‑85 (MQ=255)
                                                                                         gAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTa                                                                                                  <  2:997356/85‑1 (MQ=255)
                                                                                         gAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAg                                                     >  2:383928/1‑130 (MQ=255)
                                                                                         gAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAg                                                     <  1:383928/130‑1 (MQ=255)
                                                                                               tCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACtttttt                                      >  1:232297/1‑139 (MQ=255)
                                                                                                        aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                   <  1:1190602/69‑1 (MQ=255)
                                                                                                        aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                   >  2:1190602/1‑69 (MQ=255)
                                                                                                        aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                             <  1:1642321/139‑1 (MQ=255)
                                                                                                        aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                             <  2:1579441/139‑1 (MQ=255)
                                                                                                        aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                             <  1:1033180/139‑1 (MQ=255)
                                                                                                        aaaCAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                             <  1:493390/139‑1 (MQ=255)
                                                                                                             aGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGt                        <  2:1222921/139‑1 (MQ=255)
                                                                                                                 aGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTcc                                                                                                              <  1:1024440/49‑1 (MQ=21)
                                                                                                                 aGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTcc                                                                                                              >  2:1024440/1‑49 (MQ=21)
                                                                                                                  ggAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGat                                  <  2:1256997/124‑1 (MQ=255)
                                                                                                                  ggAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGat                                  >  1:1256997/1‑124 (MQ=255)
                                                                                                                   gAAGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTg                  >  1:1503216/1‑139 (MQ=255)
                                                                                                                     aGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCg                >  2:1136635/1‑139 (MQ=255)
                                                                                                                     aGTAAAGATAATTGGGGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACAATAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCg                <  2:537211/139‑1 (MQ=255)
                                                                                                                      gTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa               >  1:845121/1‑139 (MQ=255)
                                                                                                                      gTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa               <  2:845121/139‑1 (MQ=255)
                                                                                                                      gTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa               <  2:1503216/139‑1 (MQ=255)
                                                                                                                           gATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCAGCTTACTTTTATTGATATTGCCTCGTTGGGTGCGAggggg          >  1:192370/1‑139 (MQ=18)
                                                                                                                           gATAATTGGGGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggg          <  2:958303/139‑1 (MQ=255)
                                                                                                                             tAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg        >  1:293378/1‑139 (MQ=255)
                                                                                                                              aaTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                          >  1:770161/1‑88 (MQ=12)
                                                                                                                              aaTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                          <  2:770161/88‑1 (MQ=12)
                                                                                                                              aaTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGttt                      >  1:948800/1‑124 (MQ=255)
                                                                                                                              aaTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGttt                      <  2:948800/124‑1 (MQ=255)
                                                                                                                                   gTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTg  >  2:1046270/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                     
CAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTG  >  NC_000913/4605950‑4606220

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