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)
               aTAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGGGCGA‑GGGGGGGACTTGAAccc                                                                                                                             >  1:349545/1‑151 (MQ=255)
               aTAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCATATTACTTTTATTGATATCCTATCGTGTGGTGCGA‑GGGTGGGACTTTAAccc                                                                                                                             >  2:191404/1‑151 (MQ=17)
                    tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCTTCTTACTTTTTATGATATCCTCCCGTTGGGGGGGAgggggggggttgtcacctc                                                                                                                            >  2:671732/1‑137 (MQ=17)
                    tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGTCTTGAACCCCCACg                                                                                                                        >  1:442201/1‑151 (MQ=255)
                    tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCAc                                                                                                                         <  1:5537/150‑1 (MQ=255)
                     ggTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                        <  4:7604/150‑1 (MQ=255)
                     ggTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACg                                                                                                                        <  2:1275013/150‑1 (MQ=255)
                       tGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTcc                                                                                                                     >  4:283371/1‑151 (MQ=255)
                       tGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCATTAGTACTTTTTTTGATATTGTCGCGTAAGGGGCGA‑GGGGAGGACTTGAACCCCCACGTcc                                                                                                                     >  2:76202/1‑151 (MQ=17)
                        gctaGGGGGGGGGCTTGAACCCCCACGTCCACATCGGTAAGAACACCTGCAGCTGACGCTAGTACCAATTCCGCTACCGTCGCCCAGGCACCGTACGTGTTATGATATCGTCTTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAc                                                                                                                         <  2:425235/144‑1 (MQ=255)
                        gCGAGGGGGGGGACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATGGCCTCGTTTGGTGCGA‑GTGGGGGACTTGTACCCCCATGTCCg                                                                                                                    >  4:97131/1‑151 (MQ=21)
                        gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCg                                                                                                                    >  1:682360/1‑151 (MQ=255)
                          gAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCg                                                                                                                    <  2:445142/150‑1 (MQ=255)
                              ggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGa                                                                                                               >  2:1097751/1‑150 (MQ=255)
                                    aCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                        <  1:473790/151‑1 (MQ=255)
                                    aCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCCCGTTTGGGGCGA‑GGGGGGGACTTGAAACCCCACGTCCGTAAGAAAACTa                                                                                                         >  2:1085812/1‑150 (MQ=25)
                                         aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                                                                   <  1:167582/151‑1 (MQ=255)
                                               cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCt                                                                                             >  2:264409/1‑151 (MQ=255)
                                                 cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                                            <  1:586602/150‑1 (MQ=255)
                                                 cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                           <  1:522566/151‑1 (MQ=255)
                                                 cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                           <  2:811693/151‑1 (MQ=255)
                                                    ccGTAAGGACACTAACACTTGAAACTAGAGCGTATGCCAATTCCGCCAGCTTCGCACAGTAAGCTTACTGTATTTGATATCACCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcgc                                                                                        <  2:251357/151‑1 (MQ=21)
                                                    ccGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                         <  1:1144909/150‑1 (MQ=255)
                                                    ccGTAAGAACACTAAAACCTGATGCTGGCACGTATTCCAATTCCGCCTCGATAGCACAGGCATCTTACTTTTTTTGATATCGCCTTGTTCGGTGCGA‑GGGGGAGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCtcg                                                                                           <  2:577635/148‑3 (MQ=21)
                                                     cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgca                                                                                         <  2:1200441/149‑2 (MQ=255)
                                                     cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGGCATCTTACTTTTTGTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                                            <  2:1087502/146‑1 (MQ=255)
                                                         aGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTa                                                                                   <  1:355009/151‑1 (MQ=255)
                                                             cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                               <  1:589850/151‑1 (MQ=255)
                                                              acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGCTTGGGGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGCCTACCaa                                                                               >  2:1211436/1‑150 (MQ=255)
                                                                        gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACCCTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                                      <  4:339980/149‑1 (MQ=255)
                                                                          aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                                  <  1:324282/151‑1 (MQ=255)
                                                                            cTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                                >  1:545810/1‑151 (MQ=255)
                                                                              agaacgtattcgaaTTCCGCCTGCTTCGCACAGGCAGCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                                                                   <  2:705656/134‑1 (MQ=18)
                                                                                 acgTCTACCAATCCCGCCACCATAGCACAGTCATCTTACATATTTGGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                                                             <  2:256678/148‑1 (MQ=18)
                                                                                  cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                          <  2:409747/151‑1 (MQ=255)
                                                                                   gtatACCAATCCCGCCTCATTCGCACAGTCAGCTTACTTTATTTGATATCGACTCATTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                         <  2:643047/148‑1 (MQ=18)
                                                                                        ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                      <  1:618319/149‑1 (MQ=255)
                                                                                          aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                  >  3:250119/1‑151 (MQ=255)
                                                                                            ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                  >  3:38490/1‑149 (MQ=255)
                                                                                              ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATc                                               >  1:1057959/1‑150 (MQ=255)
                                                                                                gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                            <  1:439985/151‑1 (MQ=255)
                                                                                                gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                            <  1:699732/151‑1 (MQ=255)
                                                                                                gCCACCTTCGCACAGGCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATc                                               <  2:33738/148‑1 (MQ=255)
                                                                                                gCCACCAGCGCACAGTCATCTAACTTAAAGGGAGATCACATCGTTT‑GTGCGGAGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                            <  2:860869/151‑1 (MQ=17)
                                                                                                   aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTg                                         >  1:1211546/1‑151 (MQ=255)
                                                                                                    ccTACGCACAGTCAGCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                          <  2:768717/149‑1 (MQ=255)
                                                                                                     cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGt                                       >  1:690673/1‑151 (MQ=255)
                                                                                                        cGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtt                                      <  2:817516/149‑1 (MQ=255)
                                                                                                          cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttttt                                   <  1:194770/150‑1 (MQ=255)
                                                                                                          cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                                  >  3:350380/1‑151 (MQ=255)
                                                                                                          caaaGTCATCTTACTTTATTGGATATCGCCTCGTTTGGTCCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                                  <  1:767858/148‑1 (MQ=21)
                                                                                                           acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCTATAACAACGTGGtttttcg                                 >  2:839451/1‑149 (MQ=255)
                                                                                                            caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                                 <  1:420105/150‑1 (MQ=255)
                                                                                                              gTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGc                              >  1:410435/1‑151 (MQ=255)
                                                                                                              gTCATCTTACTTTATTTGATAACACCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGc                              <  2:458656/151‑1 (MQ=255)
                                                                                                                 aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                                 >  3:178878/1‑145 (MQ=255)
                                                                                                                 aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                            <  2:850594/150‑1 (MQ=255)
                                                                                                                 aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                            <  2:1160870/150‑1 (MQ=255)
                                                                                                                     tACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggg                         <  1:810632/149‑1 (MQ=255)
                                                                                                                      aCTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggg                      >  2:528972/1‑151 (MQ=255)
                                                                                                                       cTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggggg                     >  4:143233/1‑151 (MQ=255)
                                                                                                                                atatCGGCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATGTCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGaa             <  3:221858/150‑1 (MQ=255)
                                                                                                                                atatCGCCTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGCTTTTGGTGCGAGGCGGGGGACTTGAAc            >  3:225031/1‑151 (MQ=255)
                                                                                                                                 tttCGCCACGTTTCGTGCGA‑GGGGTGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCTCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACt                 <  2:279115/143‑1 (MQ=21)
                                                                                                                                      ccTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCAc      <  2:1001815/151‑1 (MQ=255)
                                                                                                                                       cTCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCAATATCAACGTGGTTCTTCCTCCGAGGCGGGGGCCTTGGACCCCCa       >  2:200241/1‑149 (MQ=18)
                                                                                                                                        tCGTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACCGCGAGCACTATCAACGTGGTTCTTGGTGCGAGGCGCGGGCCAGGAACCCTCACGt    >  3:286540/1‑151 (MQ=18)
                                                                                                                                          gTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc  <  2:795932/151‑1 (MQ=255)
                                                                                                                                          gTTTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTcc  <  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: ‑