Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,087,637 G→A 5.9% E213K (GAA→AAA)  yiiG → DUF3829 family lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,087,6370GA 5.9% 311.4 / 7.3 101E213K (GAA→AAA) yiiGDUF3829 family lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (39/56);  new base A (3/3);  total (42/59)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

Only 100 of 106 total aligned reads displayed.

CATTAACGAAATGGATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGTTTGATGCAGAAGCGGC  >  NC_000913/4087495‑4087785
                                                                                                                                              |                                                                                                                                                    
cATTAACGAAATGGATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTa                                                                                                                                                >  1:362458/1‑149 (MQ=255)
cATTAACGAAATGGATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTa                                                                                                                                                >  2:490382/1‑149 (MQ=255)
cATTAACGAAATGGATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTa                                                                                                                                                <  2:1137884/149‑1 (MQ=255)
cATTAACGAAATGGATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACCCTGCACCAGCCATTCTTAAAAACGCTGGAAGCCTTTGAACCTGTAACCGAATCTTATCATGCGGCGATTCAGGAAATTa                                                                                                                                                >  1:607708/1‑149 (MQ=255)
 atgaaCGAAATGGATAAATATTACAGCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTaa                                                                                                                                               <  1:527918/146‑1 (MQ=255)
     aCGAAATGGATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGAt                                                                                                                                           >  1:1261712/1‑149 (MQ=255)
        aaaTGGATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAg                                                                                                                                        >  1:733611/1‑149 (MQ=255)
        aaaTGGATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAg                                                                                                                                        >  2:813155/1‑149 (MQ=255)
        aaaTGGATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAg                                                                                                                                        >  2:4428/1‑149 (MQ=255)
        aaaTGGATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAg                                                                                                                                        >  1:1043438/1‑149 (MQ=255)
            ggATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATCAATGATAAGCGTc                                                                                                                                    >  1:323887/1‑149 (MQ=255)
                 aaTATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCtt                                                                                                                               >  2:32770/1‑149 (MQ=255)
                       acacCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGaa                                                                                                                         <  1:664040/149‑1 (MQ=255)
                           ccAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGa                                                                                                                     <  1:921506/149‑1 (MQ=255)
                            cAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTgga                                                                                                                    >  2:50129/1‑147 (MQ=255)
                            cAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGaa                                                                                                                    >  1:1201645/1‑149 (MQ=255)
                            cAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTAATGCCGGGGCGCCAAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGaa                                                                                                                    <  1:364242/149‑1 (MQ=255)
                             aGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAGTCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGaaa                                                                                                                   >  2:1159960/1‑149 (MQ=255)
                             aGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGaaa                                                                                                                   >  1:989181/1‑149 (MQ=255)
                               gAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGATGCACCAGAAATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGaaaac                                                                                                                 >  2:257658/1‑148 (MQ=255)
                                 aaaTTATAAAGACGATGCGTGTGCGAAAGGGAAAACGCTGCACCAGACAGTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAata                                                                                                               <  2:164377/149‑1 (MQ=255)
                                  aaTTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAatat                                                                                                              <  1:32770/149‑1 (MQ=255)
                                  aaTTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAatat                                                                                                              >  1:709556/1‑149 (MQ=255)
                                  aaTTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAatat                                                                                                              >  1:1235176/1‑149 (MQ=255)
                                  aaTTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACAAACTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAatat                                                                                                              <  2:1273275/149‑1 (MQ=255)
                                    gtataAAGACGATGCGGTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTg                                                                                                            <  2:506312/148‑1 (MQ=255)
                                     tataAAGACTGTGCTTTTGCGAAAGGGAAAACGCTTCATTAAATTTTCGTAAAAATTCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCTATTCATGAACTTAATTATAAGCGTCAGCTTGCGGAACTAAAAAATATTga                                                                                                           <  2:912210/149‑1 (MQ=255)
                                       ttaacaCGATGCGTTTGTAGCAGTTAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGAAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTgaag                                                                                                         >  2:1255932/6‑149 (MQ=255)
                                       taAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGAAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTgaag                                                                                                         >  2:573748/1‑149 (MQ=255)
                                       taAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGAAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTgaag                                                                                                         >  2:520766/1‑149 (MQ=255)
                                        acagagGATGCGGTTGCGCAAGGGAAAGCGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTGGAACCTGTAGCCGAATCGTAGCATGCGGGGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTgaaga                                                                                                        <  1:1202105/147‑1 (MQ=255)
                                         aaGACGGTGCGTTTGGGATAGGGAAAACGCTGCAACAGAAAATCTTAAAAAATCTGGAAGCCCTTGAACCTGTAGCCGAATCTTATCATGCGGAGATTGAGGAAAGTAATGATAAGCGTCAGCTAAAAAAAAAAAAAAATATTgaagaa                                                                                                       <  2:335207/149‑1 (MQ=255)
                                         aaGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTgaagaa                                                                                                       <  1:1153655/149‑1 (MQ=255)
                                         aaGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTgaagaa                                                                                                       <  2:1201645/149‑1 (MQ=255)
                                          aGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGaagaaa                                                                                                      >  2:433491/1‑149 (MQ=255)
                                          aGACGAGGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGaagaaa                                                                                                      <  2:426977/149‑1 (MQ=255)
                                               aggtctttGAGGAAGGGCAGTCGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCGGTAGCCGAATCGTATCAGGCGGCGCTTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACGGAAAAATATTGAAGAAAGAGaa                                                                                                 <  2:93804/144‑1 (MQ=255)
                                               aTGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTCGTTATCTTGGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGaa                                                                                                 <  1:1006450/149‑1 (MQ=255)
                                                tGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAg                                                                                                <  1:75956/149‑1 (MQ=255)
                                                tGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAg                                                                                                <  2:194118/149‑1 (MQ=255)
                                                tGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCGGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAg                                                                                                <  1:681008/149‑1 (MQ=255)
                                                ggCGTTTGCGAAAGGGAAAACGCTGCACGAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAg                                                                                                <  1:588650/148‑1 (MQ=255)
                                                      tGCGAAAGGGAAACCGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGaaaa                                                                                          <  2:1183682/149‑1 (MQ=255)
                                                        cGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCAAGCGGCGATTCAGGAAATTAATGATAATCGTCAGCGTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACa                                                                                        >  1:569121/1‑149 (MQ=255)
                                                         gAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGGAAACAt                                                                                       <  2:661442/149‑1 (MQ=255)
                                                         gAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAAGAt                                                                                       <  2:1054186/149‑1 (MQ=255)
                                                         gAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACAt                                                                                       >  2:1240481/1‑149 (MQ=255)
                                                          aaaGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACAtt                                                                                      >  2:962458/1‑149 (MQ=255)
                                                          aaaGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACAtt                                                                                      <  1:980618/149‑1 (MQ=255)
                                                          aaaGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACAtt                                                                                      >  1:767546/1‑149 (MQ=255)
                                                           aaGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGACGAAAGAGAATGGATAACATTc                                                                                     <  2:799727/149‑1 (MQ=255)
                                                           aaGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCCTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTc                                                                                     <  1:373429/149‑1 (MQ=255)
                                                           aaGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCCTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTc                                                                                     <  1:659688/149‑1 (MQ=255)
                                                             gggAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCa                                                                                   <  1:987985/149‑1 (MQ=255)
                                                               gAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCTAAACTGTAAAATATTGAAGAAAGAGAAGGGAAAACATTCCact                                                                                 >  1:564655/1‑149 (MQ=255)
                                                                   gcGCTGCACCAGACATTCTTAAAAAATTTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCactacta                                                                             <  2:79041/148‑1 (MQ=255)
                                                                             aGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACCGAAAAATATTGAAGAACGAGAAGGGAAAACATTCCACTACTACTCTCTGGCa                                                                   <  2:781007/149‑1 (MQ=255)
                                                                             aGACATCCTGAAAAAATCTGGAAGCCGTTGACCCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCa                                                                   <  2:1144546/149‑1 (MQ=255)
                                                                               aCATTCTTAAAAAACCTGAAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATCAATGATAAGCCGCAGCTTGCGGAACTGGAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGt                                                                 <  1:1090841/149‑1 (MQ=255)
                                                                                     ttAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGAt                                                           >  2:953654/1‑149 (MQ=255)
                                                                                     gtAAAAAATCAGGAAGCCTTTGAACCGGTAGCCGAATCTTAGCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAATACATTCCACTACTACTCTCTGGCAGTCATGAt                                                           <  1:444496/148‑1 (MQ=255)
                                                                                         aaaaTCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCa                                                       >  1:670059/1‑149 (MQ=255)
                                                                                         aaaaGCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCa                                                       <  2:503809/149‑1 (MQ=255)
                                                                                          aaaTCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAg                                                      >  2:333009/1‑149 (MQ=255)
                                                                                             tCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGa                                                   <  1:838763/149‑1 (MQ=255)
                                                                                             tCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAAAAAGAGAAGGGAAAACATTCCACTACTACTCTCGGGCAGTCATGATTTCAGCGa                                                   >  2:765012/1‑149 (MQ=255)
                                                                                              cTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGaa                                                  >  1:80995/1‑149 (MQ=255)
                                                                                               tGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGaaa                                                 >  1:778477/1‑149 (MQ=255)
                                                                                               gggAAGCCTTTGAACCTGTAGCCGAAGGTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATGCCACTACTACTCTCTGGCAGTCATGATTTCAGCGaaa                                                 <  2:238234/148‑1 (MQ=255)
                                                                                                ggAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAAc                                                >  2:804138/1‑149 (MQ=255)
                                                                                                      cTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATta                                          <  2:80995/149‑1 (MQ=255)
                                                                                                               tGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGa                                 <  1:1253772/149‑1 (MQ=255)
                                                                                                               tGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGa                                 <  1:741667/149‑1 (MQ=255)
                                                                                                                    ccGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAGAATATTGAAGAAAGAGCAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCCGCTAAACAAATTAATAACCTGATATCg                            >  2:360312/1‑149 (MQ=255)
                                                                                                                    ccGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCg                            >  2:104768/1‑149 (MQ=255)
                                                                                                                     cGAATCTTTTCATGCGGCGATTCAGGAAATTAATGACAAGCCTCAGCTTGCGGAACTGAAAAACATTGAAGAAAGAGCAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGc                           <  1:490382/149‑1 (MQ=255)
                                                                                                                     cGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGc                           <  2:362458/149‑1 (MQ=255)
                                                                                                                      gAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCa                          <  2:767546/149‑1 (MQ=255)
                                                                                                                       aaTCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCaa                         >  1:1247035/1‑149 (MQ=255)
                                                                                                                       aaTCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCaa                         <  2:117575/149‑1 (MQ=255)
                                                                                                                            tgtggtGCGGCGATTCAGGAAAGTAACGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAACGAGAAGGGATAACCTTCCACTACTACTCTGTGGCAGTTGTGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATaa                    <  2:1152211/144‑1 (MQ=255)
                                                                                                                            tATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATaa                    <  2:1179459/149‑1 (MQ=255)
                                                                                                                            tATCATGCGGCGATTCAGAAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATaa                    <  1:1255932/149‑1 (MQ=255)
                                                                                                                            tATCATGCGGCGATTCAGAAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATaa                    <  1:520766/149‑1 (MQ=255)
                                                                                                                            tATCATGCGGCGATTCAGAAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATaa                    <  1:573748/149‑1 (MQ=255)
                                                                                                                               caggcggtgATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATCGAAGATAGAGAAGGGAAAACATTTCACTACTATTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGtt                 <  2:440366/146‑1 (MQ=255)
                                                                                                                               cATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGTAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGtt                 <  1:1173673/149‑1 (MQ=255)
                                                                                                                               cATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAACTTAATAACCTGATACCGCAAGATAAGtt                 >  1:143585/1‑149 (MQ=255)
                                                                                                                                aTGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATAGCAAGAAACGtta                >  2:42376/1‑148 (MQ=255)
                                                                                                                                 tGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATCAGTTTg               >  2:684881/1‑149 (MQ=255)
                                                                                                                                   cggcgATTCAGGAAATGAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGTTTGAt             <  1:1243438/149‑1 (MQ=255)
                                                                                                                                    ggcgATTCAGGAAATTAATGATAAGCGTCAGCGTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGTTTGATg            <  1:454245/149‑1 (MQ=255)
                                                                                                                                          tCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGTTTGATGCAGAAg      <  2:1261712/149‑1 (MQ=255)
                                                                                                                                           cAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAACCATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGTTTGATGCAGAAGc     <  1:1056960/149‑1 (MQ=255)
                                                                                                                                            aGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGTTTGATGCAGAAGCg    >  1:1142652/1‑149 (MQ=255)
                                                                                                                                             ggAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGTTTGATGCAGAAGCgg   >  2:784378/1‑149 (MQ=255)
                                                                                                                                             ggAAATTAATGATAAGCGTAAGTTTGCGGAACTGAAAAATATTGCAGAAACAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGTTTGATGCAGAAGCgg   <  2:714629/149‑1 (MQ=255)
                                                                                                                                              gAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGTTTGATGCAGAAGCGGc  >  1:264994/1‑149 (MQ=255)
                                                                                                                                              gAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGTTTGATGCAGAAGCGGc  >  1:559834/1‑149 (MQ=255)
                                                                                                                                              gAAATAAAGGATAAGCGTCAGCTCGCGGAACTGAAAAATATTGAAGAAAGAGAACGTAAAACATTCCACTACTACTCTCTGGCAGGCATGATTTCAGCGAAACAACTTAATAACCTGATATCCCAAGATAAGTTTGATGCAGAAGCGGc  <  1:50227/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CATTAACGAAATGGATAAATATTACACCCAGGAAAATTATAAAGACGATGCGTTTGCGAAAGGGAAAACGCTGCACCAGACATTCTTAAAAAATCTGGAAGCCTTTGAACCTGTAGCCGAATCTTATCATGCGGCGATTCAGGAAATTAATGATAAGCGTCAGCTTGCGGAACTGAAAAATATTGAAGAAAGAGAAGGGAAAACATTCCACTACTACTCTCTGGCAGTCATGATTTCAGCGAAACAAATTAATAACCTGATATCGCAAGATAAGTTTGATGCAGAAGCGGC  >  NC_000913/4087495‑4087785

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