Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,041,221 G→A 7.0% A71A (GCC→GCT)  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,2210GA 7.0% 175.3 / 3.4 69A71A (GCC→GCT) mobBmolybdopterin‑guanine dinucleotide biosynthesis adaptor protein
Reads supporting (aligned to +/- strand):  ref base G (25/39);  new base A (2/3);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

AATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTC  >  NC_000913/4041080‑4041362
                                                                                                                                             |                                                                                                                                             
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)
                                                                                                                                           ctggcAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttttt   <  1:266759/143‑1 (MQ=255)
                                                                                                                                            tggcAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCGATCTTTGCCTGGCCTATCAACATCCGTATCATGGTGCGTGTGCGTAATCAGCCATGCACGGAGCCCTCTGGCGCATGATGCCGAGCTCGgttgttgg  >  2:883712/1‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                             
AATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTC  >  NC_000913/4041080‑4041362

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