Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,967,082 T→C K119R (AAA→AGA)  cheY ← chemotaxis regulator transmitting signal to flagellar motor component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,967,0820TC98.5% 249.0 / NA 65K119R (AAA→AGA) cheYchemotaxis regulator transmitting signal to flagellar motor component
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (36/28);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAA  >  NC_000913/1966934‑1967223
                                                                                                                                                      |                                                                                                                                             
aGTTCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGTGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGCTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGtctt                                                                                                                                             <  2:22410/151‑1 (MQ=255)
     ccGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTtcctc                                                                                                                                        <  2:32859/151‑1 (MQ=255)
         aaaCTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGc                                                                                                                                    <  1:1102994/151‑1 (MQ=255)
         aaaCTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGc                                                                                                                                    >  1:1023155/1‑151 (MQ=255)
          aaCTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGTTGGTTGCTCCATAGTCGCTCCCTCCCATCCCCAGTTCTCAAAAGTTTTTGTTAGTCTTTTCCTCCa                                                                                                                                      >  1:595035/1‑148 (MQ=255)
           aCTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGt                                                                                                                                  >  1:1169969/1‑151 (MQ=255)
            cTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTc                                                                                                                                 >  2:924766/1‑151 (MQ=255)
             tGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTcg                                                                                                                                >  1:763969/1‑151 (MQ=255)
              gTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCACTGATATCGCCAGCTGAACGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTc                                                                                                                                 <  1:814015/149‑1 (MQ=255)
              gTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTcgc                                                                                                                               <  2:1025846/151‑1 (MQ=255)
                  cgcAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCgcggcg                                                                                                                           <  1:1042242/151‑1 (MQ=255)
                     agcaTACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTa                                                                                                                        <  2:40879/151‑1 (MQ=255)
                      gcaTACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTa                                                                                                                        >  2:394542/1‑150 (MQ=255)
                             cgTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATTGCtt                                                                                                                >  1:620996/1‑151 (MQ=255)
                             cgTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCtt                                                                                                                <  1:53902/151‑1 (MQ=255)
                                cAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTcac                                                                                                             >  1:563392/1‑151 (MQ=255)
                                cAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTca                                                                                                              <  1:1195592/150‑1 (MQ=255)
                                     tGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACcaca                                                                                                         >  1:228041/1‑150 (MQ=255)
                                                  aaTGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGcccc                                                                                           >  1:489572/1‑151 (MQ=255)
                                                   aTGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGccccc                                                                                          <  1:458146/151‑1 (MQ=255)
                                                      atatCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCt                                                                                       <  1:421352/151‑1 (MQ=255)
                                                       tatCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCtt                                                                                      >  2:995496/1‑151 (MQ=255)
                                                       tatCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCtt                                                                                      >  2:525177/1‑151 (MQ=255)
                                                                  gAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCgcagca                                                                           <  2:618428/151‑1 (MQ=255)
                                                                  gAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCgcagc                                                                            <  1:248432/150‑1 (MQ=255)
                                                                   aaTGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCgcagca                                                                           >  2:1043907/1‑150 (MQ=255)
                                                                   aaTGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCaa                                                                          <  2:1016223/151‑1 (MQ=255)
                                                                    aTGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAat                                                                         >  2:527081/1‑151 (MQ=255)
                                                                     tGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAatg                                                                        >  2:154285/1‑151 (MQ=255)
                                                                      gCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAatg                                                                        <  2:60415/150‑1 (MQ=255)
                                                                                  gTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTtctc                                                                <  1:485631/146‑1 (MQ=255)
                                                                                  gTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTtctc                                                                >  2:485749/1‑146 (MQ=255)
                                                                                  gTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCtctttct                                                           >  2:198438/1‑151 (MQ=255)
                                                                                     tgattgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTttcttc                                                         <  1:243188/150‑1 (MQ=255)
                                                                                     tgattgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTttcttc                                                         <  1:533286/150‑1 (MQ=255)
                                                                                       attgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCt                                                      <  2:606603/151‑1 (MQ=255)
                                                                                         tgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTc                                                    >  2:537592/1‑151 (MQ=255)
                                                                                           atGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTg                                                  <  1:935037/151‑1 (MQ=255)
                                                                                              gTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCCCCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAg                                               <  1:313744/151‑1 (MQ=255)
                                                                                              gTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCTCATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAg                                               <  1:948916/151‑1 (MQ=255)
                                                                                                 gCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCTCCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCa                                            <  1:376953/151‑1 (MQ=255)
                                                                                                 gCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTc                                             >  2:164129/1‑150 (MQ=255)
                                                                                                 gCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTc                                             <  2:7839/150‑1 (MQ=255)
                                                                                                     catAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCAt                                        >  2:564777/1‑151 (MQ=255)
                                                                                                          tCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAac                                    >  2:930113/1‑150 (MQ=255)
                                                                                                            gCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACt                                 >  1:497720/1‑151 (MQ=255)
                                                                                                            gCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACt                                 <  1:1089909/151‑1 (MQ=255)
                                                                                                                  tCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAAt                           >  1:514146/1‑151 (MQ=255)
                                                                                                                      aTGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGcc                        >  1:1039533/1‑150 (MQ=255)
                                                                                                                      aTGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCg                       >  2:523327/1‑151 (MQ=255)
                                                                                                                       tGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCg                       >  1:512153/1‑150 (MQ=255)
                                                                                                                        gCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGAc                     >  2:697596/1‑151 (MQ=255)
                                                                                                                          ccAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACAt                   <  2:44286/151‑1 (MQ=255)
                                                                                                                          ccAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACAt                   >  1:141078/1‑151 (MQ=255)
                                                                                                                           cAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATc                  <  1:5310/151‑1 (MQ=255)
                                                                                                                            aGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATcg                 >  2:1030648/1‑151 (MQ=255)
                                                                                                                            aGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATc                  >  1:457936/1‑150 (MQ=255)
                                                                                                                             gTTTCTCAAGGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATcgc                >  1:871930/1‑151 (MQ=255)
                                                                                                                             gTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTTGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGATATcgc                >  1:80769/1‑151 (MQ=255)
                                                                                                                                 ctcAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGcc             >  2:853109/1‑150 (MQ=255)
                                                                                                                                     aaaGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATcc         >  2:450551/1‑151 (MQ=255)
                                                                                                                                         aTTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCAc     <  2:578175/151‑1 (MQ=255)
                                                                                                                                            ttgttgAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTTTTCGCTTTTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGa   >  1:980024/1‑150 (MQ=255)
                                                                                                                                            ttgttgAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGa   >  1:484441/1‑150 (MQ=255)
                                                                                                                                            ttgttgAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGa   >  2:504522/1‑150 (MQ=255)
                                                                                                                                             tgttgAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGaa  <  2:1189371/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
AGTTCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAA  >  NC_000913/1966934‑1967223

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