Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,463,457 C→T 5.0% A178A (GCC→GCT gspK → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,463,4570CT 5.0% 318.9 / 2.5 98A178A (GCC→GCTgspKgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (55/38);  new base T (2/3);  total (57/41)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

GCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGC  >  NC_000913/3463310‑3463605
                                                                                                                                                   |                                                                                                                                                    
gctgctGATTAAAGCCGGTATCGATCGCGGAAACACCGATTAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCaa                                                                                                                                                    <  2:1065572/150‑1 (MQ=255)
 ctgctgATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAAt                                                                                                                                                   >  2:488679/1‑150 (MQ=255)
   gctgATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCa                                                                                                                                                 <  2:1672646/150‑1 (MQ=255)
   gctgATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCa                                                                                                                                                 >  2:490012/1‑150 (MQ=255)
   gctgATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCa                                                                                                                                                 <  2:514098/150‑1 (MQ=255)
   gctgATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCa                                                                                                                                                 >  2:79139/1‑150 (MQ=255)
     tgATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGa                                                                                                                                               >  1:1676314/1‑150 (MQ=255)
      gATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGAt                                                                                                                                              <  1:1037223/150‑1 (MQ=255)
      gATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGAt                                                                                                                                              >  2:119213/1‑150 (MQ=255)
           aCGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCttt                                                                                                                                         >  1:1625082/1‑150 (MQ=255)
              ccGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTc                                                                                                                                      <  2:878868/150‑1 (MQ=255)
                   atcgatcgCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACt                                                                                                                                 <  2:68753/150‑1 (MQ=255)
                      gatcgCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGc                                                                                                                              >  1:37770/1‑150 (MQ=255)
                       atcgCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCg                                                                                                                             >  2:168212/1‑150 (MQ=255)
                       atcgCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCg                                                                                                                             >  2:1559417/1‑150 (MQ=255)
                             gAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTAc                                                                                                                       >  1:989761/1‑150 (MQ=255)
                                acacCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTc                                                                                                                    >  2:1632002/1‑150 (MQ=255)
                                acacCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTc                                                                                                                    <  1:418060/150‑1 (MQ=255)
                                  cccGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGAGGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAg                                                                                                                  <  1:657844/149‑1 (MQ=255)
                                   ccGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGa                                                                                                                 >  1:603761/1‑150 (MQ=255)
                                     gATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATa                                                                                                               >  1:397727/1‑150 (MQ=255)
                                       tGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATaaa                                                                                                             >  1:590705/1‑150 (MQ=255)
                                         agaTCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCCCGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAg                                                                                                           >  2:738695/1‑150 (MQ=255)
                                         agaTCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAg                                                                                                           >  2:1541033/1‑150 (MQ=255)
                                         agaTCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAg                                                                                                           >  1:626185/1‑150 (MQ=255)
                                           ttCGTTCAATCAATGGCTGAGAATATTGATGTCGATTACTCACCACGCTTTCACGTCGTAGAAGATAGCTTTTACCAAAGGGAGACGCCGCCCCGGCATAGTGCCAATCAGATTCTTTTTCTGAGTGGCGAATTACGTCTGATAAAAGGc                                                                                                         <  2:1508776/149‑1 (MQ=255)
                                           aTCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGc                                                                                                         <  2:626185/150‑1 (MQ=255)
                                           aTCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGc                                                                                                         >  2:851903/1‑150 (MQ=255)
                                              gTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATa                                                                                                      >  2:361716/1‑150 (MQ=255)
                                                 caatcaatCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTa                                                                                                                        <  2:1198520/129‑1 (MQ=255)
                                                 caatcaatCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTa                                                                                                                        >  1:1198520/1‑129 (MQ=255)
                                                   atcaatcGCTGACTATATTGATGTCGATGACTCCCCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGa                                                                                                 >  1:143191/1‑150 (MQ=255)
                                                   atcaatcGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTAACAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCCAAATACGTCAGATAAAAGGCATAACCGa                                                                                                 >  1:289494/1‑150 (MQ=255)
                                                      aatcGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGaaaa                                                                                              <  1:605523/150‑1 (MQ=255)
                                                       atcGCTGAGTATATTGATGTCGATGACTCACGTCGCTTTGACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATTGTGCCATTCAGATGCTTTTTCTGACTGGCGAATTACGTGAGATAATAGGCATAACCGaaatc                                                                                             <  2:1127640/150‑3 (MQ=255)
                                                       atcGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAAc                                                                                             >  2:444401/1‑150 (MQ=255)
                                                        tcGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACa                                                                                            <  1:1632002/150‑1 (MQ=255)
                                                          gCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATc                                                                                          <  1:79139/150‑1 (MQ=255)
                                                            tGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTa                                                                                        <  1:168212/150‑1 (MQ=255)
                                                             gACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATACCCGAAAACATCTAt                                                                                       >  1:719459/1‑150 (MQ=255)
                                                               cTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCa                                                                                     >  1:197104/1‑150 (MQ=255)
                                                                  tatTGATGTCGATGACTCACCACGCTTTCACGGCACAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCg                                                                                  <  2:938868/150‑1 (MQ=255)
                                                                     tGATGTCGATGACTCACCACGCTTTCATGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCt                                                                               <  2:1483451/150‑1 (MQ=255)
                                                                     tGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCt                                                                               >  1:1035184/1‑150 (MQ=255)
                                                                     tGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCt                                                                               <  2:523015/150‑1 (MQ=255)
                                                                     tGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCt                                                                               <  2:1097130/150‑1 (MQ=255)
                                                                       aTGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTa                                                                             <  1:1214181/150‑1 (MQ=255)
                                                                         gTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATc                                                                           >  1:511560/1‑150 (MQ=255)
                                                                          tCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATcc                                                                          <  1:361716/150‑1 (MQ=255)
                                                                             ttGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGt                                                                       <  2:357411/149‑1 (MQ=255)
                                                                             aTGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTTCCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCGGACTGGCGAATTACGGCAGATAAAAGGCATAACCGAACACATCTATCAGCGGCTTATCACGt                                                                       >  2:637855/1‑150 (MQ=255)
                                                                                 cTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGt                                                                   >  2:718004/1‑150 (MQ=255)
                                                                                     ccacGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCTAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGc                                                               >  1:878878/1‑150 (MQ=255)
                                                                                     ccacGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCTAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGc                                                               >  1:872884/1‑150 (MQ=255)
                                                                                     ccacGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTg                                                                <  1:234788/149‑1 (MQ=255)
                                                                                     ccacGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTg                                                                >  2:234788/1‑149 (MQ=255)
                                                                                      cacGCTTTCACGGCGCAGTAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCg                                                              <  2:731441/150‑1 (MQ=255)
                                                                                      cacGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCg                                                              <  2:37770/150‑1 (MQ=255)
                                                                                      cacGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCg                                                              >  1:1717336/1‑150 (MQ=255)
                                                                                       acGCTTTCACGGCGCAGAAGTTAGCTTTTACCAAAGCCAGACGCCGCCCCGGCTTAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGGATGTCTGCGt                                                             <  2:103842/150‑1 (MQ=255)
                                                                                         gCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTtc                                                           <  1:488679/150‑1 (MQ=255)
                                                                                          cTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTtct                                                          >  2:797417/1‑150 (MQ=255)
                                                                                            ttCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTcc                                                        >  1:764084/1‑150 (MQ=255)
                                                                                            ttCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTcc                                                        <  2:764084/150‑1 (MQ=255)
                                                                                             tCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTccc                                                       >  2:1158769/1‑150 (MQ=255)
                                                                                               aCGTCGCAGAAGATAGATTTTACCAAAGTCAGACGCCGCCCCGGATTAGTGCCAATCAGATGCTTTATCTGACTTGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGa                                                     <  2:884340/150‑1 (MQ=255)
                                                                                                cGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGAc                                                    >  1:881609/1‑150 (MQ=255)
                                                                                                cGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGAc                                                    >  2:1615835/1‑150 (MQ=255)
                                                                                                   cgcAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCAc                                                 >  1:859452/1‑150 (MQ=255)
                                                                                                       gaagaTAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGaa                                             >  1:21827/1‑150 (MQ=255)
                                                                                                        aagaTAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAAc                                            <  1:444401/150‑1 (MQ=255)
                                                                                                          gaTAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATATATCAGCGGCTTATCCCGTATGTCTGCGTTATCCCGACCACGGAACtt                                          >  2:1251875/1‑150 (MQ=255)
                                                                                                             aGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGc                                       >  2:147116/1‑150 (MQ=255)
                                                                                                               cTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGCCCACGGAACTTAGCAt                                     >  2:881704/1‑150 (MQ=255)
                                                                                                                    aCCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATc                                <  2:981472/150‑1 (MQ=255)
                                                                                                                    aCCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATc                                >  2:1500942/1‑150 (MQ=255)
                                                                                                                     tcAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCtt                              >  2:944486/2‑150 (MQ=255)
                                                                                                                     ccAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCt                               >  1:1490971/1‑150 (MQ=255)
                                                                                                                     ccAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCt                               >  1:246304/1‑150 (MQ=255)
                                                                                                                      cAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCtt                              >  2:305521/1‑150 (MQ=255)
                                                                                                                         aGCCAGACGCCGCCCCGGCATAGTGACAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAat                           <  2:509893/150‑1 (MQ=255)
                                                                                                                          gCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAata                          >  1:1505289/1‑150 (MQ=255)
                                                                                                                           ccAGACGCCGTCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAatat                         >  2:103235/1‑150 (MQ=255)
                                                                                                                             aGACGCCGCCCCGGCATAGTGCTAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGc                       <  2:872884/150‑1 (MQ=255)
                                                                                                                             aGACGCCGCCCCGGCATAGTGCTAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGc                       <  2:878878/150‑1 (MQ=255)
                                                                                                                             aGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCCCGGAACTTAGCATCAATCTTAATATGc                       >  1:136190/1‑150 (MQ=255)
                                                                                                                             aGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGc                       <  1:797417/150‑1 (MQ=255)
                                                                                                                             aGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTAAGATAAAATGCATAACCGAAAACATCTATCAGCGGGTTATCCCGTATGAAAGCTAACTCCCGCCCACCGAACTTAGCATCAATCTTAATATga                       >  2:245072/1‑149 (MQ=255)
                                                                                                                                cgccgcCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTaa                    >  1:1223416/1‑150 (MQ=255)
                                                                                                                                  ccgccCCTGCATAGTGCTAATGAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTGTGCGTTCTCGCGACCACGGAACTTAGCATCAATCTTAATATGCTTACg                  <  2:876549/150‑1 (MQ=255)
                                                                                                                                   cgccCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACgg                 <  1:1559417/150‑1 (MQ=255)
                                                                                                                                        cGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATTTCTGCGTTCTCCCGAACACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAAc            >  1:207579/1‑150 (MQ=255)
                                                                                                                                        cGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAAc            <  2:989761/150‑1 (MQ=255)
                                                                                                                                              aGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATAtt      >  1:684389/1‑150 (MQ=255)
                                                                                                                                              aGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATAtt      >  1:739735/1‑150 (MQ=255)
                                                                                                                                               gTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTc     <  2:476933/150‑1 (MQ=255)
                                                                                                                                               gTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTc     <  2:1467836/150‑1 (MQ=255)
                                                                                                                                                 gCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCg   >  1:266485/1‑150 (MQ=255)
                                                                                                                                                  ccAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGc  <  2:881609/150‑1 (MQ=255)
                                                                                                                                                  ccAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGc  <  2:1356682/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCTGCTGATTAACGCCGGTATCGATCGCGGAAACACCGATGAGATCGTGCAATCAATCGCTGACTATATTGATGTCGATGACTCACCACGCTTTCACGGCGCAGAAGATAGCTTTTACCAAAGCCAGACGCCGCCCCGGCATAGTGCCAATCAGATGCTTTTTCTGACTGGCGAATTACGTCAGATAAAAGGCATAACCGAAAACATCTATCAGCGGCTTATCCCGTATGTCTGCGTTCTCCCGACCACGGAACTTAGCATCAATCTTAATATGCTAACGGAGAACGATATTCCGC  >  NC_000913/3463310‑3463605

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