Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,000,368 C→T 5.5% L6L (CTG→CTA fliY ← cystine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,000,3680CT 5.5% 241.2 / 3.4 73L6L (CTG→CTAfliYcystine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base T (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

GAACGGCGGATAAGTTCCTTCCAGCCCTACCAGCAGCGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAATATATATCAGAAGAAGGCAGGCTGGAGGATTTCTGCACAA  >  NC_000913/2000224‑2000516
                                                                                                                                                |                                                                                                                                                    
gAACGGCGGATAAGTTCCTTCCAGCCCTACCAGCAGCGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGAt                                                                                                                                                  >  2:70472/1‑149 (MQ=255)
   cggcggATAAGTTCCTTCCAGCCCTACCAGCAGCGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGAtgtg                                                                                                                                               >  3:662287/1‑149 (MQ=255)
    ggcggATAAGTTCCTTCCAGCCCTACCAGCAGCGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGc                                                                                                                                              <  4:470663/149‑1 (MQ=255)
            aGTTCCTTCCAGCCCTACCAGCAGCGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCa                                                                                                                                      <  1:117567/149‑1 (MQ=255)
              ttccttccAGCCCTACCAGCAGCGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCata                                                                                                                                    >  4:804090/1‑149 (MQ=255)
              ttccttccAGCCCTACCAGCAGCGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCata                                                                                                                                    >  3:519120/1‑149 (MQ=255)
                             ccAGCAGCGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAAtgtt                                                                                                                     <  1:660897/149‑1 (MQ=255)
                                gcagcGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATgttgtt                                                                                                                  >  2:644744/1‑149 (MQ=255)
                                gcagcGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATgttgtt                                                                                                                  <  3:457978/149‑1 (MQ=255)
                                 cagcGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTa                                                                                                                 <  4:484234/149‑1 (MQ=255)
                                   gcGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATg                                                                                                               >  2:614316/1‑149 (MQ=255)
                                    cGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGt                                                                                                              <  1:614316/149‑1 (MQ=255)
                                    cGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGt                                                                                                              <  1:6548/149‑1 (MQ=255)
                                      tGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCt                                                                                                            <  1:212388/149‑1 (MQ=255)
                                      tGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCt                                                                                                            <  4:109759/149‑1 (MQ=255)
                                              tctTTAACTTTATTAAGCAGACCTCAATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCAATAATTGCCTTCCGGCCCAGTTGTGCTAATTTCTTATTCACCCCGATTGTTGTTATGTCTGTTTGGAg                                                                                                    >  4:671279/1‑149 (MQ=255)
                                                tttAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATAtt                                                                                                                                  <  3:505231/117‑1 (MQ=255)
                                                tttAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATAtt                                                                                                                                  >  4:505231/1‑117 (MQ=255)
                                                tttAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAgtg                                                                                                  >  2:198958/1‑149 (MQ=255)
                                                  tAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTa                                                                                                <  1:279915/149‑1 (MQ=255)
                                                    aCTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTaga                                                                                              <  3:804090/149‑1 (MQ=255)
                                                         attaAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAt                                                                                         >  2:65914/1‑149 (MQ=255)
                                                          ttaAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAtt                                                                                        >  2:130204/1‑149 (MQ=255)
                                                            aaGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgt                                                                                      >  4:513375/1‑149 (MQ=255)
                                                            aaGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgt                                                                                      <  1:644744/149‑1 (MQ=255)
                                                            aaGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgt                                                                                      >  3:700496/1‑149 (MQ=255)
                                                             aGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgtt                                                                                     >  3:523498/1‑149 (MQ=255)
                                                              gCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgttt                                                                                    >  2:495334/1‑149 (MQ=255)
                                                              gCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgttt                                                                                    <  4:662287/149‑1 (MQ=255)
                                                                aGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgtttgt                                                                                  <  4:246511/149‑1 (MQ=255)
                                                                   ccTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTaaa                                                                               <  4:523498/149‑1 (MQ=255)
                                                                   ccTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTaaa                                                                               <  2:90904/149‑1 (MQ=255)
                                                                        aTCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACaa                                                                          <  4:257318/149‑1 (MQ=255)
                                                                         tCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACaaa                                                                         >  3:802133/1‑149 (MQ=255)
                                                                           tGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACaaaaa                                                                       <  4:700496/149‑1 (MQ=255)
                                                                            gCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAaac                                                                      >  4:168957/1‑149 (MQ=255)
                                                                              aaaaCTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTa                                                                                                                 <  2:415926/104‑1 (MQ=255)
                                                                              aaaaCTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTa                                                                                                                 >  1:415926/1‑104 (MQ=255)
                                                                              aaaaCTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAaacaa                                                                    <  2:546227/149‑1 (MQ=255)
                                                                              aaaaCTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAaacaa                                                                    <  1:130204/149‑1 (MQ=255)
                                                                              aaaaCTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAaacaa                                                                    >  1:10888/1‑149 (MQ=255)
                                                                               aaaCTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAaacaac                                                                   <  4:193183/149‑1 (MQ=255)
                                                                                  cTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTcc                                                                >  2:121252/1‑149 (MQ=255)
                                                                                          gCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACAtc                                                        >  3:500943/1‑149 (MQ=255)
                                                                                               tGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACAtcttatt                                                   <  1:495334/149‑1 (MQ=255)
                                                                                                gCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgt                                                                                      <  3:279343/113‑1 (MQ=255)
                                                                                                gCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgt                                                                                      >  4:279343/1‑113 (MQ=255)
                                                                                                gCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACAtcttattc                                                  <  3:671279/149‑1 (MQ=255)
                                                                                                  ccGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCTAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgttt                                                                                    >  3:586660/1‑113 (MQ=255)
                                                                                                  ccGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCTAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgttt                                                                                    >  1:481578/1‑113 (MQ=255)
                                                                                                  ccGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCTAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgttt                                                                                    <  4:586660/113‑1 (MQ=255)
                                                                                                  ccGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCTAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgttt                                                                                    <  2:481578/113‑1 (MQ=255)
                                                                                                  ccGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCt                                                                                                            <  1:454712/89‑1 (MQ=255)
                                                                                                  ccGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCt                                                                                                            >  2:454712/1‑89 (MQ=255)
                                                                                                  ccGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGc                                                              <  1:317866/135‑1 (MQ=255)
                                                                                                  ccGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGc                                                              >  2:317866/1‑135 (MQ=255)
                                                                                                        ccAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACAtcttattcttat                                              >  4:70071/1‑145 (MQ=255)
                                                                                                        ccAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACAtcttattcttat                                              <  3:70071/145‑1 (MQ=255)
                                                                                                        ccAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACAtcttat                                                    >  2:561328/1‑139 (MQ=255)
                                                                                                        ccAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACAtcttat                                                    <  1:561328/139‑1 (MQ=255)
                                                                                                        ccAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTaa                                          >  4:301676/1‑149 (MQ=255)
                                                                                                        ccAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTaa                                          <  1:70472/149‑1 (MQ=255)
                                                                                                          aGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAata                                        <  4:519120/149‑1 (MQ=255)
                                                                                                           gcgcCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAatat                                       >  3:388288/1‑149 (MQ=255)
                                                                                                              ccACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAatatata                                    <  4:388288/149‑1 (MQ=255)
                                                                                                                    ccATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAATATATATCagaa                              <  1:198958/149‑1 (MQ=255)
                                                                                                                       tCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAATATATATCAgaagaa                           >  3:129014/1‑149 (MQ=255)
                                                                                                                             ccATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgttt                                                                                    <  1:74167/86‑1 (MQ=255)
                                                                                                                             ccATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCAttgttt                                                                                    >  2:74167/1‑86 (MQ=255)
                                                                                                                                  aaTGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAATATATATCAGAAGAAGGCAGGCTgga                <  3:513375/149‑1 (MQ=255)
                                                                                                                                   aTGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAATATATATCAGAAGAAGGCAGGCTggag               >  4:390160/1‑149 (MQ=255)
                                                                                                                                          aCGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAATATATATCAGAAGAAGGCAGGCTGGAGGATTTCt        >  4:446494/1‑149 (MQ=255)
                                                                                                                                             tCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAATATATATCAGAAGAAGGCAGGCTGGAGGATTTCTGca     <  4:802133/149‑1 (MQ=255)
                                                                                                                                              cccAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAATATATATCAGAAGAAGGCAGGCTGGAGGATTTCTGcac    >  1:191681/1‑149 (MQ=255)
                                                                                                                                                cAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAATATATATCAGAAGAAGGCAGGCTGGAGGATTTCTGCACaa  <  2:191681/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GAACGGCGGATAAGTTCCTTCCAGCCCTACCAGCAGCGTGCCGCGCTCTTTAACTTTATTAAGCAGACCTTCATCTGCAAAACTTTTAACGCTCATGCCCGCAACCAGCGCCACGGCCATCACACCCATCAATGCCTGACGTCCCAGATGTGCTAATTTCATATTCACCCCGAATGTTGTTATGTCTGTTTGCAGTGTAGAGCCATTGTTTGTAAACACAAAAACAACTCCGCTACATCTTATTCTTATTTAATATATATCAGAAGAAGGCAGGCTGGAGGATTTCTGCACAA  >  NC_000913/2000224‑2000516

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