Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,458,416 T→A 5.4% I26I (ATT→ATA gspE → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,458,4160TA 5.4% 526.1 / 15.2 148I26I (ATT→ATAgspEgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base T (71/69);  new base A (3/5);  total (74/74)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 159 total aligned reads displayed.

TTGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAG  >  NC_000913/3458270‑3458563
                                                                                                                                                  |                                                                                                                                                   
ttGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGatta                                                                                                                                                    <  4:673146/148‑1 (MQ=255)
ttGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGatta                                                                                                                                                    <  3:964949/148‑1 (MQ=255)
    cGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGACATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGAttatttat                                                                                                                                               <  1:453954/149‑1 (MQ=255)
      gTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGTGATTATTTAtc                                                                                                                                              >  3:446496/1‑148 (MQ=255)
          tGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATAAATAATCTCGcc                                                                                                                                         <  2:970266/149‑1 (MQ=255)
          tGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATAAATAATCTCGcc                                                                                                                                         <  2:387440/149‑1 (MQ=255)
             aTGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGAttattt                                                                                                                                                 >  2:23095/1‑138 (MQ=255)
             aTGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGa                                                                                                                                       >  4:643831/1‑148 (MQ=255)
              tGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGAt                                                                                                                                      <  4:567096/148‑1 (MQ=255)
                 aaaCACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGc                                                                                                                                  >  1:655678/1‑149 (MQ=255)
                  aaCACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGcc                                                                                                                                 >  1:5024/1‑149 (MQ=255)
                  aaCACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGc                                                                                                                                  >  3:925676/1‑148 (MQ=255)
                   acacGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGcc                                                                                                                                 >  3:130683/1‑148 (MQ=255)
                     acGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATAAATAATCTCGc                                                                                                                                          <  2:786314/137‑1 (MQ=255)
                     acGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATAAATAATCTCGc                                                                                                                                          >  1:786314/1‑137 (MQ=255)
                          cAACAGTCACGCCCCTGCGTCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAg                                                                                                                          <  3:268011/148‑1 (MQ=255)
                           aaCAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCg                                                                                                                        >  2:277893/1‑149 (MQ=255)
                           aaCAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCg                                                                                                                        >  2:1016029/1‑149 (MQ=255)
                             cAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCACCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCgg                                                                                                                       >  3:701290/1‑148 (MQ=255)
                                  aCGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTc                                                                                                                 >  1:548545/1‑149 (MQ=255)
                                    gCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGa                                                                                                               <  1:226719/149‑1 (MQ=255)
                                     ccccTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGa                                                                                                               <  3:667061/148‑1 (MQ=255)
                                     ccccTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGa                                                                                                               >  4:348878/1‑148 (MQ=255)
                                      cccTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGata                                                                                                             <  1:6947/149‑1 (MQ=255)
                                               aTCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCgg                                                                                                                       <  4:162339/130‑1 (MQ=255)
                                               aTCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCgg                                                                                                                       >  3:162339/1‑130 (MQ=255)
                                               aTCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGt                                                                                                     >  3:402249/1‑148 (MQ=255)
                                                  gTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAgg                                                                                                 >  2:172293/1‑149 (MQ=255)
                                                     aCGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGtt                                                                                               <  3:293216/148‑1 (MQ=255)
                                                         tGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCa                                                                                          >  2:853541/1‑149 (MQ=255)
                                                          gAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCa                                                                                          <  4:513509/148‑1 (MQ=255)
                                                             gCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGc                                                                                      <  2:548545/149‑1 (MQ=255)
                                                              cATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGc                                                                                      >  4:669674/1‑148 (MQ=255)
                                                               aTTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCa                                                                                    >  2:400199/1‑149 (MQ=255)
                                                               aTTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCa                                                                                    <  2:676101/149‑1 (MQ=255)
                                                                 tCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAg                                                                                   <  4:109956/148‑1 (MQ=255)
                                                                  cATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGa                                                                                  <  3:348878/148‑1 (MQ=255)
                                                                   atatGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGAc                                                                                 <  1:588314/148‑1 (MQ=255)
                                                                   atatGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGAc                                                                                 >  2:588314/1‑148 (MQ=255)
                                                                     atGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCAt                                                                              >  1:764248/1‑149 (MQ=255)
                                                                     atGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCAt                                                                              <  1:934524/149‑1 (MQ=255)
                                                                     atGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATAAATAATCTCGc                                                                                                                                          <  1:52825/89‑1 (MQ=255)
                                                                      tGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCAGTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATg                                                                             >  2:885069/1‑149 (MQ=255)
                                                                       gagaATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGa                                                                            >  2:575604/1‑149 (MQ=255)
                                                                        agaATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGAc                                                                           >  2:466483/1‑149 (MQ=255)
                                                                          aaTTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCa                                                                         >  2:88564/1‑149 (MQ=255)
                                                                           aTTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCa                                                                         <  4:655492/148‑1 (MQ=255)
                                                                            ttCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGc                                                                       >  2:239793/1‑149 (MQ=255)
                                                                            ttCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGc                                                                       >  2:638044/1‑149 (MQ=255)
                                                                                ctcacCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTg                                                                    <  3:788416/148‑1 (MQ=255)
                                                                                 tcacCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGt                                                                   <  4:402249/148‑1 (MQ=255)
                                                                                 tcacCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTc                                                                  <  1:61795/149‑1 (MQ=255)
                                                                                 tcacCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTc                                                                  >  2:980458/1‑149 (MQ=255)
                                                                                   acCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTc                                                                                                    >  4:412137/1‑113 (MQ=255)
                                                                                   acCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTcac                                                                >  1:455395/1‑149 (MQ=255)
                                                                                    ccGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTcac                                                                >  3:937719/1‑148 (MQ=255)
                                                                                     cGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTcaca                                                               <  4:925676/148‑1 (MQ=255)
                                                                                      gTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTACTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCACTGCCAGACCATGACCAGCTTGTCACAg                                                              <  4:878840/148‑1 (MQ=255)
                                                                                          cccGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGa                                                         >  1:737292/1‑149 (MQ=255)
                                                                                           ccGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGa                                                         >  4:1007160/1‑148 (MQ=255)
                                                                                           ccGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAt                                                        >  1:235599/1‑149 (MQ=255)
                                                                                            cGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAt                                                        >  4:454679/1‑148 (MQ=255)
                                                                                              ccAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAtttt                                                     <  2:655678/149‑1 (MQ=255)
                                                                                              ccAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAtttt                                                     >  1:909861/1‑149 (MQ=255)
                                                                                                 gTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGaa                                                  <  1:646683/149‑1 (MQ=255)
                                                                                                 gTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGaa                                                  <  2:461458/149‑1 (MQ=255)
                                                                                                 gTTGGGCGCTGGCACAACAAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGAGATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGa                                                   <  3:50989/148‑1 (MQ=255)
                                                                                                        gcTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCgg                                            <  3:454679/148‑1 (MQ=255)
                                                                                                            gCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGa                                       >  2:752932/1‑149 (MQ=255)
                                                                                                                 aCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTg                                                                    <  3:507458/115‑1 (MQ=255)
                                                                                                                 aCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTg                                                                    >  4:507458/1‑115 (MQ=255)
                                                                                                                 aCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCgg                                  >  1:338537/1‑149 (MQ=255)
                                                                                                                 aCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCgg                                  >  2:117518/1‑149 (MQ=255)
                                                                                                                 aCGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCg                                   >  3:113083/1‑148 (MQ=255)
                                                                                                                   gAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTa                                <  1:205428/149‑1 (MQ=255)
                                                                                                                   gAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTa                                >  2:593143/1‑149 (MQ=255)
                                                                                                                     aTTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTAtt                              <  1:593143/149‑1 (MQ=255)
                                                                                                                      ttGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTc                             >  2:688285/1‑149 (MQ=255)
                                                                                                                       tGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTcc                            <  2:1020638/149‑1 (MQ=255)
                                                                                                                       tGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTcc                            <  2:1073991/149‑1 (MQ=255)
                                                                                                                       tGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTcc                            <  1:172293/149‑1 (MQ=255)
                                                                                                                          ttATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATc                         <  1:117518/149‑1 (MQ=255)
                                                                                                                          ttATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATc                         <  1:159097/149‑1 (MQ=255)
                                                                                                                          ttATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATc                         >  2:715954/1‑149 (MQ=255)
                                                                                                                           tATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATc                         <  4:443185/148‑1 (MQ=255)
                                                                                                                           tATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCa                        <  1:980458/149‑1 (MQ=255)
                                                                                                                              ctctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGa                      <  4:605243/148‑1 (MQ=255)
                                                                                                                              ctctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGGTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGaa                     <  2:807712/149‑1 (MQ=255)
                                                                                                                                 tATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATAc                  >  1:131688/1‑149 (MQ=255)
                                                                                                                                  aTTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGaa                                      <  3:754429/128‑1 (MQ=255)
                                                                                                                                     cAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCgg               <  4:496481/148‑1 (MQ=255)
                                                                                                                                     cAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGT‑‑CAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGa            <  1:38626/149‑1 (MQ=255)
                                                                                                                                       gagGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGt                                                                                                <  3:549417/65‑1 (MQ=255)
                                                                                                                                        agGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAg                                                 <  3:871222/111‑1 (MQ=255)
                                                                                                                                        agGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAg                                                 >  4:871222/1‑111 (MQ=255)
                                                                                                                                              agaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCa      >  4:329401/1‑148 (MQ=255)
                                                                                                                                                 ttatttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGaa                                      >  4:468789/1‑113 (MQ=255)
                                                                                                                                                 ttatttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGaa                                      <  3:468789/113‑1 (MQ=255)
                                                                                                                                                 ttatttatCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACa   >  4:180031/1‑148 (MQ=255)
                                                                                                                                                  tatttatATCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAg  <  4:700493/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAG  >  NC_000913/3458270‑3458563

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