Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,606,215 (G)8→7 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,2080G.97.0% 274.3 / 3.5 68noncoding (79/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base . (26/39);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

AGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCC  >  NC_000913/4606059‑4606348
                                                                                                                                                     |                                                                                                                                            
aGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgg                                                                                                                                             <  2:1046444/151‑1 (MQ=255)
    aGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg                                                                                                                                         >  4:381149/1‑151 (MQ=255)
    aGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCACCTTACTTTTTTTGATATCCCATCGTTTGGTCCGAgggtgg                                                                                                                                         >  4:36696/1‑151 (MQ=255)
      gAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggggg                                                                                                                                        >  2:810470/1‑150 (MQ=255)
        aGTAAAGATAATTGGTGCGAGGGGGGGTACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTTCCAATTCCGCCACCTTCGCACAGTCATCTTTCTTTTTTTGATATCGCCTCGTTTGGTGCGAgggtggg                                                                                                                                        <  2:607155/148‑1 (MQ=255)
        aGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCTTATTACTTTTTTTGATTTCACCTCGGTTGGGGCGCgtggggg                                                                                                                                        >  2:1076948/1‑148 (MQ=18)
        aGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggact                                                                                                                                     <  1:378655/151‑4 (MQ=255)
         gTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggggactt                                                                                                                                    <  1:364622/151‑5 (MQ=255)
               aTAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGGGCGAGGGGGGGACTTGAAccc                                                                                                                             >  1:349545/1‑151 (MQ=255)
               aTAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCATATTACTTTTATTGATATCCTATCGTGTGGTGCGAGGGTGGGACTTTAAccc                                                                                                                             >  2:191404/1‑151 (MQ=17)
                    tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCTTCTTACTTTTTATGATATCCTCCCGTTGGGGGGGAgggggggggttgtcacctc                                                                                                                            >  2:671732/1‑137 (MQ=17)
                    tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGTCTTGAACCCCCACg                                                                                                                        >  1:442201/1‑151 (MQ=255)
                    tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCAc                                                                                                                         <  1:5537/150‑1 (MQ=255)
                     ggTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACg                                                                                                                        <  4:7604/150‑1 (MQ=255)
                     ggTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACg                                                                                                                        <  2:1275013/150‑1 (MQ=255)
                       tGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTcc                                                                                                                     >  4:283371/1‑151 (MQ=255)
                       tGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCATTAGTACTTTTTTTGATATTGTCGCGTAAGGGGCGAGGGGAGGACTTGAACCCCCACGTcc                                                                                                                     >  2:76202/1‑151 (MQ=17)
                        gctaGGGGGGGGGCTTGAACCCCCACGTCCACATCGGTAAGAACACCTGCAGCTGACGCTAGTACCAATTCCGCTACCGTCGCCCAGGCACCGTACGTGTTATGATATCGTCTTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAc                                                                                                                         <  2:425235/144‑1 (MQ=255)
                        gCGAGGGGGGGGACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATGGCCTCGTTTGGTGCGAGTGGGGGACTTGTACCCCCATGTCCg                                                                                                                    >  4:97131/1‑151 (MQ=21)
                        gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCg                                                                                                                    >  1:682360/1‑151 (MQ=255)
                          gAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                                                                                                                    <  2:445142/150‑1 (MQ=255)
                              ggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGa                                                                                                               >  2:1097751/1‑150 (MQ=255)
                                    aCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                        <  1:473790/151‑1 (MQ=255)
                                    aCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCCCGTTTGGGGCGAGGGGGGGACTTGAAACCCCACGTCCGTAAGAAAACTa                                                                                                         >  2:1085812/1‑150 (MQ=25)
                                         aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                                                                   <  1:167582/151‑1 (MQ=255)
                                               cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCt                                                                                             >  2:264409/1‑151 (MQ=255)
                                                 cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                                            <  1:586602/150‑1 (MQ=255)
                                                 cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                           <  1:522566/151‑1 (MQ=255)
                                                 cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                           <  2:811693/151‑1 (MQ=255)
                                                    ccGTAAGGACACTAACACTTGAAACTAGAGCGTATGCCAATTCCGCCAGCTTCGCACAGTAAGCTTACTGTATTTGATATCACCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcgc                                                                                        <  2:251357/151‑1 (MQ=21)
                                                    ccGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                         <  1:1144909/150‑1 (MQ=255)
                                                    ccGTAAGAACACTAAAACCTGATGCTGGCACGTATTCCAATTCCGCCTCGATAGCACAGGCATCTTACTTTTTTTGATATCGCCTTGTTCGGTGCGAGGGGGAGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCtcg                                                                                           <  2:577635/148‑3 (MQ=21)
                                                     cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgca                                                                                         <  2:1200441/149‑2 (MQ=255)
                                                     cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGGCATCTTACTTTTTGTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                                            <  2:1087502/146‑1 (MQ=255)
                                                         aGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTa                                                                                   <  1:355009/151‑1 (MQ=255)
                                                             cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                               <  1:589850/151‑1 (MQ=255)
                                                              acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGCTTGGGGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGCCTACCaa                                                                               >  2:1211436/1‑150 (MQ=255)
                                                                        gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACCCTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                                      <  4:339980/149‑1 (MQ=255)
                                                                          aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                                  <  1:324282/151‑1 (MQ=255)
                                                                            cTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                                >  1:545810/1‑151 (MQ=255)
                                                                              agaacgtattcgaaTTCCGCCTGCTTCGCACAGGCAGCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                                                                   <  2:705656/134‑1 (MQ=18)
                                                                                 acgTCTACCAATCCCGCCACCATAGCACAGTCATCTTACATATTTGGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                                                             <  2:256678/148‑1 (MQ=18)
                                                                                  cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                          <  2:409747/151‑1 (MQ=255)
                                                                                   gtatACCAATCCCGCCTCATTCGCACAGTCAGCTTACTTTATTTGATATCGACTCATTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                         <  2:643047/148‑1 (MQ=18)
                                                                                        ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                      <  1:618319/149‑1 (MQ=255)
                                                                                          aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                  >  3:250119/1‑151 (MQ=255)
                                                                                            ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                  >  3:38490/1‑149 (MQ=255)
                                                                                              ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATc                                               >  1:1057959/1‑150 (MQ=255)
                                                                                                gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                            <  1:439985/151‑1 (MQ=255)
                                                                                                gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                            <  1:699732/151‑1 (MQ=255)
                                                                                                gCCACCTTCGCACAGGCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATc                                               <  2:33738/148‑1 (MQ=255)
                                                                                                gCCACCAGCGCACAGTCATCTAACTTAAAGGGAGATCACATCGTTTGTGCGGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                            <  2:860869/151‑1 (MQ=17)
                                                                                                   aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTg                                         >  1:1211546/1‑151 (MQ=255)
                                                                                                    ccTACGCACAGTCAGCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                          <  2:768717/149‑1 (MQ=255)
                                                                                                     cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGt                                       >  1:690673/1‑151 (MQ=255)
                                                                                                        cGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtt                                      <  2:817516/149‑1 (MQ=255)
                                                                                                          cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttttt                                   <  1:194770/150‑1 (MQ=255)
                                                                                                          cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                                  >  3:350380/1‑151 (MQ=255)
                                                                                                          caaaGTCATCTTACTTTATTGGATATCGCCTCGTTTGGTCCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                                  <  1:767858/148‑1 (MQ=21)
                                                                                                           acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCTATAACAACGTGGtttttcg                                 >  2:839451/1‑149 (MQ=255)
                                                                                                            caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                                 <  1:420105/150‑1 (MQ=255)
                                                                                                              gTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGc                              >  1:410435/1‑151 (MQ=255)
                                                                                                              gTCATCTTACTTTATTTGATAACACCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGc                              <  2:458656/151‑1 (MQ=255)
                                                                                                                 aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                                 >  3:178878/1‑145 (MQ=255)
                                                                                                                 aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                            <  2:850594/150‑1 (MQ=255)
                                                                                                                 aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                            <  2:1160870/150‑1 (MQ=255)
                                                                                                                     tACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggg                         <  1:810632/149‑1 (MQ=255)
                                                                                                                      aCTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggg                      >  2:528972/1‑151 (MQ=255)
                                                                                                                       cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggggg                     >  4:143233/1‑151 (MQ=255)
                                                                                                                                atatCGGCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATGTCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGaa             <  3:221858/150‑1 (MQ=255)
                                                                                                                                atatCGCCTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGCTTTTGGTGCGAGGCGGGGGACTTGAAc            >  3:225031/1‑151 (MQ=255)
                                                                                                                                 tttCGCCACGTTTCGTGCGAGGGGTGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCTCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACt                 <  2:279115/143‑1 (MQ=21)
                                                                                                                                      ccTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAc      <  2:1001815/151‑1 (MQ=255)
                                                                                                                                       cTCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCAATATCAACGTGGTTCTTCCTCCGAGGCGGGGGCCTTGGACCCCCa       >  2:200241/1‑149 (MQ=18)
                                                                                                                                        tCGTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACCGCGAGCACTATCAACGTGGTTCTTGGTGCGAGGCGCGGGCCAGGAACCCTCACGt    >  3:286540/1‑151 (MQ=18)
                                                                                                                                          gTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc  <  2:795932/151‑1 (MQ=255)
                                                                                                                                          gTTTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc  <  2:1056195/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
AGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCC  >  NC_000913/4606059‑4606348

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