Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,024,286 C→T 9.1% W21* (TGG→TGA)  yccF ← DUF307 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,024,2860CT 9.1% 197.2 / 7.6 66W21* (TGG→TGA) yccFDUF307 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base T (3/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

AGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTTCTCTGACGGTT  >  NC_000913/1024140‑1024431
                                                                                                                                                  |                                                                                                                                                 
agcaCATTTTTGCCGGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCATCACAATACTGACCAGAGTCGCCAAcagcc                                                                                                                                                  <  8:222351/148‑1 (MQ=255)
   acaTTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAAcagccagc                                                                                                                                               >  3:200793/1‑148 (MQ=255)
   acaTTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAAcagccagc                                                                                                                                               >  7:41134/1‑148 (MQ=255)
        tttGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCa                                                                                                                                            >  7:113386/1‑146 (MQ=255)
           gCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAgt                                                                                                                                        <  5:73685/147‑1 (MQ=255)
                    gttcagttcaTCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAAt                                                                                                                              <  8:8655/148‑1 (MQ=255)
                       cagttcaTCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATAAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCg                                                                                                                           >  2:317361/1‑148 (MQ=255)
                         gttcaTCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACATTACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGc                                                                                                                          >  5:310023/1‑147 (MQ=255)
                         gttcaTCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGc                                                                                                                          >  1:112107/1‑147 (MQ=255)
                          ttcaTCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGTCAGCCCAGAGTGGTGGCAAATCCGcc                                                                                                                         >  4:26133/1‑147 (MQ=255)
                          ttcaTCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGcc                                                                                                                         >  5:236166/1‑147 (MQ=255)
                          ttcaTCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAACAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGTCAGCCCAGAGTGGTGGCAAATCCGcc                                                                                                                         >  2:350876/1‑147 (MQ=255)
                          ttcaTCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAACAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGTCAGCCCAGAGTGGTGGCAAATCCGcc                                                                                                                         >  2:345631/1‑147 (MQ=255)
                               cGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGca                                                                                                                    >  5:273889/1‑147 (MQ=255)
                                 aCATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACaa                                                                                                                 <  6:273889/148‑1 (MQ=255)
                                   aTGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCT‑CCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGt                                                                                                              <  4:92147/148‑1 (MQ=255)
                                     taaTAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATAAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGCGTGGTGGCAAATCCGCCAAGCACAAAGtt                                                                                                             >  6:179188/2‑148 (MQ=255)
                                     gAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGtt                                                                                                             >  3:127640/1‑148 (MQ=255)
                                       aTAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATCCTGACCAGTGTCGCCAACCGCCAGCCCAGCGTGGTGGCAAATCCGCCAAGCACAAAGtt                                                                                                             >  3:291200/1‑146 (MQ=255)
                                       aTAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGTGTCGCCAACAGCCAGCCGAGAGTGGTGGAAAATCCGCCAAGCACAAAGtt                                                                                                             >  3:211167/1‑146 (MQ=255)
                                       aTAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCTCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCa                                                                                                           <  6:85888/148‑1 (MQ=255)
                                       aTAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCa                                                                                                           <  6:90997/148‑1 (MQ=255)
                                         aGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTTGTGGCAAATCCGCCAAGCACAAAGTTCAGa                                                                                                         <  4:108252/148‑1 (MQ=255)
                                            ttCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATg                                                                                                      <  2:374499/148‑1 (MQ=255)
                                             tCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGt                                                                                                     >  7:244760/1‑148 (MQ=255)
                                                      tAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACgg                                                                                            <  1:317361/148‑1 (MQ=255)
                                                       aaGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACgg                                                                                            >  6:211482/1‑147 (MQ=255)
                                                           cACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACg                                                                                        >  2:348815/1‑147 (MQ=255)
                                                            aCCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCa                                                                                      >  7:213332/1‑148 (MQ=255)
                                                              cAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCAt                                                                                     >  6:250671/1‑147 (MQ=255)
                                                                gagaCAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGTCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATaaa                                                                                  <  3:26133/148‑1 (MQ=255)
                                                                gagaCAGTTTAGTGATCTCCCAACAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGTCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATaaa                                                                                  <  1:345631/148‑1 (MQ=255)
                                                                gagaCAGTTTAGTGATCTCCCAACAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGTCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATaaa                                                                                  <  1:350876/148‑1 (MQ=255)
                                                                 agaCAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATaaaa                                                                                 <  4:200793/148‑1 (MQ=255)
                                                                  gaCAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAAc                                                                                >  7:355756/1‑148 (MQ=255)
                                                                  gaCAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAAc                                                                                <  6:111578/148‑1 (MQ=255)
                                                                   aCAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGCCCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCCAGCACAAAGTTCAGAAGGTTCAACACGGTACGCATTACAcc                                                                               >  2:281083/1‑148 (MQ=255)
                                                                               tctcCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACtgt                                                                   >  8:367599/1‑148 (MQ=255)
                                                                                    cagcagGATCGTGTCAGCTGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACATCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGa                                                                                                         <  4:195897/105‑1 (MQ=255)
                                                                                    cagcagGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGa                                                                                                         >  3:195897/1‑105 (MQ=255)
                                                                                    cagcagGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGttt                                                              >  8:65536/1‑148 (MQ=255)
                                                                                    cagcagGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTC‑AAACGGTACGCATAAAACCTCGCTTTACTGTGGtttt                                                             <  8:113386/148‑1 (MQ=255)
                                                                                      gcagGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTc                                                            <  4:169039/148‑1 (MQ=255)
                                                                                      gcagGATCGTGTCAGAGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCTCCAACAGCCATCCCAGAGTGGTGGCAAATCCGCCAAGGACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTc                                                            <  5:179188/148‑1 (MQ=255)
                                                                                      gcagGATCGTGGCAGCGGTAAGGTAAAAATCAGCACAATACTGCCCAGAGTCGCCAACAGCCAGCCCAGAGGGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTc                                                            <  4:243455/148‑1 (MQ=255)
                                                                                      gcagGAGCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTc                                                            <  5:211482/148‑1 (MQ=255)
                                                                                         ggATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTAGGCATAAAACCTCGCTTTACTGTGGTTTTCAGt                                                         <  8:41134/148‑1 (MQ=255)
                                                                                          gATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTa                                                        <  8:355756/148‑1 (MQ=255)
                                                                                          gATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTa                                                        <  8:273162/148‑1 (MQ=255)
                                                                                          gATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTa                                                        <  1:348815/148‑1 (MQ=255)
                                                                                           aTCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTaa                                                       >  8:70729/1‑148 (MQ=255)
                                                                                              gtgtCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATg                                                     <  2:112107/147‑1 (MQ=255)
                                                                                                gtCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGggc                                                  >  2:67387/1‑148 (MQ=255)
                                                                                                     cGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGa                                                                                                         <  6:142420/88‑1 (MQ=255)
                                                                                                     cGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGa                                                                                                         >  5:142420/1‑88 (MQ=255)
                                                                                                      ggtaaggtaaAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAAt                                            >  3:163959/1‑148 (MQ=255)
                                                                                                                   tCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGtttttt                               <  8:213332/148‑1 (MQ=255)
                                                                                                                   tCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGtttttt                               <  3:326045/148‑1 (MQ=255)
                                                                                                                         caATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCt                         >  5:261731/1‑148 (MQ=255)
                                                                                                                           aTACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAAt                                                      >  7:281837/1‑117 (MQ=255)
                                                                                                                           aTACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAAt                                                      <  8:281837/117‑1 (MQ=255)
                                                                                                                           aTACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTgg                       >  7:190369/1‑148 (MQ=255)
                                                                                                                                    agagTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGttt              <  7:65536/148‑1 (MQ=255)
                                                                                                                                     gagTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGtttt             >  7:283788/1‑148 (MQ=255)
                                                                                                                                      agTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTtc            >  1:190363/1‑148 (MQ=255)
                                                                                                                                      agTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTtc            <  4:163959/148‑1 (MQ=255)
                                                                                                                                         cGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTtc            <  8:20415/145‑1 (MQ=255)
                                                                                                                                         cGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTtc            >  7:20415/1‑145 (MQ=255)
                                                                                                                                             aaCAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTTCTCTGACg     >  2:149380/1‑148 (MQ=255)
                                                                                                                                                agccagccCAGAGTGGTGGCTAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTTCTCTGACGGtt  <  5:250671/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTCAGTAATGGGCGGCAATTGTAACGTTTTTTTGGGCTGGAGCACGTTTTCTCTGACGGTT  >  NC_000913/1024140‑1024431

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