Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,024,226 G→A 5.2% C41C (TGC→TGT)  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,2260GA 5.2% 248.8 / 2.7 76C41C (TGC→TGT) yccFDUF307 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (40/32);  new base A (2/2);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

AACCACCAGCCAAAGAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTC  >  NC_000913/1024080‑1024373
                                                                                                                                                  |                                                                                                                                                   
aaCCACCAGCCAAAGAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCcag                                                                                                                                                     >  3:169039/1‑147 (MQ=255)
aaCCACCAGCCAAAGAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCcag                                                                                                                                                     >  3:108252/1‑147 (MQ=255)
     ccaGCCAAAGAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGAt                                                                                                                                               <  8:317689/148‑1 (MQ=255)
     ccaGCCAAAGAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGAt                                                                                                                                               >  4:326045/1‑148 (MQ=255)
     ccaGCCAAAGAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGAt                                                                                                                                               <  1:38584/148‑1 (MQ=255)
          cAAAGAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCgtgt                                                                                                                                          <  1:313632/148‑1 (MQ=255)
           aaaGAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCgtgt                                                                                                                                          >  7:8655/1‑147 (MQ=255)
           aaaGAAAATCAGCCAGAAAATATTCAAGACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCgtgt                                                                                                                                          >  7:10428/1‑147 (MQ=255)
             aGAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAg                                                                                                                                       >  4:256394/1‑148 (MQ=255)
              gAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGc                                                                                                                                      <  7:161641/148‑1 (MQ=255)
              gAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGc                                                                                                                                      >  8:311368/1‑148 (MQ=255)
                     aGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCggtaag                                                                                                                                <  8:310680/147‑1 (MQ=255)
                       ccAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCT‑CCAGCAGGATCGTGTCAGCggtaaggta                                                                                                                             >  3:92147/1‑147 (MQ=255)
                        cAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCGGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCggtaaggta                                                                                                                             >  7:222351/1‑147 (MQ=255)
                         aGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCggtaaggtaa                                                                                                                            >  7:28704/1‑147 (MQ=255)
                          gAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTaaa                                                                                                                           >  6:73685/1‑147 (MQ=255)
                           aaaaTATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTaaaaa                                                                                                                         >  5:111578/1‑148 (MQ=255)
                            aaaTATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGCTCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAAt                                                                                                                        >  5:85888/1‑148 (MQ=255)
                            aaaTATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAAt                                                                                                                        >  5:90997/1‑148 (MQ=255)
                             aaTATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATc                                                                                                                       <  3:181912/148‑1 (MQ=255)
                                 ttCAATCCCCTCCCGCCAGTAGTCAGCAGCACCTTTTTGCCAGACGGGTTCAGTTCATCGACAAGAATCGCTTCATTGCCATACGGCACCAGACACAGTTTAGTGCTCTCCCAGCAGGCTCGTGTCAGCGGTAAGGTAAAAATCAGc                                                                                                                    <  2:326763/147‑1 (MQ=255)
                                  tCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGca                                                                                                                   >  7:273162/1‑147 (MQ=255)
                                    aaTACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACaa                                                                                                                <  8:14253/148‑1 (MQ=255)
                                     aTACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAAt                                                                                                               <  1:79737/148‑1 (MQ=255)
                                       accgtaccgCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCg                                                                                                                                             >  5:34208/1‑116 (MQ=255)
                                       accgtaccgCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCg                                                                                                                                             <  6:34208/116‑1 (MQ=255)
                                        ccgtaccgCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATAc                                                                                                             <  8:125714/147‑1 (MQ=255)
                                         cgtaccgCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTg                                                                                                           <  1:315304/148‑1 (MQ=255)
                                              cgCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCag                                                                                                      <  7:311368/148‑1 (MQ=255)
                                              cgCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCag                                                                                                      <  8:186650/148‑1 (MQ=255)
                                                    tATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGc                                                                                                <  8:10428/148‑1 (MQ=255)
                                                     aTTCATCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGc                                                                                                <  3:256394/147‑1 (MQ=255)
                                                     aTTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGc                                                                                                >  7:233769/1‑147 (MQ=255)
                                                     aTTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGc                                                                                                >  1:374499/1‑147 (MQ=255)
                                                         agcagcACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAAc                                                                                            <  8:28704/147‑1 (MQ=255)
                                                          gcagcaCATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAAcag                                                                                          <  4:326981/148‑1 (MQ=255)
                                                            agcaCATTTTTGCCGGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCATCACAATACTGACCAGAGTCGCCAAcagcc                                                                                        <  8:222351/148‑1 (MQ=255)
                                                               acaTTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAAcagccagc                                                                                     >  7:41134/1‑148 (MQ=255)
                                                               acaTTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAAcagccagc                                                                                     >  3:200793/1‑148 (MQ=255)
                                                                   tttaGCCAGCCTGGTTCAGTTCATCTACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCTTGTCAGCGGTAAGGTTAAAATCa                                                                                                                      >  1:314936/1‑111 (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:345631/1‑147 (MQ=255)
                                                                                      ttcaTCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAACAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGTCAGCCCAGAGTGGTGGCAAATCCGcc                                                               >  2:350876/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)
                                                                                                tGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCaga                                                                                                     >  3:246054/1‑99 (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)
                                                                                                            ttGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAAt                                                                                                               <  5:262291/77‑1 (MQ=255)
                                                                                                            ttGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAAt                                                                                                               >  6:262291/1‑77 (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:350876/148‑1 (MQ=255)
                                                                                                                            gagaCAGTTTAGTGATCTCCCAACAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGTCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATaaa                        <  1:345631/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)
                                                                                                                                                  |                                                                                                                                                   
AACCACCAGCCAAAGAAAATCAGCCAGAAAATATTCAATACCGTACCGCCAGTATTCAGCAGCACATTTTTGCCAGCCGGGTTCAGTTCATCGACATGAATAGCTTCATTGCCATAAGGCACCAGAGACAGTTTAGTGATCTCCCAGCAGGATCGTGTCAGCGGTAAGGTAAAAATCAGCACAATACTGACCAGAGTCGCCAACAGCCAGCCCAGAGTGGTGGCAAATCCGCCAAGCACAAAGTTCAGAATGTTCAAAACGGTACGCATAAAACCTCGCTTTACTGTGGTTTTC  >  NC_000913/1024080‑1024373

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: ‑