Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,707,226 T→A 10.3% L103* (TTA→TAA) ‡ rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,707,2260TA10.3% 181.5 / 11.7 58L103* (TTA→TAA) ‡rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base T (30/22);  new base A (2/4);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-01

AAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGCTGGGCGGTGCAGGATTC  >  NC_000913/1707079‑1707374
                                                                                                                                                   |                                                                                                                                                    
aaGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTa                                                                                                                                                     <  4:99510/149‑1 (MQ=255)
      aGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTAAAGCTGaa                                                                                                                                               <  2:139451/149‑1 (MQ=255)
         tgtgCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTAAAGCTGaagga                                                                                                                                            <  2:309783/149‑4 (MQ=255)
         tgtgCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTAAAGCTGaagga                                                                                                                                            <  3:379361/149‑4 (MQ=255)
            gCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGc                                                                                                                                         <  2:290742/149‑1 (MQ=255)
                        gNGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGAt                                                                                                                             >  3:259872/1‑149 (MQ=255)
                          gATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGc                                                                                                                           >  3:384505/1‑149 (MQ=255)
                           aTAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGcc                                                                                                                          <  1:521584/149‑1 (MQ=255)
                                 tATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGt                                                                                                                    >  2:519104/1‑149 (MQ=255)
                                       gcgGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGAc                                                                                                              >  2:300284/1‑149 (MQ=255)
                                         ggCCAGTCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGActg                                                                                                            >  4:639662/1‑149 (MQ=255)
                                         ggCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGat                                                                                                                                    >  4:203747/1‑125 (MQ=255)
                                         ggCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGat                                                                                                                                    <  3:203747/125‑1 (MQ=255)
                                             aGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTgg                                                                                                        >  3:282911/1‑149 (MQ=255)
                                               ccGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGc                                                                                                                                                   <  3:193672/104‑1 (MQ=255)
                                               ccGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGc                                                                                                                                                   >  4:193672/1‑104 (MQ=255)
                                                               gcgGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTgg                                                                                      >  4:651479/1‑149 (MQ=255)
                                                                 ggCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGAGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGc                                                                                    >  4:585725/1‑149 (MQ=255)
                                                                   cAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCg                                                                                  >  4:265661/1‑149 (MQ=255)
                                                                   cAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCg                                                                                  <  1:519104/149‑1 (MQ=255)
                                                                    aaaaTGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGa                                                                                 >  4:500755/1‑149 (MQ=255)
                                                                    aaaaTGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGa                                                                                 >  3:268792/1‑149 (MQ=255)
                                                                    aaaaTGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGa                                                                                 >  4:248408/1‑149 (MQ=255)
                                                                         gctgcCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATc                                                                            <  2:271500/149‑1 (MQ=255)
                                                                         gctgcCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATc                                                                            <  4:268792/149‑1 (MQ=255)
                                                                           tgcCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTAGTGATGCCGATGGTGGAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGc                                                                          >  4:40990/1‑149 (MQ=255)
                                                                                 ttCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGc                                                                    >  2:190310/1‑149 (MQ=255)
                                                                                  tCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCa                                                                   <  3:211569/149‑1 (MQ=255)
                                                                                    aCGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGt                                                                 >  2:556892/1‑149 (MQ=255)
                                                                                      gcgcCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTcg                                                               >  2:633402/1‑149 (MQ=255)
                                                                                          ccACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGaa                                                           <  4:259872/149‑1 (MQ=255)
                                                                                          ccACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGaa                                                           <  4:77298/149‑1 (MQ=255)
                                                                                           cACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAag                                                          >  1:504813/1‑149 (MQ=255)
                                                                                            aCCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAaga                                                         <  1:347374/149‑1 (MQ=255)
                                                                                            aCCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAaga                                                         <  3:639662/149‑1 (MQ=255)
                                                                                                cGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTa                                                     <  3:585725/149‑1 (MQ=255)
                                                                                                          tACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTGATCGAGCGCa                                           <  3:651479/149‑1 (MQ=255)
                                                                                                                 aTTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACTTCa                                    <  4:588365/149‑1 (MQ=255)
                                                                                                                       ccccACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCa                                           >  4:376578/1‑136 (MQ=255)
                                                                                                                       ccccACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCa                                           <  3:376578/136‑1 (MQ=255)
                                                                                                                       ccccACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTgg                              <  3:665949/149‑1 (MQ=255)
                                                                                                                        cccACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAagag                                                        >  3:23180/1‑122 (MQ=255)
                                                                                                                        cccACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAagag                                                        <  4:23180/122‑1 (MQ=255)
                                                                                                                        ccaaCTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGAGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACTGCTGGGCCGATTATCGCACTCGCCGTCGCGAAGAGCTAATCGAGCGCATACATCAGTTTGgt                             >  1:69649/4‑149 (MQ=255)
                                                                                                                           aCTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGtt                          >  4:195780/1‑149 (MQ=255)
                                                                                                                             tctACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGc                        <  3:248408/149‑1 (MQ=255)
                                                                                                                             tctACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGc                        >  2:180493/1‑149 (MQ=255)
                                                                                                                                aCGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCgg                     >  4:217303/1‑149 (MQ=255)
                                                                                                                                aCGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCgg                     >  2:125637/1‑149 (MQ=255)
                                                                                                                                  ggCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGc                   >  3:508801/1‑149 (MQ=255)
                                                                                                                                     tCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGCTgg                >  3:506918/1‑149 (MQ=255)
                                                                                                                                     tCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGCTgg                >  4:411739/1‑149 (MQ=255)
                                                                                                                                       ttaattCAGCTAAAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAaga                                                         >  1:192312/5‑106 (MQ=255)
                                                                                                                                       ttaattCAGCTAAAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAaga                                                         <  2:192312/102‑1 (MQ=255)
                                                                                                                                       ttaattCAGCTAAAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGCTGGGc              >  3:156445/5‑149 (MQ=255)
                                                                                                                                           ttCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGAGCTGGGCGGTg          >  4:584314/1‑149 (MQ=255)
                                                                                                                                             cAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGCTGGGCGGTGCa        <  1:248039/149‑1 (MQ=255)
                                                                                                                                             cAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGAGCTGGGCGGTGCa        <  3:584314/149‑1 (MQ=255)
                                                                                                                                                 tttAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGtt                          <  1:575744/127‑1 (MQ=255)
                                                                                                                                                 tttAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGtt                          >  2:575744/1‑127 (MQ=255)
                                                                                                                                                 tttAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGCTGGGCGGTGCAGGAt    <  1:190310/149‑1 (MQ=255)
                                                                                                                                                   tAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGCTGGGCGGTGCAGGATTc  <  4:506918/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGCTGGGCGGTGCAGGATTC  >  NC_000913/1707079‑1707374

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