Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,041,217 G→A 5.8% Q73* (CAG→TAG) ‡ mobB ← molybdopterin‑guanine dinucleotide biosynthesis adaptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,041,2170GA 5.8% 189.0 / 3.0 69Q73* (CAG→TAG) ‡mobBmolybdopterin‑guanine dinucleotide biosynthesis adaptor protein
Reads supporting (aligned to +/- strand):  ref base G (25/40);  new base A (2/2);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

CACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATC  >  NC_000913/4041077‑4041353
                                                                                                                                            |                                                                                                                                        
cacaATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGTGAAACTGTAGATGCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTgctggc                                                                                                                                        >  1:512177/1‑143 (MQ=255)
cacaATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTGGTCTGGCGCTGCCGTCCTCAAGGCCCCTCGTTgctggc                                                                                                                                        >  1:749230/1‑143 (MQ=255)
cacaATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCATGGCCCATCGTTgttggc                                                                                                                                        >  1:464069/1‑143 (MQ=255)
 acaATCTTTCCGATCTCTTCATGCTTAAACCCTTCGACCAGAAGCAAATCCAGCTTTGCGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTgctggct                                                                                                                                       <  1:285605/143‑1 (MQ=255)
   aaTCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGctggctgg                                                                                                                                     >  1:167151/1‑143 (MQ=255)
     tCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCGAGGCCCATCGTTCCtggctggct                                                                                                                                   >  1:394960/1‑142 (MQ=255)
     tCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCACGGCCCATCGTTGCTGGCTGGCa                                                                                                                                   >  2:266759/1‑143 (MQ=255)
         tGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCGTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGa                                                                                                                               <  2:440161/143‑1 (MQ=255)
            gATCTCTTCATGCTTAAGCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCACCAAGGCCCATCGTTGCCAGCCAGCAACGATgg                                                                                                                            <  2:579345/143‑1 (MQ=255)
            gATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATgg                                                                                                                            <  1:234354/143‑1 (MQ=255)
             aTCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTCTCCATTCGACTCGCGCGCAACTGTAGATCCCAATCTGGATGCTCTGGTGGTTCTGGCGTCCAGGCACATCGGTGGTGGCTGGCGACggttgt                                                                                                                           >  1:252620/1‑138 (MQ=255)
                tctTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTtg                                                                                                                        <  1:1000575/143‑1 (MQ=255)
                 ctTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGCCTGGTGTTTCCGCCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTtgt                                                                                                                       <  1:309657/143‑1 (MQ=255)
                     aTGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAAAGACGGACTGTGCcg                                                                                                                   >  2:246107/1‑143 (MQ=255)
                      ttcctAAACCCTTCGACCAGAATCAACTCCAGCTTTGAGGTATCCATTCGCCTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCcgc                                                                                                                  <  1:783727/139‑1 (MQ=255)
                      ggCTTAAGCCCTTCGACCAGAATCAAATCCAGCTGTGAGGTATCCAGTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCcgc                                                                                                                  <  1:219246/142‑1 (MQ=255)
                       gCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCCCCGCTGGTTTGCGCcgcg                                                                                                                 >  1:833207/1‑143 (MQ=255)
                         ttAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGcc                                                                                                               <  1:438119/143‑1 (MQ=255)
                             aCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTGGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGcc                                                                                                           >  1:9979/1‑143 (MQ=255)
                             aCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCGATTGTTGCTGCGCGCCAACGATGGTTTGTGCCGCGCCAGcc                                                                                                           >  2:11033/1‑143 (MQ=255)
                                cTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTg                                                                                                        >  1:154262/1‑143 (MQ=255)
                                  tCGACCCGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcg                                                                                                      <  1:1000325/143‑1 (MQ=255)
                                         gAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAACGCCCATCGTTGCTGCCTGGCCCCGGTCGTTTGTCCCGCGCCCGCCTCGCGCatccca                                                                                               >  2:851652/1‑138 (MQ=255)
                                          aaTCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCACAt                                                                                              >  2:412996/1‑143 (MQ=255)
                                            tCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCGAGCCTTGCGCAGCTCATAg                                                                                            >  1:552914/1‑143 (MQ=255)
                                              aaaTCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGCCCCATCGTTGCTGGCTGGCAACGCTGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCCTCGCt                                                                                          >  1:498037/1‑143 (MQ=255)
                                                  ccAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCt                                                                                      <  1:617202/143‑1 (MQ=255)
                                                      cgttGAGGTATGCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTGATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGc                                                                                  <  2:186112/141‑1 (MQ=255)
                                                      cTTTGGGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGc                                                                                  <  1:402346/143‑1 (MQ=255)
                                                         tGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTg                                                                               <  1:721877/143‑1 (MQ=255)
                                                          gatgTATCCATTCGACTTGCGAGCAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTgg                                                                              <  1:62973/140‑1 (MQ=255)
                                                          gAGGTATCCATTCGACTTGCGGGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTgg                                                                              <  1:14467/143‑1 (MQ=255)
                                                           aGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGc                                                                             <  2:98795/143‑1 (MQ=255)
                                                           aGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGc                                                                             <  1:829359/143‑1 (MQ=255)
                                                               acccACTCGCCTGGCGAGAAACTGTAGAGCCGGCTCTGCTTCGTCCGGTTTTTCCGTCATCAAGGCCCACCCTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTAt                                                                         <  2:438393/141‑1 (MQ=255)
                                                                tCCAGTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                        <  1:73974/143‑1 (MQ=255)
                                                                 ccATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCa                                                                       <  1:938584/143‑1 (MQ=255)
                                                                   aTTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc                                                                     <  1:25629/143‑1 (MQ=255)
                                                                   aTTCGACGTGCGAGTAGCTGTGGATCCAGCTCTGCTGCGGCTGGTGTTTCCGTCATCAAGGCCCACCGTGGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc                                                                     <  2:87137/143‑1 (MQ=255)
                                                                     tCGACTTGCGAGAAACTGTAGATCCAGCCCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACAt                                                                   <  1:529064/143‑1 (MQ=255)
                                                                      cGACTTGCGAGAAACTGTAGATCCAGCTCTTCGTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATc                                                                  <  2:290618/143‑1 (MQ=255)
                                                                       gACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCCTATCAACATcc                                                                 <  1:1075492/143‑1 (MQ=255)
                                                                        aCTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCa                                                                <  2:795368/143‑1 (MQ=255)
                                                                             cGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGAGGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCAGAGCTATCTTTGCCTGGCTTATCAACGTCCATATCa                                                           >  2:696599/1‑143 (MQ=255)
                                                                               tgaAACGGTAGGGCGAGCTCTTCGTGGTCGGGTGTTTCGGTCATCAAGGGCCATCGTGGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATg                                                         <  2:743835/142‑1 (MQ=255)
                                                                                gaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTGATAGCTATCTTTGCCTGGCTTATCAACATCCATATCCTgg                                                        >  1:687236/1‑143 (MQ=255)
                                                                                 aaaCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGCTGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCCTCTCACGGt                                                       >  2:1052879/1‑143 (MQ=255)
                                                                                      gTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTa                                                  <  1:485641/143‑1 (MQ=255)
                                                                                       tAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTAt                                                 >  1:704107/1‑143 (MQ=255)
                                                                                          gtCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCt                                              <  1:23112/142‑1 (MQ=255)
                                                                                          aTCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCAATATCATGGTGCGTATgat                                              >  1:1048806/1‑141 (MQ=255)
                                                                                            ccAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTa                                            >  1:878447/1‑143 (MQ=255)
                                                                                             cAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTaa                                           <  1:280590/143‑1 (MQ=255)
                                                                                                  cttcttCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGc                                      >  1:965014/1‑143 (MQ=255)
                                                                                                  cttcttCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGc                                      <  2:167151/143‑1 (MQ=255)
                                                                                                         ggcTGGTGGTCCCGCCATCAAGGCCCGTCGTTGCTGGCTGGCAACGATGGTTTGTCCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGAc                               <  2:997141/141‑1 (MQ=255)
                                                                                                         gTCTGTTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGAc                               <  1:108836/143‑1 (MQ=255)
                                                                                                               tttttCCGTAATCAAGGCCTATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATTTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATcc                         <  2:464069/141‑1 (MQ=255)
                                                                                                               tgtCTCCGCTATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATcc                         <  2:749230/143‑1 (MQ=255)
                                                                                                                 tttCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGAc                               <  1:949987/135‑1 (MQ=255)
                                                                                                                 tttCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGAc                               >  2:949987/1‑135 (MQ=255)
                                                                                                                   tCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATTGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCctct                     <  1:737869/143‑1 (MQ=255)
                                                                                                                     cGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTgg                   >  1:496104/1‑143 (MQ=255)
                                                                                                                         atcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCa               >  2:1046342/1‑143 (MQ=255)
                                                                                                                           cacaaaccATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGTTCATAGCTATCTTTGCCTGGCTTATCAACATCCATCTCATGGTGCGTATGCTTGATCCGCCCTGGACGGATCCCTCTGGGGCCTa             >  1:73517/7‑143 (MQ=255)
                                                                                                                           caAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCaca             >  1:1072906/1‑141 (MQ=255)
                                                                                                                            agggCCCATGGTTGCGGGCTGGCAACGGTGGTTTGTGCCGCGCCAGCCGTGCGCAGCTCAGAGCGATCTTTGCCTGGCTTAGGAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGTCGGATCCCTCTGGCGCATaa            <  2:255063/141‑1 (MQ=255)
                                                                                                                            aaGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCACTCTGGGGCATaa            >  1:682430/1‑143 (MQ=255)
                                                                                                                                 ccagccgtGCTGGCTAGAAACGTTGGCTTGGGCCGCCCCAGCGTTGCGCAGCTCATAGCGATCTTTGCCTGGCTGATCAACAGCCATATCAGGGTGCGTATGCATAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCg       <  2:468720/136‑1 (MQ=255)
                                                                                                                                  cATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCgg      <  1:121020/143‑1 (MQ=255)
                                                                                                                                  cATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCCTAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGAACCCTCTGGAGCAAAATGCCgg      >  1:11621/1‑143 (MQ=255)
                                                                                                                                   aTCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATCATGCCggg     >  1:405539/1‑143 (MQ=255)
                                                                                                                                    tCGTTGCTGGCTGGCAACGATGTTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGa    <  2:9979/143‑1 (MQ=255)
                                                                                                                                     cGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCTCAGCTCATAGCTATCTTCGAGTAAAACGCCAACAGCAGCAACTTGAAAAAGTGGCTCAGTCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGAt   <  1:515417/143‑1 (MQ=255)
                                                                                                                                      gTTGCTGGCTGGCAACGATGGTTTGTGCCGCACCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCATATGCTTAATCAGCCATGAACGGATCCCTCTGGCGCATAATACCGGGATc  >  2:28281/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                        
CACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATC  >  NC_000913/4041077‑4041353

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