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

AGTTCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑TTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAA  >  NC_000913/1966934‑1967223
                                                                                                                                                      |                                                                                                                                             
aGTTCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGTGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGCTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGt‑ctt                                                                                                                                             <  2:22410/151‑1 (MQ=255)
     ccGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTtcctc                                                                                                                                        <  2:32859/151‑1 (MQ=255)
         aaaCTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGc                                                                                                                                    <  1:1102994/151‑1 (MQ=255)
         aaaCTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGc                                                                                                                                    >  1:1023155/1‑151 (MQ=255)
          aaCTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGTTGGTTGCTCCATAGTCGCTCCCTCCCATCCCCAG‑TTCTCAAAAGTTTTTGTT‑AGTCTTTTCCTCCa                                                                                                                                      >  1:595035/1‑148 (MQ=255)
           aCTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGt                                                                                                                                  >  1:1169969/1‑151 (MQ=255)
            cTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTc                                                                                                                                 >  2:924766/1‑151 (MQ=255)
             tGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTcg                                                                                                                                >  1:763969/1‑151 (MQ=255)
              gTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCACTGATATCGCCAGCTGAACGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTc                                                                                                                                 <  1:814015/149‑1 (MQ=255)
              gTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTcgc                                                                                                                               <  2:1025846/151‑1 (MQ=255)
                  cgcAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCgcggcg                                                                                                                           <  1:1042242/151‑1 (MQ=255)
                     agcaTACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTa                                                                                                                        <  2:40879/151‑1 (MQ=255)
                      gcaTACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTa                                                                                                                        >  2:394542/1‑150 (MQ=255)
                             cgTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATTGCtt                                                                                                                >  1:620996/1‑151 (MQ=255)
                             cgTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCtt                                                                                                                <  1:53902/151‑1 (MQ=255)
                                cAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTcac                                                                                                             >  1:563392/1‑151 (MQ=255)
                                cAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTca                                                                                                              <  1:1195592/150‑1 (MQ=255)
                                     tGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACcaca                                                                                                         >  1:228041/1‑150 (MQ=255)
                                                  aaTGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGcccc                                                                                           >  1:489572/1‑151 (MQ=255)
                                                   aTGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGccccc                                                                                          <  1:458146/151‑1 (MQ=255)
                                                      atatCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCt                                                                                       <  1:421352/151‑1 (MQ=255)
                                                       tatCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCtt                                                                                      >  2:995496/1‑151 (MQ=255)
                                                       tatCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCtt                                                                                      >  2:525177/1‑151 (MQ=255)
                                                                  gAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCgcagca                                                                           <  2:618428/151‑1 (MQ=255)
                                                                  gAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCgcagc                                                                            <  1:248432/150‑1 (MQ=255)
                                                                   aaTGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCgcagca                                                                           >  2:1043907/1‑150 (MQ=255)
                                                                   aaTGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCaa                                                                          <  2:1016223/151‑1 (MQ=255)
                                                                    aTGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAat                                                                         >  2:527081/1‑151 (MQ=255)
                                                                     tGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAatg                                                                        >  2:154285/1‑151 (MQ=255)
                                                                      gCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAatg                                                                        <  2:60415/150‑1 (MQ=255)
                                                                                  gTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTtctc                                                                <  1:485631/146‑1 (MQ=255)
                                                                                  gTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTtctc                                                                >  2:485749/1‑146 (MQ=255)
                                                                                  gTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCtctttct                                                           >  2:198438/1‑151 (MQ=255)
                                                                                     tgattgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTttcttc                                                         <  1:243188/150‑1 (MQ=255)
                                                                                     tgattgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTttcttc                                                         <  1:533286/150‑1 (MQ=255)
                                                                                       attgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCt                                                      <  2:606603/151‑1 (MQ=255)
                                                                                         tgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTc                                                    >  2:537592/1‑151 (MQ=255)
                                                                                           atGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTg                                                  <  1:935037/151‑1 (MQ=255)
                                                                                              gTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCCCCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAg                                               <  1:313744/151‑1 (MQ=255)
                                                                                              gTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCTCATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAg                                               <  1:948916/151‑1 (MQ=255)
                                                                                                 gCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCTCCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCa                                            <  1:376953/151‑1 (MQ=255)
                                                                                                 gCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTc                                             >  2:164129/1‑150 (MQ=255)
                                                                                                 gCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTc                                             <  2:7839/150‑1 (MQ=255)
                                                                                                     catAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCAt                                        >  2:564777/1‑151 (MQ=255)
                                                                                                          tCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAac                                    >  2:930113/1‑150 (MQ=255)
                                                                                                            gCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACt                                 >  1:497720/1‑151 (MQ=255)
                                                                                                            gCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACt                                 <  1:1089909/151‑1 (MQ=255)
                                                                                                                  tCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAAt                           >  1:514146/1‑151 (MQ=255)
                                                                                                                      aTGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGcc                        >  1:1039533/1‑150 (MQ=255)
                                                                                                                      aTGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCg                       >  2:523327/1‑151 (MQ=255)
                                                                                                                       tGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCg                       >  1:512153/1‑150 (MQ=255)
                                                                                                                        gCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGAc                     >  2:697596/1‑151 (MQ=255)
                                                                                                                          ccAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACAt                   <  2:44286/151‑1 (MQ=255)
                                                                                                                          ccAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACAt                   >  1:141078/1‑151 (MQ=255)
                                                                                                                           cAGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATc                  <  1:5310/151‑1 (MQ=255)
                                                                                                                            aGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATcg                 >  2:1030648/1‑151 (MQ=255)
                                                                                                                            aGTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATc                  >  1:457936/1‑150 (MQ=255)
                                                                                                                             gTTTCTC‑AAGGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATcgc                >  1:871930/1‑151 (MQ=255)
                                                                                                                             gTTTCTC‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTTGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGATATcgc                >  1:80769/1‑151 (MQ=255)
                                                                                                                                 ctc‑AAAGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGcc             >  2:853109/1‑150 (MQ=255)
                                                                                                                                     aaaGATTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATcc         >  2:450551/1‑151 (MQ=255)
                                                                                                                                         aTTTTGTTGAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCAc     <  2:578175/151‑1 (MQ=255)
                                                                                                                                            ttgttgAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTTTTCGCTTTTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGa   >  1:980024/1‑150 (MQ=255)
                                                                                                                                            ttgttgAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGa   >  1:484441/1‑150 (MQ=255)
                                                                                                                                            ttgttgAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGa   >  2:504522/1‑150 (MQ=255)
                                                                                                                                             tgttgAGT‑CTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGaa  <  2:1189371/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
AGTTCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTC‑AAAGATTTTGTTGAGT‑TTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAA  >  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: ‑