Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,707,225 T→A 10.3% L103I (TTA→ATA) ‡ 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,2250TA10.3% 178.0 / 11.3 58L103I (TTA→ATA) ‡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 = 6.54e-01

AAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGCTGGGCGGTGCAGGAT  >  NC_000913/1707079‑1707372
                                                                                                                                                  |                                                                                                                                                   
aaGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTa                                                                                                                                                   <  4:99510/149‑1 (MQ=255)
      aGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTAAAGCTGaa                                                                                                                                             <  2:139451/149‑1 (MQ=255)
         tgtgCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTAAAGCTGaagga                                                                                                                                          <  3:379361/149‑4 (MQ=255)
         tgtgCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTAAAGCTGaagga                                                                                                                                          <  2:309783/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                                                                                                                                                 >  4:193672/1‑104 (MQ=255)
                                               ccGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGc                                                                                                                                                 <  3:193672/104‑1 (MQ=255)
                                                               gcgGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTgg                                                                                    >  4:651479/1‑149 (MQ=255)
                                                                 ggCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGAGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGc                                                                                  >  4:585725/1‑149 (MQ=255)
                                                                   cAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCg                                                                                <  1:519104/149‑1 (MQ=255)
                                                                   cAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCg                                                                                >  4:265661/1‑149 (MQ=255)
                                                                    aaaaTGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGa                                                                               >  4:248408/1‑149 (MQ=255)
                                                                    aaaaTGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGa                                                                               >  3:268792/1‑149 (MQ=255)
                                                                    aaaaTGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGa                                                                               >  4:500755/1‑149 (MQ=255)
                                                                         gctgcCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATc                                                                          <  4:268792/149‑1 (MQ=255)
                                                                         gctgcCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATc                                                                          <  2:271500/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:77298/149‑1 (MQ=255)
                                                                                          ccACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGaa                                                         <  4:259872/149‑1 (MQ=255)
                                                                                           cACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAag                                                        >  1:504813/1‑149 (MQ=255)
                                                                                            aCCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAaga                                                       <  3:639662/149‑1 (MQ=255)
                                                                                            aCCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAaga                                                       <  1:347374/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              >  4:411739/1‑149 (MQ=255)
                                                                                                                                     tCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGCTgg              >  3:506918/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)
                                                                                                                                                  |                                                                                                                                                   
AAGGTGAGTTGTGCGTTAGCGTCGGCGATAAAGTATTGCGCGGCCAGCCGCTTACCCGTGGTCGCGGCAAAATGCTGCCTGTTCACGCGCCCACCTCGGGTACCGTTACGGCTATTGCGCCCCACTCTACGGCTCATCCTTCAGCTTTAGCTGAATTAAGCGTGATTATTGATGCCGATGGTGAAGACTGCTGGATCCCGCGCGACGGCTGGGCCGATTATCGCACTCGCAGTCGCGAAGAGTTAATCGAGCGCATACATCAGTTTGGTGTTGCCGGGCTGGGCGGTGCAGGAT  >  NC_000913/1707079‑1707372

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: