Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,987,103 C→T 6.7% intergenic (+27/‑433) yqeG → / → yqeH putative transporter/putative LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,987,1030CT 6.7% 311.5 / 8.7 105intergenic (+27/‑433)yqeG/yqeHputative transporter/putative LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (45/51);  new base T (3/4);  total (50/55)
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.13e-01

Only 100 of 107 total aligned reads displayed.

CTTTTCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACT  >  NC_000913/2986954‑2987251
                                                                                                                                                     |                                                                                                                                                    
cTTTTCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGc                                                                                                                                                      >  1:140405/1‑150 (MQ=255)
attttCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTTTGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGc                                                                                                                                                      <  1:1094412/149‑1 (MQ=255)
 ttttCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCt                                                                                                                                                     <  2:733594/150‑1 (MQ=255)
   ttCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTAc                                                                                                                                                   >  1:338592/1‑150 (MQ=255)
          tGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTg                                                                                                                                            <  1:407044/150‑1 (MQ=255)
           gCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGa                                                                                                                                           <  1:1096989/150‑1 (MQ=255)
           gCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGa                                                                                                                                           >  1:584370/1‑150 (MQ=255)
            ccTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTTATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGaa                                                                                                                                          >  2:629286/1‑150 (MQ=255)
            ccTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGaa                                                                                                                                          <  2:1080896/150‑1 (MQ=255)
             cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAc                                                                                                                                         <  1:1080655/150‑1 (MQ=255)
             cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAc                                                                                                                                         <  2:1051334/150‑1 (MQ=255)
              tACGCTGTCTACGTATCTCTTCCCGGCGCTTAATCACTGGCTTTGCATCTGAAATCTGATTAGGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTTCTACAATTTTGAAcc                                                                                                                                        <  2:538731/150‑1 (MQ=255)
              tACGCTGTCAACGTATCTCATCCCGGCGCGTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAcc                                                                                                                                        <  1:1071316/150‑1 (MQ=255)
              tACGCTGTCAACGTATCTCATCCCGGCGCCTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAcc                                                                                                                                        <  1:434393/150‑1 (MQ=255)
                  cTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccatta                                                                                                                                    >  1:305549/1‑149 (MQ=255)
                   tGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattcc                                                                                                                                   <  1:1189526/150‑1 (MQ=255)
                      cAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATCGTTCAAAATTGCTACAATTTTGAAccattcaatt                                                                                                                                >  2:185285/1‑150 (MQ=255)
                       aaCGTATCTCCTCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGAACGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTAGTGGATGGTTCAAAATTGCTAAAATATTGAACCCTTCCAttt                                                                                                                               >  2:805078/1‑150 (MQ=255)
                           tATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCa                                                                                                                           <  2:1106962/150‑1 (MQ=255)
                              ctcATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCataa                                                                                                                        >  1:567466/1‑150 (MQ=255)
                                        cgcTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCAtttt                                                                                                              >  2:1272420/1‑149 (MQ=255)
                                          cTTAAACCCTGGCGTTCCATCGGAAATCTGGTTACGCTGATCGTGGGTATCCTGTGCGTATCGGTTATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGATACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtt                                                                                                            >  2:83020/1‑150 (MQ=255)
                                           ttAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtta                                                                                                           <  1:623521/150‑1 (MQ=255)
                                           ttAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtta                                                                                                           <  2:609985/150‑1 (MQ=255)
                                            tAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttat                                                                                                          >  1:402027/1‑150 (MQ=255)
                                             aaaCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTAATTTGttatt                                                                                                         >  2:1297892/1‑150 (MQ=255)
                                              aaCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttatta                                                                                                        <  1:889256/150‑1 (MQ=255)
                                               aCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGTGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTtattat                                                                                                       >  1:950489/1‑150 (MQ=255)
                                                cccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTattatt                                                                                                      >  2:1227190/1‑150 (MQ=255)
                                                cccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTattatt                                                                                                      >  1:59484/1‑150 (MQ=255)
                                                  cTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTcc                                                                                                    <  2:409339/150‑1 (MQ=255)
                                                    ggCTTTCCATCGTTAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATAAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAc                                                                                                  <  1:116228/150‑1 (MQ=255)
                                                    ggCGTTCCATCTGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCTTATCGGAAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAc                                                                                                  <  1:243807/150‑1 (MQ=255)
                                                    ggCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAc                                                                                                  <  2:399990/150‑1 (MQ=255)
                                                     gCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAcc                                                                                                 >  1:915397/1‑150 (MQ=255)
                                                          ccATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATc                                                                                            >  1:280496/1‑150 (MQ=255)
                                                           cATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCa                                                                                           >  2:718602/1‑150 (MQ=255)
                                                            aTCGGAAATCTGATTACTCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCAc                                                                                          <  1:1247656/150‑1 (MQ=255)
                                                              cGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACtt                                                                                        >  1:277577/1‑150 (MQ=255)
                                                                 aaaTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctc                                                                                     >  2:566133/1‑150 (MQ=255)
                                                                   aTCTGATTACTCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctctc                                                                                   <  1:1310183/150‑1 (MQ=255)
                                                                    tCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGATACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACACTATCACTTCTCTCg                                                                                  >  2:1193907/1‑150 (MQ=255)
                                                                      tGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTc                                                                                <  1:880489/150‑1 (MQ=255)
                                                                       gATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGGTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAAAATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCa                                                                               <  2:555535/150‑1 (MQ=255)
                                                                          tACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCatat                                                                            <  1:780910/150‑1 (MQ=255)
                                                                            cacTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATcc                                                                          <  1:532542/148‑1 (MQ=255)
                                                                            cGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATcc                                                                          >  2:802608/1‑150 (MQ=255)
                                                                               tGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAg                                                                       <  2:321965/150‑1 (MQ=255)
                                                                                gATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGt                                                                      >  2:1065230/1‑150 (MQ=255)
                                                                                  tCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtat                                                                    <  1:1053327/150‑1 (MQ=255)
                                                                                  tCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtat                                                                    <  2:685244/150‑1 (MQ=255)
                                                                                    gTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtatta                                                                  <  1:106581/150‑1 (MQ=255)
                                                                                      gggTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATc                                                                >  2:769732/1‑150 (MQ=255)
                                                                                        gTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCag                                                              <  2:265185/150‑1 (MQ=255)
                                                                                        gTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCag                                                              >  2:932335/1‑150 (MQ=255)
                                                                                          aTCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtt                                                            <  2:915397/150‑1 (MQ=255)
                                                                                            ccTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttat                                                          >  2:495072/1‑149 (MQ=255)
                                                                                              tgtgAGTATCGGTAATTTTCTTTAGCTAATATCAGTTGTGGTTGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagatagtt                                                        <  2:1155054/150‑1 (MQ=255)
                                                                                               gtgCGTATAGGTAATGTTCTTTAGCTTATATCAGTTGTTGATTGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCAATTCTCTCGTCATATCCTAGTATTATCAgttagtta                                                       <  1:1068828/150‑1 (MQ=255)
                                                                                                 gCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaaa                                                     <  1:1242024/150‑1 (MQ=255)
                                                                                                  cGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAac                                                    <  2:584370/150‑1 (MQ=255)
                                                                                                   gTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAaca                                                   >  2:1147335/1‑150 (MQ=255)
                                                                                                   gTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAaca                                                   <  2:1231332/150‑1 (MQ=255)
                                                                                                     aTCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCAGTCCATTTTTCATAATATTCATTTGTTATTATTACACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAacaca                                                 >  2:344848/1‑150 (MQ=255)
                                                                                                       cGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAAcacaca                                               <  2:3838/150‑1 (MQ=255)
                                                                                                       cGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAAcacaca                                               >  1:925322/1‑150 (MQ=255)
                                                                                                             tGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAgaga                                         >  1:1099744/1‑150 (MQ=255)
                                                                                                              gTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCATAGTATTATCAGTTAGTTAAACACACAAAGAGAt                                        >  2:927051/1‑150 (MQ=255)
                                                                                                               ttCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATc                                       <  2:983645/150‑1 (MQ=255)
                                                                                                                 cTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTg                                     <  2:338592/150‑1 (MQ=255)
                                                                                                                 atttATCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTg                                     <  1:426360/149‑1 (MQ=255)
                                                                                                                   ttAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGaa                                   <  2:140405/150‑1 (MQ=255)
                                                                                                                    tAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAAc                                  <  2:363853/150‑1 (MQ=255)
                                                                                                                     aGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACa                                 >  2:125361/1‑150 (MQ=255)
                                                                                                                      gCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACAt                                >  1:176734/1‑150 (MQ=255)
                                                                                                                       cTAATATCAGTTGTGGATGGTTCAAAATTTCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCt                              <  1:984323/150‑1 (MQ=255)
                                                                                                                          atatCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTTTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTaa                           <  2:1306536/150‑1 (MQ=255)
                                                                                                                           tatCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTaa                           <  2:402027/150‑1 (MQ=255)
                                                                                                                             tCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGCACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTAGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGAACTGAACAACTTAAAt                         >  2:801731/1‑150 (MQ=255)
                                                                                                                               aGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATag                       >  1:692046/1‑149 (MQ=255)
                                                                                                                                gTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAg                      >  1:1174620/1‑150 (MQ=255)
                                                                                                                                gTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAg                      <  2:1268163/150‑1 (MQ=255)
                                                                                                                                  tgtgGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGAc                    <  2:1069942/150‑1 (MQ=255)
                                                                                                                                   gtgGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACACTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACt                   >  2:927775/1‑150 (MQ=255)
                                                                                                                                    tgGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTa                  >  2:1307565/1‑150 (MQ=255)
                                                                                                                                         ggTTCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTaa             >  1:831271/1‑150 (MQ=255)
                                                                                                                                          gTTCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg            >  1:1057531/1‑150 (MQ=255)
                                                                                                                                          gTTCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg            >  1:118139/1‑150 (MQ=255)
                                                                                                                                          gTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTaat            >  2:201821/1‑149 (MQ=255)
                                                                                                                                          gTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg            >  2:20494/1‑150 (MQ=255)
                                                                                                                                          gTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg            >  2:967533/1‑150 (MQ=255)
                                                                                                                                          gTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGACATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg            >  2:1134888/1‑150 (MQ=255)
                                                                                                                                           ttCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa           <  1:964649/150‑1 (MQ=255)
                                                                                                                                            tCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGAt          <  1:315528/150‑1 (MQ=255)
                                                                                                                                            tCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCAACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACTCAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGAt          <  2:33501/150‑1 (MQ=255)
                                                                                                                                             cAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATg         <  1:688232/150‑1 (MQ=255)
                                                                                                                                             cAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTAATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGACCATCTTAAATAAGCATATTTAATATg         >  2:531387/1‑150 (MQ=255)
                                                                                                                                                aaTTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAGTCTTATTTAAGATgttc      <  2:850293/150‑4 (MQ=255)
                                                                                                                                                 aTTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATa     <  2:280496/150‑1 (MQ=255)
                                                                                                                                                    gCTACAATTTTGAACCATTCCATTTTTCATAATATTAATTTGTTATTATTCCACCCTATACCTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTTAACACACATAGCGATCTGAACATCTTAAATATGCCTATTTAAGATGGCTAACt  >  2:185279/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTTTCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACT  >  NC_000913/2986954‑2987251

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