Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,602,905 T→C 6.0% I363T (ATT→ACT)  hyfB → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,602,9050TC 6.0% 203.3 / 2.9 67I363T (ATT→ACT) hyfBhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (35/28);  new base C (2/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

TTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGG  >  NC_000913/2602758‑2603050
                                                                                                                                                   |                                                                                                                                                 
ttattTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGattatt                                                                                                                                                  >  1:83330/1‑149 (MQ=255)
     ttGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCCAAGGGCTGCTATGTCTCGGCGCGGGAGCGATTATTTcgcg                                                                                                                                             >  4:87266/1‑149 (MQ=255)
           gCAGTCGGTGTGGCGTTGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGc                                                                                                                                        <  7:316128/148‑1 (MQ=255)
                  gtggGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCGCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCATAGGGCTGCTATTGCTCGGCGCGGTAGCAATTATTTCGCGTTTGCATACCCAc                                                                                                                                <  6:387356/149‑1 (MQ=255)
                       gCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACAt                                                                                                                           >  6:110117/1‑149 (MQ=255)
                          aTGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTCCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATgg                                                                                                                         >  2:177518/1‑148 (MQ=255)
                            ggtcggtcTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGaa                                                                                                                       >  1:387482/1‑148 (MQ=255)
                              tcggtcTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTCTTCAAAGGGCTGCTATTTCTCGGCGCGGGACCGATTATTTCGCGTATGAATACCCACGACATGGaaaa                                                                                                                     >  1:361224/1‑148 (MQ=255)
                                 gtcTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATg                                                                                                                 <  4:72005/149‑1 (MQ=255)
                                 gtcTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATg                                                                                                                 <  3:292791/149‑1 (MQ=255)
                                 gtcTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATg                                                                                                                 >  7:72317/1‑149 (MQ=255)
                                  tcTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTCCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATgg                                                                                                                <  1:177518/149‑1 (MQ=255)
                                  gcTGTCACTGCACGACCCGCTGCTCAGCGTGGGTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGTCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAATAAATgg                                                                                                                <  2:109811/148‑1 (MQ=255)
                                    tGTCACTGCACGACCCGCTGCTCACCGGGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATgggg                                                                                                              >  6:127170/1‑149 (MQ=255)
                                      tCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATggg                                                                                                               <  5:278996/146‑1 (MQ=255)
                                           gCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATggg                                                                                                               <  6:64284/141‑1 (MQ=255)
                                           gCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATggg                                                                                                               >  5:64284/1‑141 (MQ=255)
                                            cACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGc                                                                                                      <  2:387482/149‑1 (MQ=255)
                                                   cGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACg                                                                                                >  4:147190/1‑148 (MQ=255)
                                                     ctgctCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGAt                                                                                             <  5:110117/149‑1 (MQ=255)
                                                           aCCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTg                                                                                       <  2:230374/149‑1 (MQ=255)
                                                                ggTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAg                                                                                  >  1:300240/1‑149 (MQ=255)
                                                                  ttGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGcc                                                                                <  1:151218/149‑1 (MQ=255)
                                                                   tGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCg                                                                               >  6:385039/1‑149 (MQ=255)
                                                                    ggACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTAGTTCTCGGCGCGGGAGCGCTTATTTCGCGTTTGCATACCCACGACTTGGCAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCg                                                                               >  3:284528/1‑148 (MQ=255)
                                                                       cTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCgcag                                                                            >  3:126826/1‑148 (MQ=255)
                                                                          cTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCgcagca                                                                          <  5:127170/147‑1 (MQ=255)
                                                                          cTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATg                                                                        >  5:109610/1‑149 (MQ=255)
                                                                          cTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATg                                                                        >  8:53171/1‑149 (MQ=255)
                                                                             ggCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCt                                                                     <  8:72545/149‑1 (MQ=255)
                                                                                 cACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGAtt                                                                 <  5:385039/149‑1 (MQ=255)
                                                                                   cTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTgg                                                               >  6:266170/1‑149 (MQ=255)
                                                                                    tGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGt                                                              >  8:197037/1‑149 (MQ=255)
                                                                                    tGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGt                                                              >  3:150524/1‑149 (MQ=255)
                                                                                      tttCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGCCGCGGGAGCGATTACTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTg                                                            >  8:11584/1‑149 (MQ=255)
                                                                                      tttCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGCCGCGGGAGCGATTACTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTg                                                            >  8:11577/1‑149 (MQ=255)
                                                                                              gCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGa                                                    >  8:44198/1‑149 (MQ=255)
                                                                                              gCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTACCTCGCGa                                                    >  7:388307/1‑149 (MQ=255)
                                                                                                     cATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGTGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGc                                             <  2:83330/149‑1 (MQ=255)
                                                                                                     cATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGc                                             >  1:119602/1‑149 (MQ=255)
                                                                                                       tGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCa                                           <  5:136740/149‑1 (MQ=255)
                                                                                                          gcTATTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTAGCCTCGCGATATCAGCCATtc                                        >  2:49319/1‑149 (MQ=255)
                                                                                                                cataGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATTGAAAAATTGGGGGGACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCg                                   <  2:172972/145‑1 (MQ=255)
                                                                                                                cAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAAc                                                                                                 >  8:370528/1‑86 (MQ=255)
                                                                                                                cAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAAc                                                                                                 <  7:370528/86‑1 (MQ=255)
                                                                                                                 aaaGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATtc                                        >  2:34544/1‑142 (MQ=255)
                                                                                                                 aaaGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATtc                                        <  1:34544/142‑1 (MQ=255)
                                                                                                                 aaaGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCt                                 >  4:375854/1‑149 (MQ=255)
                                                                                                                   aGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTg                                >  5:100888/1‑148 (MQ=255)
                                                                                                                     ggCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAAt                             >  8:188371/1‑149 (MQ=255)
                                                                                                                       ctgctATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTcgc                                                      <  1:245454/122‑1 (MQ=255)
                                                                                                                       ctgctATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTcgc                                                      >  2:245454/1‑122 (MQ=255)
                                                                                                                         gctATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCTCGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGtt                         <  4:342120/149‑1 (MQ=255)
                                                                                                                          ctATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGttt                        >  4:191708/1‑149 (MQ=255)
                                                                                                                          ctATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCAGCAGCATGAATGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGttt                        >  7:369405/1‑149 (MQ=255)
                                                                                                                            aTTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTa                      >  3:291234/1‑149 (MQ=255)
                                                                                                                             tttCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGtttt                       <  6:100888/147‑1 (MQ=255)
                                                                                                                               tctcGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATc                    <  4:291234/148‑1 (MQ=255)
                                                                                                                                ctcGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTATCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAg                  >  7:170408/1‑149 (MQ=255)
                                                                                                                                 tcGGCGCGGGAGCGATTATTTCGCGTTTGCAAACCCACGACATTGAAAAAATCCGGCCACTAGCGAAACTGATGCCGTGGACAGCCTCAGCATGCCTGATACGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGc                 >  8:296075/1‑149 (MQ=255)
                                                                                                                                    gcgcTGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGaa              <  3:87266/149‑1 (MQ=255)
                                                                                                                                      gcgGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAATAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATg            >  3:40102/1‑149 (MQ=255)
                                                                                                                                       cgGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATgg           >  8:28087/1‑149 (MQ=255)
                                                                                                                                            gCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTacac      <  4:126826/149‑1 (MQ=255)
                                                                                                                                             cGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACAcc     <  5:266170/149‑1 (MQ=255)
                                                                                                                                             cGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACAcc     <  8:59880/149‑1 (MQ=255)
                                                                                                                                             cGATTACTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACAcc     <  7:11584/149‑1 (MQ=255)
                                                                                                                                             cGATTACTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACAcc     <  7:11577/149‑1 (MQ=255)
                                                                                                                                                ttattTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTgg  <  7:44198/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTATTTTGCTGGCAGTCGGTGTGGCGATGGTCGGTCTGTCACTGCACGACCCGCTGCTCACCGTGGTTGGACTGCTCGGCGCACTGTTTCATCTGCTCAACCATGCGCTGTTCAAAGGGCTGCTATTTCTCGGCGCGGGAGCGATTATTTCGCGTTTGCATACCCACGACATGGAAAAAATGGGGGCACTAGCGAAACGGATGCCGTGGACAGCCGCAGCATGCCTGATTGGTTGCCTCGCGATATCAGCCATTCCTCCGCTGAATGGTTTTATCAGCGAATGGTACACCTGG  >  NC_000913/2602758‑2603050

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