Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,196,219 A→C 6.6% intergenic (‑120/+113) thiC ← / ← rsd phosphomethylpyrimidine synthase/regulator of sigma D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,196,2190AC 6.6% 221.1 / 2.2 87intergenic (‑120/+113)thiC/rsdphosphomethylpyrimidine synthase/regulator of sigma D
Reads supporting (aligned to +/- strand):  ref base A (44/35);  new base C (3/4);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCG  >  NC_000913/4196078‑4196359
                                                                                                                                             |                                                                                                                                            
gggtgtattttGGTGCAGACATGGCGCATTCCAAAAGGTTAAGGACGTGGCTTGTCAGGCGACGGATGAAGCAAGAGACGATGGCCGCGGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGAt                                                                                                                                             <  2:286856/136‑1 (MQ=255)
gggTCAGTTTTTTTGCAGTCATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGAt                                                                                                                                             <  1:21722/143‑1 (MQ=255)
gggTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTACGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGAt                                                                                                                                             <  1:231921/143‑1 (MQ=255)
 ggTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCCGATTAACTCTTGTTCCCTTCGCAGGCATTAGCCTGATc                                                                                                                                            >  2:744247/1‑143 (MQ=255)
  gTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGGTAAGGACGTGGCTTGTCAGGCGACGGAGGAAGCAAGTGATGTGCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCa                                                                                                                                           <  2:139121/143‑1 (MQ=255)
    cAGTTTTGTTGCAGACATAGCTCATTCCAAAAGGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGACCGTGTAGCGGTGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAgg                                                                                                                                         >  1:615899/1‑143 (MQ=255)
      gTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATGAGGtt                                                                                                                                       <  1:697065/143‑1 (MQ=255)
       ttttGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCCGGTATTAGCCTGAACAGGTTc                                                                                                                                      >  2:51269/1‑143 (MQ=255)
       ttttGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCGTGTTACCTTCGCAGGTATTAGCCTGATCAGGTTc                                                                                                                                      >  2:322274/1‑143 (MQ=255)
              gCAGACATAGCTCATTCCACCAAGTTCAGGACGTGCCTCGTGAGACGACGGATGAAGCAAGAGACGATCGCAGCTGGGGGCAAGAGACAGCGGATTAACTCTTGTTCGCTTCGCAGGGATTAGCCTGAGCAGGCTCCGCgggt                                                                                                                               >  2:214083/1‑141 (MQ=255)
                    aTAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGaa                                                                                                                         <  1:640236/143‑1 (MQ=255)
                        cTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATGGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCGGGTATTAGCCTGAGCAGGTTCCGCGGATCCCgaagaac                                                                                                                     >  1:256259/1‑142 (MQ=255)
                         tCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAAc                                                                                                                    >  2:116193/1‑143 (MQ=255)
                          cATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCACGAGACGATCGCCGCAGGGGCGAGGCGATAGCAAATTAACGCTTGCTCCCTTCGCGGGCAGTAGCCCGATCAGGTTCCGCGGCTACCGACAAAACg                                                                                                                   >  2:72860/1‑143 (MQ=255)
                          cATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACg                                                                                                                   <  1:204048/143‑1 (MQ=255)
                           aTTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACgg                                                                                                                  >  1:216395/1‑143 (MQ=255)
                            ttCCAAATAGTTAAGGAGTTGGCTGGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGTTTCCGCGGATCCCGAATAAACGGt                                                                                                                 <  2:514095/143‑1 (MQ=255)
                            tgccAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGAGGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGt                                                                                                                 <  2:267142/141‑1 (MQ=255)
                             tCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCGGATCAGGTTCCGCGGATCCCGAATAAACGGtc                                                                                                                >  1:669919/1‑143 (MQ=255)
                               ccagaaGTTGAGGGCGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtctc                                                                                                              <  2:749837/139‑1 (MQ=255)
                                aaaaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCa                                                                                                             <  1:471855/143‑1 (MQ=255)
                                  aaaGGTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTGCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGc                                                                                                           <  2:134955/143‑1 (MQ=255)
                                    aGTTAAGGGCGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCa                                                                                                         <  2:263388/143‑1 (MQ=255)
                                    aGTTAAGGCCGGGGCGTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCa                                                                                                         <  2:26541/143‑1 (MQ=255)
                                     ggtAAGGACCTGGCTTGTCGCGCGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACCCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAGACGGTCTCAGCCAg                                                                                                        <  2:341597/141‑1 (MQ=255)
                                         aGGACGTCGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATCATCTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGTATCCCGAATAAACGGTCTCAGCCAGTTaa                                                                                                    <  2:782827/143‑1 (MQ=255)
                                           gACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGGTGCGATAGCAGAGGAACTCTTGTTCCCTTCGCAGGGATTAGCCTGATCAGGTGCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTACgg                                                                                                  >  1:119713/1‑143 (MQ=255)
                                                      tCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACaa                                                                                       >  1:539299/1‑143 (MQ=255)
                                                      tCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGACACTCCGACaa                                                                                       >  2:33224/1‑143 (MQ=255)
                                                        aGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGa                                                                                     >  2:291636/1‑143 (MQ=255)
                                                           cggcgGATGAAGCAAGAGACCAGCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAAt                                                                                  <  2:674856/140‑1 (MQ=255)
                                                            gacgGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATACCGAATAAACGGTCTCAGCCAGTTAAGGCACACCGACAAGaattc                                                                                 >  2:497437/1‑140 (MQ=255)
                                                              gggATGAAGCAATGGACGGTCGCCGCATGGGCGATGCGATAGCAGATTAACTCTTGCTCCCTTCGCAGCTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAg                                                                               <  2:773239/142‑1 (MQ=255)
                                                              cgGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAg                                                                               >  1:294069/1‑143 (MQ=255)
                                                              cgGATCAAGCAAGAGCCGATCGCCGCAGGGGCGCTGCGATAGCAGATTAACCCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAg                                                                               <  1:716752/143‑1 (MQ=255)
                                                               ggATGAAGCAAGAGACGCTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGc                                                                              <  2:36848/143‑1 (MQ=255)
                                                               ggATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATCAGCCTGCTACGGTTCCACGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGc                                                                              <  1:801102/143‑1 (MQ=255)
                                                                gATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGcc                                                                             >  1:558352/1‑143 (MQ=255)
                                                                    aaGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGc                                                                         <  1:844556/143‑1 (MQ=255)
                                                                    aaGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTACATTAGCAGGTATTAGCCTGATCAGATTACGCGGGTCCAGAGTTAACTGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGc                                                                         >  1:312526/1‑143 (MQ=255)
                                                                     aGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCa                                                                        >  2:577629/1‑143 (MQ=255)
                                                                      gCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCaa                                                                       >  1:770612/1‑143 (MQ=255)
                                                                       cAAGAGACGATCGCCGCAGGGGCGATGCGATCGCACATTAACTCTTGTTCCCTTCGCAGGTATTGGCCTGATCCGGTTCCGCGGATCCCGAATAAACGGTCTCCGCCAGGTAAGGCACTCCGGCACGCAATCAGCCCCCCacc                                                                      >  2:253843/1‑141 (MQ=255)
                                                                        aggagaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAAt                                                                     <  2:538849/141‑1 (MQ=255)
                                                                           tggcGATCGCCGCAGGGGCGATGCGATAGCTGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGACCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATggg                                                                  <  2:753632/140‑1 (MQ=255)
                                                                           agaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATggg                                                                  >  1:18071/1‑143 (MQ=255)
                                                                            gaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAACTAGCCCCGCAAATgggg                                                                 >  2:108678/1‑143 (MQ=255)
                                                                            gaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTAAGCCAGTTAAGGCACTGCGACAAGAAATTAGCCCAGCAAAAgggg                                                                 >  1:341720/1‑143 (MQ=255)
                                                                             aCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTCTTAGCCTGATCAGGTTCAGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGc                                                                >  2:100341/1‑143 (MQ=255)
                                                                               gATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAt                                                              >  1:540536/1‑143 (MQ=255)
                                                                                aTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCCACGGGGGCCtt                                                             >  2:480278/1‑143 (MQ=255)
                                                                                    ctgcAGGGGGGATGCGATAGCGCAGGACCTCTTGTCCCCTTCGCAGGTACTAGCGTGTGCAGGTTGCGCGGATCCCGAAGACGCGGTCTCAGCCAGTGAAGGCACTCCGACAAGAAAGTAGCCCCGCAAATGGGGCATTGAAt                                                         <  2:328722/141‑1 (MQ=255)
                                                                                    ccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAAt                                                         >  1:605936/1‑143 (MQ=255)
                                                                                    ccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAAt                                                         <  1:307713/143‑1 (MQ=255)
                                                                                      gcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATCGAAtca                                                       >  2:557903/1‑141 (MQ=255)
                                                                                       gcggggCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGCTCCCCAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTa                                                      <  1:284395/141‑1 (MQ=255)
                                                                                               atgcgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATGAACGGTCTCAGCCAGTTAAGGCACTCCGACCAGAAATTAGCCCCGCACATGGGGCATTGTGCCTAAATTAcgc                                              >  1:308934/1‑143 (MQ=255)
                                                                                                    gtAGGAGATTAACTCTTGCTGCCTTCCAAGGTATCAGCCTGATCAGGTGCCGCGGAGCCCGAATAAACTGTCCCAGCCAGTTAAGGCACGCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTaa                                         <  2:79944/142‑1 (MQ=255)
                                                                                                     tAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAACGAATGTGAATTACGCGtgagc                                        >  1:792676/1‑139 (MQ=255)
                                                                                                      aGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACa                                       <  1:639417/143‑1 (MQ=255)
                                                                                                        cAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAgc                                     <  2:435458/143‑1 (MQ=255)
                                                                                                         aGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCATCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAgcg                                    <  1:14804/143‑1 (MQ=255)
                                                                                                          gATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAGGAAATTAGCCCCGCAAAGGGGGCATTGAGTGTAAAGTACGCGTGACCGgcgc                                   >  1:33852/1‑143 (MQ=255)
                                                                                                          gATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGATTGTCAATTACGCGTTATCCgcgc                                   >  1:597803/1‑143 (MQ=255)
                                                                                                           aTTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCa                                  <  2:518768/143‑1 (MQ=255)
                                                                                                             tAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCCTAGAGGGGAAATTCGACGATAACAGCGCAGa                                >  2:590780/1‑143 (MQ=255)
                                                                                                              aaCTCTTGTTCCCTTCGCCGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTGAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACCGCGCAgca                               >  2:19678/1‑141 (MQ=255)
                                                                                                              aaCTCTTGCGCCCGTCGCAGGTATGAGCCTGCTCAGGGTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGaa                               <  2:237736/143‑1 (MQ=255)
                                                                                                                   ttGTTCCCTTCGCAGGCATTAGCCTGATCAGGTTCCGCCGATCCCGAATAAACGCTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGCGCATTCAATGTAAATTACGCCTTAACAGCGCAGAACTCaa                          <  2:429317/143‑1 (MQ=255)
                                                                                                                     gTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGc                        >  1:467729/1‑143 (MQ=255)
                                                                                                                      ttCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGTACTCCGACAAGAAATTAGCCCCGCACATGGGGAACTGAATGTCAATTACGCGATAACAGCGCAGAACTCAAGCa                       >  2:290262/1‑143 (MQ=255)
                                                                                                                      ttCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAACGGGGCACTGAATGCAAATTACGCGCTAACAGCGCAGCACACAAGCa                       >  2:226519/1‑143 (MQ=255)
                                                                                                                       tCCCTTCGCAGGTATTAGCCTGATCAGGTGCCGCGGATCCCGAATAAACGGGCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAg                      <  1:682306/143‑1 (MQ=255)
                                                                                                                        cccTTCGCAGGTAGTAGCCTGATCAGGTTCCGCGGATCCCGAAGAAACGGTCTCAGCCAGGTAAGGCACGCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAgg                     <  1:344908/143‑1 (MQ=255)
                                                                                                                           ttCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCCGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAGTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATg                  >  2:166389/1‑143 (MQ=255)
                                                                                                                            tCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGt                 >  1:490572/1‑143 (MQ=255)
                                                                                                                                aGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATGAGCCCCGCAAATGGGGCATTGAATGTACATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTgc             >  2:764374/1‑142 (MQ=255)
                                                                                                                                 ggTGTTTGCCTGCTCAGGTTCCGCGGATCCCGAATTAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATGAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGAc            <  2:241558/143‑1 (MQ=255)
                                                                                                                                  gTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGAcg           >  2:780483/1‑143 (MQ=255)
                                                                                                                                   tATTGGCCTGATCAGGTGCAGCGGATCCCGAATAAACGGTCGCAGCCAGCTAAGGCACTCCGACAAGAAATCAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGAcgc          <  1:462480/143‑1 (MQ=255)
                                                                                                                                     ttAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcgg        <  2:382641/143‑1 (MQ=255)
                                                                                                                                       aGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGACCTCAAGCAGGATGTTTGACgcgggc      >  1:331976/1‑143 (MQ=255)
                                                                                                                                       aGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCCCAGAACTCAAGCAGGATGTTTGACgcgggc      >  2:368635/1‑143 (MQ=255)
                                                                                                                                        gCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcgggcg     >  2:527328/1‑143 (MQ=255)
                                                                                                                                          cTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGggcggc   >  1:588279/1‑143 (MQ=255)
                                                                                                                                          cTGATCAGGTGCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGggcggc   <  1:590780/143‑1 (MQ=255)
                                                                                                                                           tGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAAGAGCGCAGAACTGAAGCAGGATGTGTGGCGCGGgcggcg  >  2:844429/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                            
GGGTCAGTTTTGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCG  >  NC_000913/4196078‑4196359

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