Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,319,371 G→A 6.2% T169I (ACC→ATC)  rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,319,3710GA 6.2% 261.0 / 5.6 81T169I (ACC→ATC) rcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base G (36/40);  new base A (2/3);  total (38/43)
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.33e-01

CAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAG  >  NC_000913/2319232‑2319513
                                                                                                                                           |                                                                                                                                              
cAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGt                                                                                                                                        <  1:700629/148‑1 (MQ=255)
 aGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTc                                                                                                                                      >  3:584741/1‑149 (MQ=255)
   aTAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGATACCGCGttgaa                                                                                                                                    <  3:919666/149‑4 (MQ=255)
   aTAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGATACCGCGttga                                                                                                                                     <  2:520925/148‑3 (MQ=255)
       aCTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTgg                                                                                                                                <  4:874227/149‑1 (MQ=255)
        cTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGc                                                                                                                               >  4:786856/1‑149 (MQ=255)
        cTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGc                                                                                                                               <  3:871584/149‑1 (MQ=255)
        cTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGc                                                                                                                               <  3:803816/149‑1 (MQ=255)
           gCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCata                                                                                                                            <  4:43360/149‑1 (MQ=255)
              tCAACTCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATc                                                                                                                          >  2:684191/1‑148 (MQ=255)
               cAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCAc                                                                                                                        <  3:786856/149‑1 (MQ=255)
                 aCGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACg                                                                                                                       >  1:612281/1‑148 (MQ=255)
                  cgcgTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCa                                                                                                                     <  3:63151/149‑1 (MQ=255)
                    cgTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAg                                                                                                                    <  1:181724/148‑1 (MQ=255)
                     gTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGa                                                                                                                   <  1:285385/148‑1 (MQ=255)
                      tAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGAcc                                                                                                                 >  4:754571/1‑149 (MQ=255)
                      tAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGAc                                                                                                                  >  1:663901/1‑148 (MQ=255)
                       aaaaaTACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGTGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCg                                                                                                                >  4:745609/1‑149 (MQ=255)
                       aaaaaTACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCata                                                                                                                            >  1:775447/1‑137 (MQ=255)
                       aaaaaTACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCata                                                                                                                            <  2:775447/137‑1 (MQ=255)
                         aaaTACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGaa                                                                                                              >  4:734395/1‑149 (MQ=255)
                            tACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAAt                                                                                                            <  1:954892/148‑1 (MQ=255)
                             aCCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAAtt                                                                                                           <  2:875172/148‑1 (MQ=255)
                               ccGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGc                                                                                                        <  3:457916/149‑1 (MQ=255)
                                cGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGTAGACCGAAATTGGcc                                                                                                       <  4:913130/149‑1 (MQ=255)
                                cGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAAt                                                                                                            >  1:875504/1‑144 (MQ=255)
                                cGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAAt                                                                                                            <  2:875504/144‑1 (MQ=255)
                                cGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGcc                                                                                                       <  3:925276/149‑1 (MQ=255)
                                cGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGc                                                                                                        >  1:558840/1‑148 (MQ=255)
                                     ccAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGc                                                                                                   <  1:8956/148‑1 (MQ=255)
                                       aGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCag                                                                                                 <  1:964152/148‑1 (MQ=255)
                                       aGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCag                                                                                                 >  1:568938/1‑148 (MQ=255)
                                        gggCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCaga                                                                                                <  2:458577/148‑1 (MQ=255)
                                                      gCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGa                                                                                 >  3:836217/1‑149 (MQ=255)
                                                         gATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATACGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGAtt                                                                               >  2:305968/1‑148 (MQ=255)
                                                          aTCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTaa                                                                             >  3:245555/1‑149 (MQ=255)
                                                           tCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTaa                                                                             <  2:683887/148‑1 (MQ=255)
                                                           tCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGATACCGCGttgaa                                                                                                                                    <  2:288458/93‑4 (MQ=255)
                                                           tCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGATACCGCGttgaa                                                                                                                                    >  1:288458/1‑90 (MQ=255)
                                                                  gagGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATAccc                                                                     >  4:499904/1‑149 (MQ=255)
                                                                           gCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCCGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGAt                                                            <  3:587829/149‑1 (MQ=255)
                                                                                tATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTa                                                        <  2:663901/148‑1 (MQ=255)
                                                                                 atcatcTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGaaact                                                                                                           >  2:490491/1‑94 (MQ=255)
                                                                                 atcatcTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGaaact                                                                                                           <  1:490491/96‑3 (MQ=255)
                                                                                 atcatcTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTaa                                                       >  1:434762/1‑148 (MQ=255)
                                                                                 ataatcTTGTGGTGCATTTCGATATATATCCATGCGTTCAACCAAAGCTTTCAACGCGGTACGGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCAGGCCGAGGTTATCGCTGCAGATTAAAAATACCCGGTTTAGATAGTaa                                                       <  2:608793/145‑1 (MQ=255)
                                                                                   catcTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGc                                                     >  2:226963/1‑148 (MQ=255)
                                                                                    atcTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGGTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGcc                                                    >  2:613576/1‑148 (MQ=255)
                                                                                     tcTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAAtt                                                                                                           >  4:315494/1‑92 (MQ=255)
                                                                                     tcTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAAtt                                                                                                           <  3:315494/92‑1 (MQ=255)
                                                                                           ggtgCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGaaa                                             <  1:930532/148‑1 (MQ=255)
                                                                                            gtgCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGaaaa                                            <  1:226963/148‑1 (MQ=255)
                                                                                             tgCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGaaaaa                                           >  1:753262/1‑148 (MQ=255)
                                                                                             tgCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAAt                                          >  4:793298/1‑149 (MQ=255)
                                                                                              gCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAAt                                          >  1:375243/1‑148 (MQ=255)
                                                                                                 tttCGATAGTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTAt                                       <  2:39968/148‑1 (MQ=255)
                                                                                                   tCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATcgc                                    <  4:895207/149‑1 (MQ=255)
                                                                                                     gATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATcgc                                    >  4:183349/1‑147 (MQ=255)
                                                                                                     gATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATcgc                                    <  3:183349/147‑1 (MQ=255)
                                                                                                     gATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCgcgc                                  <  3:754571/149‑1 (MQ=255)
                                                                                                             ttGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAgcg                          >  4:790655/1‑149 (MQ=255)
                                                                                                                aTGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCAt                       >  3:279085/1‑149 (MQ=255)
                                                                                                                 tGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCaga                                                                                                <  2:946192/75‑1 (MQ=255)
                                                                                                                 tGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCaga                                                                                                >  1:946192/1‑75 (MQ=255)
                                                                                                                 tGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATa                      >  4:275894/1‑149 (MQ=255)
                                                                                                                    gTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAAc                   <  4:836217/149‑1 (MQ=255)
                                                                                                                     ttCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAAc                   <  2:105112/148‑1 (MQ=255)
                                                                                                                       cATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCa                 >  2:378346/1‑148 (MQ=255)
                                                                                                                           cAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGcc             >  1:683501/1‑148 (MQ=255)
                                                                                                                            aaaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCaa           >  4:173562/1‑149 (MQ=255)
                                                                                                                            aaaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCaa           >  3:77501/1‑149 (MQ=255)
                                                                                                                            aaaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCa            >  1:112651/1‑148 (MQ=255)
                                                                                                                             aaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCaa           >  1:553317/1‑148 (MQ=255)
                                                                                                                             aaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCaa           <  2:965436/148‑1 (MQ=255)
                                                                                                                             aaGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCAGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAAt          <  4:275272/149‑1 (MQ=255)
                                                                                                                               gCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATg         <  1:378346/148‑1 (MQ=255)
                                                                                                                                 tttCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGAc       >  1:760776/1‑148 (MQ=255)
                                                                                                                                   tCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTc     <  2:375243/148‑1 (MQ=255)
                                                                                                                                   ggaaCGCGATACCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTcc    >  3:510541/3‑149 (MQ=255)
                                                                                                                                    cAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTcc    <  1:727288/148‑1 (MQ=255)
                                                                                                                                     aaCGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCag  <  3:918037/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
CAGATAAACTGGCGTCAACGCGTAAAAATACCCGACGCCAGGGCGCGGACCTTCGCTGATCCAGTAGAGGTTACTGCCGCTATCATCTTGTGGTGCATTTCGATATTTATTGATGCGTTCATGCAAAGCTTTCAACGCGGTATCGCGTTCCACTGGCATATCACGCAGACCGAAATTGGCCATGCAGAGGTTATCGCTGCCGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAG  >  NC_000913/2319232‑2319513

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