Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,406,380 C→G 6.9% A88P (GCC→CCC) ‡ lrhA ← transcriptional repressor of flagellar, motility and chemotaxis genes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,406,3800CG 6.9% 297.1 / 10.2 98A88P (GCC→CCC) ‡lrhAtranscriptional repressor of flagellar, motility and chemotaxis genes
Reads supporting (aligned to +/- strand):  ref base C (36/55);  new base G (3/4);  total (39/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TACGTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGA  >  NC_000913/2406233‑2406528
                                                                                                                                                   |                                                                                                                                                    
tACGTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCt                                                                                                                                                    >  2:359716/1‑150 (MQ=255)
   gTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCtcat                                                                                                                                                 >  2:1158412/1‑150 (MQ=255)
    tttAACTCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCAAGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCtcatc                                                                                                                                                <  2:448447/150‑1 (MQ=255)
     ttAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCtcatca                                                                                                                                               <  2:764912/150‑1 (MQ=255)
     ttAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCtcatca                                                                                                                                               >  1:430165/1‑150 (MQ=255)
      tAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGTAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGATCATTAATGATGAGCAGGCCTcatcat                                                                                                                                              <  2:278690/150‑1 (MQ=255)
        aCGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTa                                                                                                                                            <  2:1111565/150‑1 (MQ=255)
        aCGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTa                                                                                                                                            >  1:930207/1‑150 (MQ=255)
        aCGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTCAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAAAGATGAGCAGGCCTCATCATTa                                                                                                                                            >  1:311250/1‑150 (MQ=255)
         cgcgGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTATCAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTaa                                                                                                                                           >  1:340211/1‑150 (MQ=255)
         cgcgGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTATGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTaa                                                                                                                                           >  1:450088/1‑150 (MQ=255)
         cgcgGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGAAAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTaa                                                                                                                                           >  2:1270522/1‑150 (MQ=255)
          gcgTACAACCAGCGATAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAATATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGGCCTGCTCATCATTaat                                                                                                                                          <  2:860623/150‑2 (MQ=255)
          gcgGACATCCAGCGCTAATTTCGTATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGGCCTGCTCATCATTaat                                                                                                                                          <  2:249065/150‑2 (MQ=255)
           cgGACATCCAGCGCTAATTTCGTATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGGCCTGCTCATCATTaatg                                                                                                                                         <  2:554892/150‑3 (MQ=255)
           cgGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTTATGATGAGCAGGCCTCATCATTaaaa                                                                                                                                         >  2:674790/1‑150 (MQ=255)
           cgGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTaaaa                                                                                                                                         >  2:1090233/1‑150 (MQ=255)
            ggACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGGCCTGCTCATCATTaatg                                                                                                                                         <  1:569774/149‑3 (MQ=255)
            ggACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGGCCTGCTCATCATTaatg                                                                                                                                         >  2:569774/1‑147 (MQ=255)
            ggACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAAc                                                                                                                                        <  2:542956/150‑1 (MQ=255)
             gACATCCAGCGCTAATTTCGGATAAACCGAACTAACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACg                                                                                                                                       <  2:633007/150‑1 (MQ=255)
                 tCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAATGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAgg                                                                                                                                   <  1:22502/150‑1 (MQ=255)
                 accACCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAgg                                                                                                                                   <  1:68568/149‑1 (MQ=255)
                  ccAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGa                                                                                                                                  >  1:359645/1‑150 (MQ=255)
                     gcgcTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGAttt                                                                                                                               >  1:1173390/1‑150 (MQ=255)
                        cTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTcc                                                                                                                            <  2:164653/150‑1 (MQ=255)
                        cTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTcc                                                                                                                            <  2:74960/150‑1 (MQ=255)
                        cTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCt                                                                                                                           <  1:799322/150‑1 (MQ=255)
                          aaTTTCGGATAACCCGAAATCACGCTATTTAACAGGAAAGTTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATAAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTg                                                                                                                          <  2:935485/150‑1 (MQ=255)
                           aTTTCGGATAAACCGAACTCACGCTATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATAAGACCTGATGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTgg                                                                                                                         <  2:341608/150‑1 (MQ=255)
                           aTTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTgg                                                                                                                         <  2:307814/150‑1 (MQ=255)
                                ggATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTaa                                                                                                                    >  2:517374/1‑150 (MQ=255)
                                 gATAAACCGAACTCACGCGAGTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCGCTTAAGATTAGTGAACATTAATGATGCGCAGGCCTCATCTTTAAAACGCCGGATTTTCCTGGCGTAAc                                                                                                                   >  2:64769/1‑150 (MQ=255)
                                  aTAAACCGAACTCACGCGATTTATCAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTCGCGTAAcc                                                                                                                  >  1:166270/1‑150 (MQ=255)
                                   aaaaCCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCa                                                                                                                 <  1:926156/149‑1 (MQ=255)
                                       ccaaaCTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCaagaa                                                                                                             <  2:1158850/147‑1 (MQ=255)
                                         gAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGt                                                                                                           <  2:15359/150‑1 (MQ=255)
                                          aaCTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGtt                                                                                                          <  1:594152/150‑1 (MQ=255)
                                           aCTCACGCGATTTAACAGGAAAGCTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCGGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTg                                                                                                         <  2:1002711/150‑1 (MQ=255)
                                            cTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGa                                                                                                        <  1:46613/150‑1 (MQ=255)
                                             tCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCTCCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGaa                                                                                                       <  1:257216/150‑1 (MQ=255)
                                              cACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAAt                                                                                                      >  1:242724/1‑150 (MQ=255)
                                               aCGCTATTTAAAAGGAAAGGTAAGATCGAATCGGCAGATTCATCTGAAGCACCTATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATGTTCCTGGCGTAACCAAGAAGTTGAATg                                                                                                     <  2:1199431/150‑1 (MQ=255)
                                               aCGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATg                                                                                                     <  1:498107/150‑1 (MQ=255)
                                                cgcgATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGc                                                                                                    >  2:389815/1‑150 (MQ=255)
                                                  cgATTTAACAGGAAAGGTAAGATCGTATCGGCATATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCa                                                                                                  <  1:1017497/150‑1 (MQ=255)
                                                       tAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTc                                                                                             >  2:1317696/1‑150 (MQ=255)
                                                          cAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGt                                                                                          <  1:720404/150‑1 (MQ=255)
                                                              aaaGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAAc                                                                                      <  1:50558/150‑1 (MQ=255)
                                                              aaaGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAAc                                                                                      <  1:524954/150‑1 (MQ=255)
                                                                  gTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAAAGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGtt                                                                                  <  2:66119/150‑1 (MQ=255)
                                                                      gATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTttgtt                                                                              >  2:1103920/1‑150 (MQ=255)
                                                                       ttCGTATCGGCAGATTCATCTGAATCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTttgttg                                                                             <  2:595654/149‑1 (MQ=255)
                                                                       aTCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTACTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTttgttg                                                                             >  1:1219145/1‑150 (MQ=255)
                                                                        tCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGc                                                                            <  1:477284/150‑1 (MQ=255)
                                                                        tCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGc                                                                            <  2:430165/150‑1 (MQ=255)
                                                                            aTCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCTTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGccc                                                                        >  2:471948/1‑150 (MQ=255)
                                                                            aTCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGAcc                                                                        <  1:1041980/150‑1 (MQ=255)
                                                                             tCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCg                                                                       <  2:928615/150‑1 (MQ=255)
                                                                              cGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGt                                                                      <  1:359716/150‑1 (MQ=255)
                                                                                   gATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTAGCGACCGTGACGa                                                                 >  2:1200371/1‑150 (MQ=255)
                                                                                        aTCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGaa                                                            >  2:949069/1‑150 (MQ=255)
                                                                                         tCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAAc                                                           <  2:923443/150‑1 (MQ=255)
                                                                                            gAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGt                                                        <  1:183845/150‑1 (MQ=255)
                                                                                            gAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGt                                                        <  2:1234192/150‑1 (MQ=255)
                                                                                             aaGCACCGATAGTTAAAACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGtt                                                       <  1:282124/150‑1 (MQ=255)
                                                                                                  ccGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTc                                                  <  2:750472/150‑1 (MQ=255)
                                                                                                  ccGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTc                                                  >  1:407762/1‑150 (MQ=255)
                                                                                                  ccGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTATc                                                  >  1:401914/1‑150 (MQ=255)
                                                                                                   cGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTcc                                                 <  2:523701/150‑1 (MQ=255)
                                                                                                    gATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTccc                                                <  2:930207/150‑1 (MQ=255)
                                                                                                    gATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTccc                                                <  1:395676/150‑1 (MQ=255)
                                                                                                         ttAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGa                                           >  2:328403/1‑150 (MQ=255)
                                                                                                                 cccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGCACAGTTCTTTCCCAACGAGTTGTTcc                                   >  2:144383/1‑150 (MQ=255)
                                                                                                                 cccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTcc                                   <  1:994369/150‑1 (MQ=255)
                                                                                                                 cccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTcc                                   >  2:1132182/1‑150 (MQ=255)
                                                                                                                  ccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCa                                  >  2:1087398/1‑150 (MQ=255)
                                                                                                                   cTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAg                                 >  2:1058690/1‑150 (MQ=255)
                                                                                                                           tACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCa                         <  1:1090233/150‑1 (MQ=255)
                                                                                                                            aCTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCAt                        <  2:1199689/150‑1 (MQ=255)
                                                                                                                               gcaCATTAATGAAGCGCAGGCCTGAGCATTAAAACTCAGGATTTTCCTGGCGTAACAAATAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTg                     <  1:788029/148‑1 (MQ=255)
                                                                                                                               gAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGATTGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTg                     >  1:1255252/1‑150 (MQ=255)
                                                                                                                               gAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTg                     >  2:904098/1‑150 (MQ=255)
                                                                                                                                aaCATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGc                    <  1:911239/150‑1 (MQ=255)
                                                                                                                                   aTTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGa                 <  1:174540/150‑1 (MQ=255)
                                                                                                                                    ttAATGATGAGGCCTGCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGAc                >  2:210316/1‑150 (MQ=255)
                                                                                                                                    ttAATGATGAGGCCTGCTCATAATTAAAACGCAGGATTTACAGGGCGTAACCAAGATGATGAATGCAATGTTCACTTAACAGTTTGTTGCGACCGTGAAGAGCGAACAGTACTTTCCCAACGAGTTGTTCCAGACGCTGCATTtgatgaa                >  2:210320/1‑149 (MQ=255)
                                                                                                                                     tAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACt               <  1:929279/150‑1 (MQ=255)
                                                                                                                                     tAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACt               <  2:299752/150‑1 (MQ=255)
                                                                                                                                     tAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACt               <  2:242724/150‑1 (MQ=255)
                                                                                                                                      aaTGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGCCGCTGCATTTGCTGACtt              >  1:1129739/1‑150 (MQ=255)
                                                                                                                                      aaTGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACtt              >  2:582438/1‑150 (MQ=255)
                                                                                                                                      aaTGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTATTTCCCAACGAGGTGTTCCAGACGCTGCATTTGCTGCCtt              >  2:835757/1‑150 (MQ=255)
                                                                                                                                       atgatgAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTa             <  1:293894/150‑1 (MQ=255)
                                                                                                                                             aGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTAcggcgg       <  2:579187/150‑1 (MQ=255)
                                                                                                                                             aGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTAcggcgg       >  2:348318/1‑150 (MQ=255)
                                                                                                                                               cAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGAc     <  2:990791/150‑1 (MQ=255)
                                                                                                                                                  gCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGa  <  2:1042929/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TACGTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGA  >  NC_000913/2406233‑2406528

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