Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,179,966 T→A 5.5% intergenic (+383/‑30) rplA → / → rplJ 50S ribosomal subunit protein L1/50S ribosomal subunit protein L10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,179,9660TA 5.5% 653.1 / 16.3 197intergenic (+383/‑30)rplA/rplJ50S ribosomal subunit protein L1/50S ribosomal subunit protein L10
Reads supporting (aligned to +/- strand):  ref base T (83/103);  new base A (3/8);  total (86/111)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 199 total aligned reads displayed.

AGACCGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAATGACT  >  NC_000913/4179819‑4180112
                                                                                                                                                   |                                                                                                                                                  
aGACCGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGATGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTtc                                                                                                                                                  <  1:693322/150‑1 (MQ=255)
 gACCGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTtct                                                                                                                                                 <  1:710716/150‑1 (MQ=255)
  aCCGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTtctc                                                                                                                                                <  1:992329/150‑1 (MQ=255)
   ccGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTctct                                                                                                                                               >  1:515665/1‑150 (MQ=255)
    cGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTTTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTg                                                                                                                                              <  2:600272/150‑1 (MQ=255)
    cGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTg                                                                                                                                              >  1:834014/1‑150 (MQ=255)
    cGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTg                                                                                                                                              <  1:170571/150‑1 (MQ=255)
      cAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGc                                                                                                                                            <  1:830161/150‑1 (MQ=255)
      cAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGc                                                                                                                                            >  2:475922/1‑150 (MQ=255)
        ggAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGAAAATCTCTggaac                                                                                                                                          <  1:261102/150‑4 (MQ=255)
        ggAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGAAAATCTCTggaac                                                                                                                                          <  2:1412256/150‑4 (MQ=255)
        ggAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGAAAATCTCTggaac                                                                                                                                          <  1:1063207/150‑4 (MQ=255)
        ggAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGAAAATCTCTggaac                                                                                                                                          <  2:1093224/150‑4 (MQ=255)
        ggAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGAAAATCTCTggaac                                                                                                                                          <  2:1067464/150‑4 (MQ=255)
           gTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACa                                                                                                                                       <  2:635553/150‑1 (MQ=255)
                gCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAg                                                                                                                                  >  2:825338/1‑150 (MQ=255)
                   aGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAg                                                                                                                               >  1:115718/1‑150 (MQ=255)
                   aGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAg                                                                                                                               >  2:24539/1‑150 (MQ=255)
                   aGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAg                                                                                                                               <  2:112018/150‑1 (MQ=255)
                    taaaCTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGc                                                                                                                              <  1:632098/149‑1 (MQ=255)
                      aaCTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCaa                                                                                                                            >  2:994431/1‑150 (MQ=255)
                        cTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAg                                                                                                                          >  1:892713/1‑150 (MQ=255)
                        cTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTCAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAg                                                                                                                          <  2:1165078/150‑1 (MQ=255)
                         ttAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGc                                                                                                                         >  2:489840/1‑150 (MQ=255)
                          tAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCt                                                                                                                        >  2:1023434/1‑150 (MQ=255)
                              ccccTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATg                                                                                                                    <  1:1192630/150‑1 (MQ=255)
                              ccccTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATg                                                                                                                    <  2:892713/150‑1 (MQ=255)
                                  tGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGAttt                                                                                                                                                    <  1:871709/114‑1 (MQ=255)
                                  tGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGAttt                                                                                                                                                    >  2:871709/1‑114 (MQ=255)
                                   gCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCttt                                                                                                               <  2:682968/150‑1 (MQ=255)
                                   gCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAAATAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCttt                                                                                                               >  1:961652/1‑150 (MQ=255)
                                      tAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTaaa                                                                                                            <  1:994431/150‑1 (MQ=255)
                                       aGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAAt                                                                                                           >  2:1202678/1‑150 (MQ=255)
                                        gACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCttt                                                                                                               <  2:30707/145‑1 (MQ=255)
                                        gACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATc                                                                                                          <  1:20296/150‑1 (MQ=255)
                                           ggTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTc                                                                                                       <  2:1291063/150‑1 (MQ=255)
                                            gTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCa                                                                                                      <  1:615813/150‑1 (MQ=255)
                                            gTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCa                                                                                                      <  2:20453/150‑1 (MQ=255)
                                             tGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCaa                                                                                                     >  2:213263/1‑150 (MQ=255)
                                             tGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCaa                                                                                                     <  1:24539/150‑1 (MQ=255)
                                              gACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAg                                                                                                    <  2:134837/150‑1 (MQ=255)
                                                cAGAACGCTTAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGac                                                                                                  <  1:571279/150‑1 (MQ=255)
                                                 aGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGaca                                                                                                 >  1:1309024/1‑150 (MQ=255)
                                                  gAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGacaa                                                                                                >  1:582596/1‑150 (MQ=255)
                                                  gAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGa                                                                                                   <  1:1387912/147‑1 (MQ=255)
                                                  gAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGa                                                                                                   >  2:1387912/1‑147 (MQ=255)
                                                            attattCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGc                                                                                          >  1:163016/1‑146 (MQ=255)
                                                              tattCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGAttgt                                                                                    <  1:138570/150‑1 (MQ=255)
                                                              tattCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGAttgt                                                                                    <  1:316928/150‑1 (MQ=255)
                                                                 tCTTTTATATTCTGGCTTGTTTCTGCTCACCTTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGc                                                                                 <  2:1244724/150‑1 (MQ=255)
                                                                 tCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGc                                                                                 <  1:47804/150‑1 (MQ=255)
                                                                  cTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCt                                                                                >  2:62439/1‑150 (MQ=255)
                                                                   ttttATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTg                                                                               <  2:694930/150‑1 (MQ=255)
                                                                    tttATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCGAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGa                                                                              <  2:1339050/150‑1 (MQ=255)
                                                                      tataTTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAg                                                                            <  2:715556/150‑1 (MQ=255)
                                                                       atatTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGAttgttg                                                                                  >  1:645224/1‑143 (MQ=255)
                                                                         atTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCa                                                                         <  2:58440/150‑1 (MQ=255)
                                                                            cTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAg                                                                                                                          <  1:228173/98‑1 (MQ=255)
                                                                            cTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCg                                                                      >  2:703687/1‑150 (MQ=255)
                                                                                cTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGt                                                                  <  1:213263/150‑1 (MQ=255)
                                                                                   gTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGAt                                                                                       <  1:384152/126‑1 (MQ=255)
                                                                                    tttCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGacaaaca                                                                                             <  2:1359504/119‑1 (MQ=255)
                                                                                    tttCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGcc                                                              <  2:515665/150‑1 (MQ=255)
                                                                                     ttCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCa                                                             <  2:582596/150‑1 (MQ=255)
                                                                                     ttCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCa                                                             <  2:76739/150‑1 (MQ=255)
                                                                                      tCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCt                                                                                <  1:692542/130‑1 (MQ=255)
                                                                                      tCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCaa                                                            >  2:523991/1‑150 (MQ=255)
                                                                                        tGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAg                                                          >  1:1245506/1‑150 (MQ=255)
                                                                                         gCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAgg                                                         >  1:740634/1‑150 (MQ=255)
                                                                                            cACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCACACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgc                                                      >  2:993431/1‑150 (MQ=255)
                                                                                            cACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgc                                                      >  2:808910/1‑150 (MQ=255)
                                                                                             aCCGTAATTAATACTCTCTCTCCGTTTGGAGGAGTTAAGTTAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAATCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgcg                                                     <  1:1399708/150‑1 (MQ=255)
                                                                                             aCCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgcg                                                     <  1:62439/150‑1 (MQ=255)
                                                                                              ccGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgcgc                                                    <  2:834014/150‑1 (MQ=255)
                                                                                              ccGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGgcgcgc                                                    >  1:113877/1‑150 (MQ=255)
                                                                                               cGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCt                                                   <  1:1393598/150‑1 (MQ=255)
                                                                                                 tAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGt                                                 >  1:787379/1‑150 (MQ=255)
                                                                                                     tAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGc                                             <  1:537925/150‑1 (MQ=255)
                                                                                                      aaGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCa                                            >  2:830530/1‑150 (MQ=255)
                                                                                                      aaGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCa                                            >  2:139223/1‑150 (MQ=255)
                                                                                                       aGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCag                                           >  2:385064/1‑150 (MQ=255)
                                                                                                               tctcCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCg                                   <  1:116621/150‑1 (MQ=255)
                                                                                                               tctcCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCg                                   <  2:323818/150‑1 (MQ=255)
                                                                                                                ctcCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCgg                                  <  2:172392/150‑1 (MQ=255)
                                                                                                                 tcCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGa                                 >  2:476477/1‑150 (MQ=255)
                                                                                                                        ggaggaGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGt                          <  2:410494/150‑1 (MQ=255)
                                                                                                                        ggaggaGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGt                          >  1:865341/1‑150 (MQ=255)
                                                                                                                         gaggagTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTg                         <  1:1098546/150‑1 (MQ=255)
                                                                                                                                gaagtgaGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAAc                  >  2:366482/1‑150 (MQ=255)
                                                                                                                                gaagtgaGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAAc                  <  1:1023434/150‑1 (MQ=255)
                                                                                                                                   gtgaGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGt               >  1:411748/1‑150 (MQ=255)
                                                                                                                                       ttgccAGAGAAAATCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGAt           >  1:90349/4‑150 (MQ=255)
                                                                                                                                         gccAGAGAAAATCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATaa         >  1:288862/2‑150 (MQ=255)
                                                                                                                                          ccAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATaaa        <  2:856011/150‑1 (MQ=255)
                                                                                                                                           cAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAg             <  2:157684/144‑1 (MQ=255)
                                                                                                                                           cAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAg             >  1:157684/1‑144 (MQ=255)
                                                                                                                                            agagATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAAt      <  1:101818/150‑1 (MQ=255)
                                                                                                                                            agagATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAAt      >  2:905294/1‑150 (MQ=255)
                                                                                                                                             gagaTTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAATg     <  1:1263760/150‑1 (MQ=255)
                                                                                                                                                aTTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAATGACt  <  1:1428616/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGACCGCAGGAGTTTCGCAAGAAACTTAATCCCCTGCGTAGACGGTGACAGAACGCTAAGATTATTCTTTTATATTCTGGCTTGTTTCTGCTCACCGTAATTAAGACGCTCTCTCCGTTTGGAGGAGTGAAGTGAGTTCCAGAGATTTTCTCTGGCAAACATCCAGGAGCAAAGCTAATGGCTTTAAATCTTCAAGACAAACAAGCGATTGTTGCTGAAGTCAGCGAAGTAGCCAAAGGCGCGCTGTCTGCAGTAGTTGCGGATTCCCGTGGCGTAACTGTAGATAAAATGACT  >  NC_000913/4179819‑4180112

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