Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,257,602 A→T N163I (AAT→ATT)  lexA → transcriptional repressor of SOS regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,257,6020AT97.9% 371.2 / ‑0.8 96N163I (AAT→ATT) lexAtranscriptional repressor of SOS regulon
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (46/48);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCC  >  NC_000913/4257455‑4257749
                                                                                                                                                   |                                                                                                                                                    
gcgTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGCTGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTa                                                                                                                                                    >  3:267295/1‑150 (MQ=255)
gcgTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTaa                                                                                                                                                   <  3:879539/151‑1 (MQ=255)
 cgTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTaaa                                                                                                                                                  <  2:5999/151‑1 (MQ=255)
 cgTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTaa                                                                                                                                                   >  3:257078/1‑150 (MQ=255)
  gTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAg                                                                                                                                                 <  2:553580/151‑1 (MQ=255)
    cAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTc                                                                                                                                               <  2:1463720/151‑1 (MQ=255)
     aGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCg                                                                                                                                              <  3:272868/151‑1 (MQ=255)
      gCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGa                                                                                                                                             <  2:1604201/151‑1 (MQ=255)
       cGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGa                                                                                                                                             >  2:1398678/1‑150 (MQ=255)
        gggATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAAc                                                                                                                                           <  3:79091/151‑1 (MQ=255)
         ggATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACt                                                                                                                                          >  3:1350608/1‑151 (MQ=255)
             gTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTg                                                                                                                                      >  2:765619/1‑151 (MQ=255)
             gTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTg                                                                                                                                      <  2:295365/151‑1 (MQ=255)
               cGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGc                                                                                                                                     <  3:723588/150‑1 (MQ=255)
               cGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGc                                                                                                                                     <  3:1166460/150‑1 (MQ=255)
                gATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCGGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGcc                                                                                                                                    <  3:1163121/150‑1 (MQ=255)
                gATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCa                                                                                                                                   <  3:7350/151‑1 (MQ=255)
                 aTGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGcc                                                                                                                                    >  4:10387/1‑149 (MQ=255)
                   gAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGa                                                                                                                                 <  3:1569247/150‑1 (MQ=255)
                   gAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGa                                                                                                                                 >  3:721978/1‑150 (MQ=255)
                    aaaGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGaa                                                                                                                                >  2:1399348/1‑150 (MQ=255)
                      aGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGaaaa                                                                                                                              >  3:183226/1‑150 (MQ=255)
                       gATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATa                                                                                                                            <  2:600672/151‑1 (MQ=255)
                          atCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCg                                                                                                                         >  2:897931/1‑151 (MQ=255)
                          atCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCg                                                                                                                         >  2:1161899/1‑151 (MQ=255)
                              gCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGtt                                                                                                                     >  2:1185483/1‑151 (MQ=255)
                                   atggatggTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAAc                                                                                                                <  2:318068/151‑1 (MQ=255)
                                       atggTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAAt                                                                                                            <  2:754965/151‑1 (MQ=255)
                                            gACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCg                                                                                                       <  3:321287/151‑1 (MQ=255)
                                             aCTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGt                                                                                                      >  2:1197384/1‑151 (MQ=255)
                                                tGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTg                                                                                                    >  2:301518/1‑150 (MQ=255)
                                                  cTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGAcc                                                                                                 >  3:983393/1‑151 (MQ=255)
                                                  cTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGAcc                                                                                                 >  2:1446975/1‑151 (MQ=255)
                                                   tGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCt                                                                                                >  2:1046097/1‑151 (MQ=255)
                                                   tGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCt                                                                                                >  2:1145235/1‑151 (MQ=255)
                                                      cAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCg                                                                                             <  3:959534/151‑1 (MQ=255)
                                                      cAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCg                                                                                             <  3:820767/151‑1 (MQ=255)
                                                      cAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCg                                                                                             <  3:271784/151‑1 (MQ=255)
                                                      cAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCg                                                                                             >  2:1632495/1‑151 (MQ=255)
                                                       aGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGt                                                                                            <  2:329408/151‑1 (MQ=255)
                                                         tGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTc                                                                                           >  2:1163829/1‑150 (MQ=255)
                                                         tGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTc                                                                                           >  3:1213945/1‑150 (MQ=255)
                                                          gCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTca                                                                                          <  3:1206442/150‑1 (MQ=255)
                                                          gCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTca                                                                                          <  2:569699/150‑1 (MQ=255)
                                                              aaaaCTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTcagc                                                                                        >  3:278222/1‑148 (MQ=255)
                                                              aaaaCTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCaga                                                                                     <  2:381426/151‑1 (MQ=255)
                                                               aaaCTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAACAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTcagcag                                                                                      >  3:3830/1‑149 (MQ=255)
                                                                    cAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCtt                                                                                 <  3:516499/149‑1 (MQ=255)
                                                                    cAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTc                                                                                <  3:257838/150‑1 (MQ=255)
                                                                     aGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCAc                                                                              <  3:220959/151‑1 (MQ=255)
                                                                     aGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCAc                                                                              >  2:871899/1‑151 (MQ=255)
                                                                      ggATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCAcc                                                                             >  2:387727/1‑151 (MQ=255)
                                                                      ggATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCAcc                                                                             <  3:109492/151‑1 (MQ=255)
                                                                       gATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCAcc                                                                             >  2:1606242/1‑150 (MQ=255)
                                                                       gATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCAcc                                                                             >  2:1099684/1‑150 (MQ=255)
                                                                          gTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTg                                                                         <  2:1470661/151‑1 (MQ=255)
                                                                            aCGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGa                                                                        <  2:330503/150‑1 (MQ=255)
                                                                              gTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAgg                                                                     <  3:705022/151‑1 (MQ=255)
                                                                                aaCGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGc                                                                   <  2:217751/151‑1 (MQ=255)
                                                                                 aCGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCt                                                                  >  3:638801/1‑151 (MQ=255)
                                                                                   ggTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTgg                                                                <  2:1230158/151‑1 (MQ=255)
                                                                                      cAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCgg                                                             >  3:659964/1‑151 (MQ=255)
                                                                                      cAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCg                                                              <  3:889690/150‑1 (MQ=255)
                                                                                       aGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGt                                                           >  2:631692/1‑151 (MQ=255)
                                                                                        ggTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGtt                                                          >  2:414200/1‑151 (MQ=255)
                                                                                           cGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTgg                                                        <  2:257095/150‑1 (MQ=255)
                                                                                             ttGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGATTgggg                                                      >  3:23136/1‑151 (MQ=255)
                                                                                             ttGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTgggg                                                      >  3:1193322/1‑150 (MQ=255)
                                                                                                tCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTAt                                                  >  2:188186/1‑151 (MQ=255)
                                                                                                 cGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTAt                                                  <  2:1009038/150‑1 (MQ=255)
                                                                                                  gCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTAt                                                  >  2:1260045/1‑149 (MQ=255)
                                                                                                   cACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTTTTCg                                               >  3:116110/1‑151 (MQ=255)
                                                                                                   cACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCg                                               <  2:210695/151‑1 (MQ=255)
                                                                                                    aCGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCg                                               >  2:448182/1‑150 (MQ=255)
                                                                                                      gTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTAtt                                                 >  3:1339395/1‑146 (MQ=255)
                                                                                                           gatgaCGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCg                                       <  2:1443986/151‑1 (MQ=255)
                                                                                                           gatgaCGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCg                                       <  2:116114/151‑1 (MQ=255)
                                                                                                              gaCGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACt                                    <  2:1085726/151‑1 (MQ=255)
                                                                                                               aCGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTg                                   <  3:1185376/151‑1 (MQ=255)
                                                                                                                   aGTTACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTg                               <  2:148951/151‑1 (MQ=255)
                                                                                                                   aGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTgg                                  >  3:376975/1‑148 (MQ=255)
                                                                                                                      tACCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTa                             >  3:1133448/1‑150 (MQ=255)
                                                                                                                       aCCGTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTAGCCACGGCGACTGGCTGTaa                            >  3:1147811/1‑150 (MQ=255)
                                                                                                                          gTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACata                        <  2:1433463/151‑1 (MQ=255)
                                                                                                                          gTTAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACata                        <  2:1320421/151‑1 (MQ=255)
                                                                                                                            tAAGCGCCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTAACata                        >  2:951879/1‑149 (MQ=255)
                                                                                                                               gcgcCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATctct                   >  2:317670/1‑151 (MQ=255)
                                                                                                                                cgcCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAAAATATctct                   >  2:1615658/1‑150 (MQ=255)
                                                                                                                                 gcCTGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTga                 <  2:183237/151‑1 (MQ=255)
                                                                                                                                    tGAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTTGGGGTATTTGCAACGGCGACCGGCTGTTAAATAtct                     >  3:572849/1‑144 (MQ=255)
                                                                                                                                     gAAAAAACAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGAcc             <  3:875872/151‑1 (MQ=255)
                                                                                                                                         aaaCAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGGTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACcgcg          >  2:458436/1‑150 (MQ=255)
                                                                                                                                          aaCAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGAt        <  3:897874/151‑1 (MQ=255)
                                                                                                                                            cAGGGCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATgc      <  2:1407694/151‑1 (MQ=255)
                                                                                                                                                gCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATgccgcc  <  2:926579/151‑1 (MQ=255)
                                                                                                                                                gCATTAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCCGTAACATATCTCTGAGACCGCGATgccgcc  >  2:1642313/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGTCGTTGACCTTCGTCAGCAGAGCTTCACCATTGAAGGGCTGGCGGTTGGGGTTATTCGCAACGGCGACTGGCTGTAACATATCTCTGAGACCGCGATGCCGCC  >  NC_000913/4257455‑4257749

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