Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,116,854 A→C 5.8% K16N (AAA→AAC)  yceA → UPF0176 protein YceA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,116,8540AC 5.8% 224.1 / 5.1 69K16N (AAA→AAC) yceAUPF0176 protein YceA
Reads supporting (aligned to +/- strand):  ref base A (25/40);  new base C (2/2);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGG  >  NC_000913/1116712‑1116997
                                                                                                                                              |                                                                                                                                               
gcAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTgg                                                                                                                                          >  2:706044/1‑150 (MQ=255)
gcAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTgg                                                                                                                                          <  1:1254607/150‑1 (MQ=255)
gcAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCT‑AAATCCAAAATGTTGGc                                                                                                                                         <  2:21875/150‑1 (MQ=255)
  aaaTGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCt                                                                                                                                        <  1:1262437/150‑1 (MQ=255)
  aaaTGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCt                                                                                                                                        <  1:623650/150‑1 (MQ=255)
   aaTGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTg                                                                                                                                       <  2:907823/150‑1 (MQ=255)
   aaTGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTg                                                                                                                                       <  1:1374277/150‑1 (MQ=255)
     tGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTgag                                                                                                                                     >  1:1371456/1‑150 (MQ=255)
      gTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTgaga                                                                                                                                    <  1:250953/150‑1 (MQ=255)
       tAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAACATTTTGGCttttag                                                                                                                                   <  1:397670/150‑6 (MQ=255)
       tAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGagag                                                                                                                                   <  1:784822/150‑1 (MQ=255)
        aGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgc                                                                                                                                  <  1:500921/150‑1 (MQ=255)
         gCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAACATTTTGGCttttagcg                                                                                                                                 <  1:546555/150‑8 (MQ=255)
           gTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAgcgaa                                                                                                                               <  2:872031/150‑1 (MQ=255)
                     tgGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAacca                                                                                                                     >  2:478977/1‑150 (MQ=255)
                      ggATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAaccac                                                                                                                    >  2:885844/1‑150 (MQ=255)
                      ggATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAaccac                                                                                                                    >  2:404898/1‑150 (MQ=255)
                       gATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAaccacc                                                                                                                   <  1:514282/150‑1 (MQ=255)
                        aggTTAATTATCGAGAATGGCTATATCATGCCGCGGAGTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGTGAGCGAACCGCGAAccacca                                                                                                                  <  2:1311593/148‑1 (MQ=255)
                            tAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAGCCGCGGAGCACCATTTc                                                                                                              >  2:314252/1‑150 (MQ=255)
                             aaTTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCg                                                                                                             >  2:425466/1‑150 (MQ=255)
                              aTTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGt                                                                                                            <  1:1030112/150‑1 (MQ=255)
                               ttATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGtt                                                                                                           <  1:1085241/150‑1 (MQ=255)
                                 aTCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGCCGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGtttt                                                                                                         <  2:248271/150‑1 (MQ=255)
                                 aTCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGtttt                                                                                                         <  1:404898/150‑1 (MQ=255)
                                 aTCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGtttt                                                                                                         <  1:709556/150‑1 (MQ=255)
                                  tCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTa                                                                                                        <  2:1371456/150‑1 (MQ=255)
                                              atatCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCca                                                                                            >  1:1315233/1‑150 (MQ=255)
                                              atatCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCca                                                                                            <  1:433161/150‑1 (MQ=255)
                                               tatCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCcac                                                                                           >  1:1346473/1‑150 (MQ=255)
                                                atCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCcaca                                                                                          >  1:775583/1‑150 (MQ=255)
                                                    tGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGc                                                                                      <  1:706044/150‑1 (MQ=255)
                                                    tGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGc                                                                                      >  2:121559/1‑150 (MQ=255)
                                                     gCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTAAACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGcc                                                                                     >  2:531666/1‑150 (MQ=255)
                                                            aTTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTa                                                                              <  1:1382807/150‑1 (MQ=255)
                                                              ttttACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAg                                                                            <  2:293454/150‑1 (MQ=255)
                                                              ttttACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAg                                                                            <  1:79154/150‑1 (MQ=255)
                                                               tttACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAgg                                                                           <  1:885844/150‑1 (MQ=255)
                                                                      tCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCg                                                                    <  2:695687/150‑1 (MQ=255)
                                                                        ccATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTg                                                                  <  1:902545/150‑1 (MQ=255)
                                                                        ccATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTg                                                                  <  1:792706/150‑1 (MQ=255)
                                                                             tcGCAGGATCCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCt                                                             <  1:951570/150‑1 (MQ=255)
                                                                             tcGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCG‑ACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCtt                                                            >  2:398263/1‑150 (MQ=255)
                                                                              cGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCtt                                                            >  1:374260/1‑150 (MQ=255)
                                                                               gCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCttt                                                           <  2:1315233/150‑1 (MQ=255)
                                                                                cAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTta                                                          >  2:507927/1‑150 (MQ=255)
                                                                                cAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTta                                                          <  2:742057/150‑1 (MQ=255)
                                                                                  ggAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCCAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGCCCCGTGACGCTTtata                                                        >  1:817893/1‑150 (MQ=255)
                                                                                   gAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTatat                                                       <  2:775583/150‑1 (MQ=255)
                                                                                     aCCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCa                                                     >  2:555858/1‑150 (MQ=255)
                                                                                                       gtTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGaa                                   <  2:904173/150‑1 (MQ=255)
                                                                                                          acacAACCGAATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGt                                >  2:843115/1‑150 (MQ=255)
                                                                                                           cacaACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGtt                               >  1:1212881/1‑150 (MQ=255)
                                                                                                               aCCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCCCCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTg                           >  2:46384/1‑150 (MQ=255)
                                                                                                                ccGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTgg                          >  1:911301/1‑150 (MQ=255)
                                                                                                                 cGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTggg                         <  1:843115/150‑1 (MQ=255)
                                                                                                                 cGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTggg                         >  1:609806/1‑150 (MQ=255)
                                                                                                                 cGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTggg                         <  1:507927/150‑1 (MQ=255)
                                                                                                                      ttCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAgt                    <  1:761181/150‑1 (MQ=255)
                                                                                                                         cAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTa                 >  2:836792/1‑150 (MQ=255)
                                                                                                                             gacgCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTg             >  1:189966/1‑150 (MQ=255)
                                                                                                                                gcgcTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGgcg          <  2:374260/150‑1 (MQ=255)
                                                                                                                                  gcTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCa        <  1:836792/150‑1 (MQ=255)
                                                                                                                                   ctctcaGCCAACATTTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCAt       >  1:683711/6‑150 (MQ=255)
                                                                                                                                     ctcaGCCAACATTTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGa     >  2:258164/4‑150 (MQ=255)
                                                                                                                                     aaaaGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGa     <  1:1078357/150‑1 (MQ=255)
                                                                                                                                        aGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAggg  <  2:189966/150‑1 (MQ=255)
                                                                                                                                          ccAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCg                                                                    <  2:1292732/82‑1 (MQ=255)
                                                                                                                                          ccAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCg                                                                    >  1:1292732/1‑82 (MQ=255)
                                                                                                                                              |                                                                                                                                               
GCAAATGTAGCGTAAAATGTGTGGATGTTAATTATCGATAATTGCTATATCATGCCGCGGATTTTTACTTTCCCATCTCGCAGGAACCGTACACCATGCCAGTGTTACACAACCGCATTTCCAACGACGCGCTAAAAGCCAAAATGTTGGCTGAGAGCGAACCGCGAACCACCATTTCGTTTTACAAGTATTTCCACATCGCCGATCCTAAGGCGACCCGTGACGCTTTATATCAGCTGTTTACCGCGCTGAATGTTTTTGGGCGAGTGTATCTGGCGCATGAGGG  >  NC_000913/1116712‑1116997

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