Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,948,179 T→A 10.4% intergenic (‑28/+1) yebC ← / ← nudB putative transcriptional regulator YebC/dihydroneopterin triphosphate diphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,948,1790TA10.4% 239.2 / 13.0 85intergenic (‑28/+1)yebC/nudBputative transcriptional regulator YebC/dihydroneopterin triphosphate diphosphatase
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base A (6/3);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

CCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAA  >  NC_000913/1948030‑1948326
                                                                                                                                                     |                                                                                                                                                   
ccACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGAGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGt                                                                                                                                                     >  2:109951/1‑150 (MQ=255)
 cACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGtt                                                                                                                                                    <  2:1402787/150‑1 (MQ=255)
   ccGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCa                                                                                                                                                  <  2:1181040/150‑1 (MQ=255)
     gCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAgg                                                                                                                                                <  1:1185809/150‑1 (MQ=255)
        cAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAg                                                                                                                                             <  2:409176/150‑1 (MQ=255)
         aGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGc                                                                                                                                            >  2:1370730/1‑150 (MQ=255)
         aGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAgg                                                                                                                                            <  2:1122508/150‑2 (MQ=255)
          gCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCg                                                                                                                                           >  1:1254742/1‑150 (MQ=255)
          gCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCg                                                                                                                                           <  2:18321/150‑1 (MQ=255)
          gCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCg                                                                                                                                           <  2:1130965/150‑1 (MQ=255)
           cTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGt                                                                                                                                          <  1:495870/150‑1 (MQ=255)
             tAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTa                                                                                                                                        <  1:48439/150‑1 (MQ=255)
               gCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAt                                                                                                                                      >  1:972139/1‑150 (MQ=255)
               gCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAt                                                                                                                                      <  1:406456/150‑1 (MQ=255)
                ccGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAtt                                                                                                                                     <  1:802844/150‑1 (MQ=255)
                  gcgGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTAc                                                                                                                                   >  2:1113022/1‑150 (MQ=255)
                     gTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACaaa                                                                                                                                <  2:30730/150‑1 (MQ=255)
                       tACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACt                                                                                                                              >  1:1139415/1‑150 (MQ=255)
                          caccTCACCACTGATCTTAGTGAATATTTTACCGCGCCTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGtt                                                                                                                           <  2:305638/147‑1 (MQ=255)
                          cAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGtt                                                                                                                           >  1:435752/1‑150 (MQ=255)
                          cAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGtt                                                                                                                           <  2:1317081/150‑1 (MQ=255)
                           aGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTc                                                                                                                          <  1:865353/150‑1 (MQ=255)
                           aGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCATGCCTGTACAGGCAGAGTTAATTACAAAATGATc                                                                                                                          >  2:1322436/1‑150 (MQ=255)
                            gCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCa                                                                                                                         >  1:339828/1‑150 (MQ=255)
                                aCGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCg                                                                                                                     >  2:1241535/1‑150 (MQ=255)
                                 cGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGc                                                                                                                    >  1:899466/1‑150 (MQ=255)
                                   aaTGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCt                                                                                                                  >  1:30778/1‑150 (MQ=255)
                                         tttAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGt                                                                                                            <  1:194623/150‑1 (MQ=255)
                                          ttAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTTGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTAGAAACTGTTCAATCGCCTGCCGGtt                                                                                                           <  2:508170/150‑1 (MQ=255)
                                           tAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCAGTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTg                                                                                                          >  2:419025/1‑150 (MQ=255)
                                            aGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGc                                                                                                         >  2:1227874/1‑150 (MQ=255)
                                               gAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTcc                                                                                                      >  2:1071446/1‑150 (MQ=255)
                                               gAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTcc                                                                                                      >  2:246820/1‑150 (MQ=255)
                                               gAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTcc                                                                                                      <  1:1113022/150‑1 (MQ=255)
                                               gAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAg                                                                                                                                             <  1:493357/111‑1 (MQ=255)
                                               gAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAg                                                                                                                                             >  2:493357/1‑111 (MQ=255)
                                                  caTTTTCCCGCGCTTAGACTCCTGCGCAGCTTTACGATGTCTTGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAgg                                                                                                   <  2:721124/149‑1 (MQ=255)
                                                    ttttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGAc                                                                                                 >  1:729898/1‑150 (MQ=255)
                                                        aCCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAg                                                                                             <  1:1284276/150‑1 (MQ=255)
                                                          cgcgCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTg                                                                                           <  1:260984/150‑1 (MQ=255)
                                                           gcgcTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGa                                                                                          <  2:1344013/150‑1 (MQ=255)
                                                              cTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAgcg                                                                                       <  2:1139415/150‑1 (MQ=255)
                                                              cTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAgcg                                                                                       <  1:560900/150‑1 (MQ=255)
                                                                tAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGcgcc                                                                                     >  2:618197/1‑150 (MQ=255)
                                                                       cTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcag                                                                              >  2:667520/1‑150 (MQ=255)
                                                                       cTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATATACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGGTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcag                                                                              <  1:1039260/150‑1 (MQ=255)
                                                                        tGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcagc                                                                             >  1:1425115/1‑150 (MQ=255)
                                                                           gcAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAgg                                                                          >  1:806158/1‑150 (MQ=255)
                                                                              gCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGgcgc                                                                       <  2:983414/150‑1 (MQ=255)
                                                                              gCTTTACGATGTCTGGTGTGGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGgcgc                                                                       <  1:520648/150‑1 (MQ=255)
                                                                                  tACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTGCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCa                                                                   <  1:1103552/150‑1 (MQ=255)
                                                                                      aTGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGcc                                                               <  2:1129489/150‑1 (MQ=255)
                                                                                          cTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACtt                                                           >  1:202105/1‑150 (MQ=255)
                                                                                                ttGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGc                                                     <  1:1227874/150‑1 (MQ=255)
                                                                                                 tGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGcc                                                    <  1:481892/150‑1 (MQ=255)
                                                                                                  ggCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTCTAAGCCa                                                   >  2:588316/1‑150 (MQ=255)
                                                                                                   gCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAg                                                  >  1:1156258/1‑150 (MQ=255)
                                                                                                    cccATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGa                                                 <  2:72328/150‑1 (MQ=255)
                                                                                                    cccATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGa                                                 >  2:1383165/1‑150 (MQ=255)
                                                                                                    cccATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGa                                                 >  2:41506/1‑150 (MQ=255)
                                                                                                         ttACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTc                                            <  2:899466/150‑1 (MQ=255)
                                                                                                           aCTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAg                                          <  1:618197/150‑1 (MQ=255)
                                                                                                            cTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGt                                         >  2:63806/1‑150 (MQ=255)
                                                                                                              aTGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGa                                       >  2:468312/1‑150 (MQ=255)
                                                                                                                 aCCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGaaaa                                    >  1:836905/1‑150 (MQ=255)
                                                                                                                 aCCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAAACGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGaaaa                                    >  2:1432904/1‑150 (MQ=255)
                                                                                                                         tAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGc                            >  2:966009/1‑150 (MQ=255)
                                                                                                                         tAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGc                            <  2:435752/150‑1 (MQ=255)
                                                                                                                         tAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGc                            <  1:667520/150‑1 (MQ=255)
                                                                                                                           aaaaTATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCg                          >  1:1074449/1‑150 (MQ=255)
                                                                                                                                 tctcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTg                    <  1:419025/150‑1 (MQ=255)
                                                                                                                                  ctcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGc                   >  2:469658/1‑150 (MQ=255)
                                                                                                                                   tcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCg                  >  1:1630109/1‑150 (MQ=255)
                                                                                                                                   tcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCg                  <  1:1241535/150‑1 (MQ=255)
                                                                                                                                   tcCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCg                  >  1:189099/1‑150 (MQ=255)
                                                                                                                                       attaacgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAg              >  1:1536523/8‑150 (MQ=255)
                                                                                                                                       attaacgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAg              >  1:1015133/8‑150 (MQ=255)
                                                                                                                                        ttaacgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAgc             >  1:1612592/7‑150 (MQ=255)
                                                                                                                                         tAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAgcg            >  1:336446/1‑150 (MQ=255)
                                                                                                                                           gCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCa          >  2:1083640/1‑150 (MQ=255)
                                                                                                                                           gCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCa          <  2:1156258/150‑1 (MQ=255)
                                                                                                                                           gCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCa          >  1:1180308/1‑150 (MQ=255)
                                                                                                                                           acgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGc                         <  1:309599/132‑1 (MQ=255)
                                                                                                                                           acgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGc                         >  2:309599/4‑135 (MQ=255)
                                                                                                                                           acgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCa          >  2:1354571/4‑150 (MQ=255)
                                                                                                                                              ctgcctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGa       <  2:1206965/150‑1 (MQ=255)
                                                                                                                                                   tgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGaa  >  1:573439/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAA  >  NC_000913/1948030‑1948326

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