Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,509,243 C→T 8.2% intergenic (‑179/‑43) yncJ ← / → hicA uncharacterized protein/mRNA interferase toxin of the HicAB toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,509,2430CT 8.2% 230.4 / 10.0 73intergenic (‑179/‑43)yncJ/hicAuncharacterized protein/mRNA interferase toxin of the HicAB toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (34/33);  new base T (3/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-01

TCGTTTCTTGAGGTTATTATTCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTC  >  NC_000913/1509105‑1509387
                                                                                                                                          |                                                                                                                                                
tCGTTTCTTGAGGTTATTATTCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTa                                                                                                                                         <  3:199347/148‑1 (MQ=255)
   tttCTTGAGGTTATTATTCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGt                                                                                                                                            <  4:716477/142‑1 (MQ=255)
   tttCTTGAGGTTATTATTCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGt                                                                                                                                            >  3:716477/1‑142 (MQ=255)
      cTTGAGGTTATTATTCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttgt                                                                                                                                   <  4:691309/148‑1 (MQ=255)
          aGGTTATTATTCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTa                                                                                                                                         >  4:1098461/1‑138 (MQ=255)
          aGGTTATTATTCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTa                                                                                                                                         <  3:1098461/138‑1 (MQ=255)
             ttattaTTCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttgt                                                                                                                                   >  3:628229/1‑141 (MQ=255)
             tcattattCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttgt                                                                                                                                   <  4:628229/139‑1 (MQ=255)
                 tattCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTTCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACTCACGGAGGTAAACACTAATGTTTACTttgt                                                                                                                                   <  4:699963/137‑1 (MQ=255)
                 tattCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTCCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttgt                                                                                                                                   >  3:699963/1‑137 (MQ=255)
                    tCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttgt                                                                                                                                   >  1:157019/1‑134 (MQ=255)
                    tCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttgt                                                                                                                                   <  2:157019/134‑1 (MQ=255)
                    tCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTc                                                                                                                         >  2:119395/1‑144 (MQ=255)
                    tCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTc                                                                                                                         <  1:119395/144‑1 (MQ=255)
                              aaaTTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTaa                                                                                                                                               >  3:761725/1‑112 (MQ=255)
                              aaaTTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTaa                                                                                                                                               >  3:696894/1‑112 (MQ=255)
                              aaaTTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTaa                                                                                                                                               <  4:696894/112‑1 (MQ=255)
                              aaaTTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTaa                                                                                                                                               <  4:761725/112‑1 (MQ=255)
                                         aaGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttgttg                                                                                                                                 >  2:937291/1‑115 (MQ=255)
                                         aaGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttgttg                                                                                                                                 <  1:937291/115‑1 (MQ=255)
                                           gAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGAtg                                                                                                       >  1:1052592/1‑139 (MQ=255)
                                           gAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGAtg                                                                                                       <  2:1052592/139‑1 (MQ=255)
                                                 cTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTc                                                                                                                         >  4:662053/1‑115 (MQ=255)
                                                 cTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTc                                                                                                                         <  3:662053/115‑1 (MQ=255)
                                                   ggAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGa                                                                                                              <  2:1084792/124‑1 (MQ=255)
                                                   ggAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGa                                                                                                              >  1:1084792/1‑124 (MQ=255)
                                                      aTCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGaa                                                                                                   >  1:937144/1‑132 (MQ=255)
                                                      aTCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGaa                                                                                                   <  2:937144/132‑1 (MQ=255)
                                                                   aTTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAg                                                                                             <  1:896188/125‑1 (MQ=255)
                                                                   aTTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAg                                                                                             >  2:896188/1‑125 (MQ=255)
                                                                      ttgcattgcaTTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttg                                                                                                                                    >  4:524783/1‑83 (MQ=255)
                                                                      ttgcattgcaTTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTttg                                                                                                                                    <  3:524783/83‑1 (MQ=255)
                                                                           ttgcattCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACATTAGTGTTTAcctc                                                                                                                                     >  2:785751/1‑74 (MQ=39)
                                                                           ttgcattCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACATTAGTGTTTAcctc                                                                                                                                     <  3:314611/77‑4 (MQ=39)
                                                                           ttgcattCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACATTAGTGTTTAcctc                                                                                                                                     >  4:314611/1‑74 (MQ=39)
                                                                           ttgcattCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACATTAGTGTTTAcctc                                                                                                                                     <  1:785751/77‑4 (MQ=39)
                                                                           ttgcattCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGa                                                                                                                                <  2:537038/82‑1 (MQ=255)
                                                                           ttgcattCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGa                                                                                                                                >  1:537038/1‑82 (MQ=255)
                                                                                    gCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTa                                                     >  2:32696/1‑148 (MQ=255)
                                                                                        ttGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGaa                                                                                                   >  2:538239/1‑98 (MQ=255)
                                                                                        ttGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGaa                                                                                                   <  1:538239/98‑1 (MQ=255)
                                                                                         tGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAAtctc                                                                   >  1:361239/1‑129 (MQ=255)
                                                                                         tGCGGCGATTTAGTGCTATTGTCGAACAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCATGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAAtctc                                                                   <  2:361239/129‑1 (MQ=255)
                                                                                              cgATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGa                                                 <  4:570624/142‑1 (MQ=255)
                                                                                              cgATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGa                                                 >  3:570624/1‑142 (MQ=255)
                                                                                                  ttAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGt                                                                                        <  1:884288/99‑1 (MQ=255)
                                                                                                  ttAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGt                                                                                        >  2:884288/1‑99 (MQ=255)
                                                                                                      tGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCAt                                   <  1:32696/148‑1 (MQ=255)
                                                                                                       gCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGa                                                                                                              >  2:338064/1‑72 (MQ=255)
                                                                                                       gCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGa                                                                                                              <  1:338064/72‑1 (MQ=255)
                                                                                                       gCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAAtct                                                                    >  1:169242/1‑114 (MQ=255)
                                                                                                       gCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAAtct                                                                    <  2:169242/114‑1 (MQ=255)
                                                                                                          aTTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGa                               >  4:127336/1‑148 (MQ=255)
                                                                                                                       ttACACACGGAGGTAAACACTAATGTTTATTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCa                         >  4:907396/1‑141 (MQ=255)
                                                                                                                       ttACACACGGAGGTAAACACTAATGTTTATTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCa                         <  3:907396/141‑1 (MQ=255)
                                                                                                                       ttACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACaa                                                                                               <  2:457433/71‑1 (MQ=255)
                                                                                                                       ttACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACaa                                                                                               >  1:457433/1‑71 (MQ=255)
                                                                                                                       ttACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGaa                              <  3:413922/136‑1 (MQ=255)
                                                                                                                       ttACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGaa                              >  4:413922/1‑136 (MQ=255)
                                                                                                                              cGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGa                                                                                   >  1:468280/1‑76 (MQ=255)
                                                                                                                              cGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGa                                                                                   <  2:468280/76‑1 (MQ=255)
                                                                                                                              cGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGt                                                            >  1:708543/1‑99 (MQ=255)
                                                                                                                              cGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGt                                                            <  2:708543/99‑1 (MQ=255)
                                                                                                                              cGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATg                <  4:143670/143‑1 (MQ=255)
                                                                                                                              cGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATg                >  3:143670/1‑143 (MQ=255)
                                                                                                                                gAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGAAACCATTTGAAACTCAGGTTTCATGGGAGgcg         >  4:684248/1‑148 (MQ=255)
                                                                                                                                aaagTAAACATTAGTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATg                >  3:976947/4‑141 (MQ=255)
                                                                                                                                aaactAAACATTAGTGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATg                <  4:976947/137‑1 (MQ=255)
                                                                                                                                 aGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCAtt                                  >  3:706939/1‑122 (MQ=255)
                                                                                                                                 aGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCAtt                                  <  4:706939/122‑1 (MQ=255)
                                                                                                                                 aGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATg                <  4:134768/140‑1 (MQ=255)
                                                                                                                                 aGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATg                >  3:134768/1‑140 (MQ=255)
                                                                                                                                  ggTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTGTCATGGGTGacgcg       >  4:504702/1‑147 (MQ=255)
                                                                                                                                       acacTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATg                <  4:635765/134‑1 (MQ=255)
                                                                                                                                       acacTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATg                >  3:635765/1‑134 (MQ=255)
                                                                                                                                       acacTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTc  >  3:1095843/1‑148 (MQ=255)
                                                                                                                                          cTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAAtct                                                                    >  4:335984/1‑79 (MQ=255)
                                                                                                                                          cTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAAtct                                                                    <  3:335984/79‑1 (MQ=255)
                                                                                                                                          cTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAAt                                               <  1:1012162/100‑1 (MQ=255)
                                                                                                                                          cTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAAt                                               >  2:1012162/1‑100 (MQ=255)
                                                                                                                                          cTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCa       >  1:372704/1‑140 (MQ=255)
                                                                                                                                          cTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCa       <  2:372704/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                
TCGTTTCTTGAGGTTATTATTCAGTTTTGCAAATTAGCGCAAAGAAATTCTGGAATCTTCCTTCCTGATTTTGCATTGCATTCTGCCGTTGCGGCGATTTAGTGCTATTTTCGAGCAAATTACACACGGAGGTAAACACTAATGTTTACTTTGTTGATATACTCAGCGGCAGGGAGGCGATGTGAAACAAAGCGAGTTCAGACGTTGGCTCGAATCTCAGGGCGTCGATGTAGCGAATGGCAGCAACCATTTGAAACTCAGGTTTCATGGGAGGCGCAGTGTC  >  NC_000913/1509105‑1509387

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