Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,082,066 C→T 5.5% D232N (GAC→AAC)  fdoH ← formate dehydrogenase‑O, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,082,0660CT 5.5% 234.4 / 2.3 73D232N (GAC→AAC) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (38/30);  new base T (2/2);  total (41/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 = 5.89e-01

GGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGC  >  NC_000913/4081918‑4082213
                                                                                                                                                    |                                                                                                                                                   
ggACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTc                                                                                                                                                     <  3:116718/149‑1 (MQ=255)
         cccACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTg                                                                                                                                            <  6:225307/149‑1 (MQ=255)
          ccACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTg                                                                                                                                            >  5:281220/1‑148 (MQ=255)
           cACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTgca                                                                                                                                          <  1:274177/149‑1 (MQ=255)
           cACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTgca                                                                                                                                          <  6:159178/149‑1 (MQ=255)
            aCGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTgcag                                                                                                                                         >  4:5904/1‑149 (MQ=255)
                    gAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTAc                                                                                                                                  <  8:340733/148‑1 (MQ=255)
                    gAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTAc                                                                                                                                  <  8:299301/148‑1 (MQ=255)
                    gAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACa                                                                                                                                 >  4:9329/1‑149 (MQ=255)
                    gAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCATTTCCAGATTTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACa                                                                                                                                 >  8:297791/1‑149 (MQ=255)
                       gaaagtGGCCGGAAAGGTAGCCGCATAGCCAACAGCCGCCAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGa                                                                                                                              >  6:237090/6‑149 (MQ=255)
                          aCTGCCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATATCACTGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACgt                                                                                                                           <  4:80275/149‑1 (MQ=255)
                           cTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACgtg                                                                                                                          <  6:281220/149‑1 (MQ=255)
                            tGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCATGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACgtgt                                                                                                                         >  7:279552/1‑149 (MQ=255)
                            tGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTgca                                                                                                                                          <  3:255474/132‑1 (MQ=255)
                            tGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTgca                                                                                                                                          >  4:255474/1‑132 (MQ=255)
                              gCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACgtgt                                                                                                                         <  8:174398/147‑1 (MQ=255)
                              gCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATATCACGGTATCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTCGCTTGTCAGCATGGTGCAGCACGTACATGACgtgtg                                                                                                                        >  1:78366/1‑148 (MQ=255)
                                cGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTac                                                                                                                     >  8:100593/1‑149 (MQ=255)
                                 gCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTac                                                                                                                     >  8:60509/1‑148 (MQ=255)
                                   aaaGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTTCAGCACGTACATGACGTCTGTaccac                                                                                                                  >  3:44832/1‑149 (MQ=255)
                                   aaaGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTaccac                                                                                                                  >  8:2972/1‑149 (MQ=255)
                                   aaaGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCAGGCAAGCCATGATACAGATTTGGCTTGTAAGCATGGTGCAGCACGTACATGACGAGAGTaccac                                                                                                                  >  7:124560/1‑149 (MQ=255)
                                       gTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGAc                                                                                                              >  2:263964/1‑149 (MQ=255)
                                             gCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGc                                                                                                        >  8:250794/1‑149 (MQ=255)
                                               aaaGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCg                                                                                                      >  1:289726/1‑149 (MQ=255)
                                                   ccAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATc                                                                                                  <  5:299199/149‑1 (MQ=255)
                                                     aaCAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCg                                                                                                 >  5:177512/1‑148 (MQ=255)
                                                      aCAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTa                                                                                               <  3:5904/149‑1 (MQ=255)
                                                        aGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACa                                                                                             <  4:225392/149‑1 (MQ=255)
                                                                        cAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGt                                                                              >  3:176183/1‑148 (MQ=255)
                                                                         aGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTTAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTaa                                                                            >  7:241763/1‑149 (MQ=255)
                                                                         aGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTTAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTaa                                                                            >  7:241764/1‑149 (MQ=255)
                                                                          gATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAAc                                                                           >  7:134452/1‑149 (MQ=255)
                                                                             gCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACcgc                                                                        >  6:182780/1‑149 (MQ=255)
                                                                               cTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCgg                                                                      <  1:179734/149‑1 (MQ=255)
                                                                               cTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGCTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCgg                                                                      <  7:323302/149‑1 (MQ=255)
                                                                                tttCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCggg                                                                     >  8:174953/1‑149 (MQ=255)
                                                                                 ttCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGt                                                                    >  2:156894/1‑149 (MQ=255)
                                                                                    cAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGtttt                                                                 >  4:253614/1‑149 (MQ=255)
                                                                                    cAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGtttt                                                                 >  4:47235/1‑149 (MQ=255)
                                                                                     aGAATTTCTCGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGGACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTtc                                                                <  2:78366/149‑1 (MQ=255)
                                                                                         tttCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAACATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTtcagct                                                            >  2:172647/1‑149 (MQ=255)
                                                                                           tCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTtcagct                                                            >  1:46345/1‑147 (MQ=255)
                                                                                             aCGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTtcagctcag                                                         >  1:36280/1‑148 (MQ=255)
                                                                                                     gCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGcgct                                                >  6:61860/1‑149 (MQ=255)
                                                                                                         aTTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGcgctcgct                                            >  6:297045/1‑149 (MQ=255)
                                                                                                           ttCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGCCCCGAATTGTCGTAACCGCGGGTGTTCAGCTCAGCAACCCGctcgctcg                                          >  4:234952/1‑149 (MQ=255)
                                                                                                            tCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGGGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGctcgctcg                                          <  7:100593/148‑1 (MQ=255)
                                                                                                              cGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCtcgctcgc                                         >  3:22346/1‑147 (MQ=255)
                                                                                                               gggTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAg                                      >  8:751/1‑149 (MQ=255)
                                                                                                                 gTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCg                                    <  6:177512/149‑1 (MQ=255)
                                                                                                                     tcCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGtttt                                >  8:338405/1‑149 (MQ=255)
                                                                                                                       cGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCa                              >  1:214570/1‑149 (MQ=255)
                                                                                                                       cGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCa                              <  2:214570/149‑1 (MQ=255)
                                                                                                                        ggCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCAt                             <  2:289726/149‑1 (MQ=255)
                                                                                                                        ggCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCAt                             <  5:182780/149‑1 (MQ=255)
                                                                                                                           aaGCCATGATACAGATTTGGCTTGTTAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGa                          <  8:241763/149‑1 (MQ=255)
                                                                                                                           aaGCCATGATACAGATTTGGCTTGGTAGCATGGTGCAGAACGGACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCGCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGa                          <  8:241764/149‑1 (MQ=255)
                                                                                                                             gCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGAc                         <  3:320947/148‑1 (MQ=255)
                                                                                                                             gCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGAc                         >  4:320947/1‑148 (MQ=255)
                                                                                                                              ccATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTc                                                                                <  4:138220/92‑1 (MQ=255)
                                                                                                                              ccATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTc                                                                                >  3:138220/1‑92 (MQ=255)
                                                                                                                              ccATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGActc                       >  2:141778/1‑149 (MQ=255)
                                                                                                                                aTGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCtt                     <  5:242852/149‑1 (MQ=255)
                                                                                                                                     aCAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTa                >  4:150058/1‑149 (MQ=255)
                                                                                                                                      cAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTa                <  7:2972/148‑1 (MQ=255)
                                                                                                                                      cAGATTTGGCTTGTCAGCAGGGTGCAGCACGTACATGACGTGTGTCCCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCGCGCTCGCCAGCGTTTTCATCGACTCTTTCGTa                <  6:25238/148‑1 (MQ=255)
                                                                                                                                             ggCTTGTCAGCATGGTGCAGCACGTACATGACGTGGGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGt         <  7:297791/148‑1 (MQ=255)
                                                                                                                                              gCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTg        <  3:152366/148‑1 (MQ=255)
                                                                                                                                              gCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGa       <  5:61860/149‑1 (MQ=255)
                                                                                                                                                 tGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTAAAGACCCGCATTGTCGTAACCGCGGGTTTTCAGATCAGCACCGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATc    >  8:280496/1‑149 (MQ=255)
                                                                                                                                                   tCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATcgc  <  7:338405/149‑1 (MQ=255)
                                                                                                                                                   tCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATcgc  <  7:174953/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTCGCTGATTTCCGGGTTCTCCGGCAAGCCATGATACAGATTTGGCTTGTCAGCATGGTGCAGCACGTACATGACGTGTGTACCACCGACGCCTGCCGGATCGTACAGACCCGCATTGTCGTAACCGCGGGTTTTCAGCTCAGCAACGCGCTCGCTCGCCAGCGTTTTCATCGACTCTTTCGTACCGAAGTGAATCGC  >  NC_000913/4081918‑4082213

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