Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 75,712 C→T 5.7% G530S (GGC→AGC) ‡ sgrR ← transcriptional DNA‑binding transcriptional activator of sgrS sRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,7120CT 5.7% 292.6 / 6.5 88G530S (GGC→AGC) ‡sgrRtranscriptional DNA‑binding transcriptional activator of sgrS sRNA
Reads supporting (aligned to +/- strand):  ref base C (42/41);  new base T (2/3);  total (44/44)
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.21e-01

GCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAAT  >  NC_000913/75578‑75853
                                                                                                                                      |                                                                                                                                              
gcACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGt                                                                                                                                  <  2:610250/148‑1 (MQ=255)
               gCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGt                                                                                                                                  <  2:993847/133‑1 (MQ=255)
               gCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGt                                                                                                                                  <  4:859383/133‑1 (MQ=255)
               gCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGt                                                                                                                                  >  1:993847/1‑133 (MQ=255)
               gCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGt                                                                                                                                  >  3:859383/1‑133 (MQ=255)
                          agaACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCa                                                                                                              <  3:1135607/142‑1 (MQ=255)
                          agaACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCa                                                                                                              >  4:1135607/1‑142 (MQ=255)
                           gaACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCGGCCTGCGCAt                                                                                                                                         <  3:631386/114‑1 (MQ=255)
                           gaACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCGGCCTGCGCAt                                                                                                                                         >  4:631386/1‑114 (MQ=255)
                            aaCGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAgt                                                                                                            <  2:294915/142‑1 (MQ=255)
                            aaCGGCGTTAGCGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAgt                                                                                                            >  1:294915/1‑142 (MQ=255)
                                           gtgATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCg                                                                                       >  1:931906/1‑148 (MQ=255)
                                            tgATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGcc                                                                                                               <  4:195099/123‑1 (MQ=255)
                                            tgATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGcc                                                                                                               >  3:195099/1‑123 (MQ=255)
                                                tttGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaataa                                                                                                  >  3:414596/1‑132 (MQ=255)
                                                tttGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaataa                                                                                                  <  4:414596/132‑1 (MQ=255)
                                                  tGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCt                                                                                <  2:931906/148‑1 (MQ=255)
                                                    ttATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCCGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaataat                                                                                                 >  1:261471/1‑129 (MQ=255)
                                                    ttATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaataat                                                                                                 >  3:723354/1‑129 (MQ=255)
                                                    ttATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaataat                                                                                                 <  4:723354/129‑1 (MQ=255)
                                                    ttATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATCCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaataat                                                                                                 <  2:261471/129‑1 (MQ=255)
                                                     tATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGt                                                                                                                                  <  1:15377/95‑1 (MQ=255)
                                                     tATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGt                                                                                                                                  >  2:15377/1‑95 (MQ=255)
                                                     tATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGtaataa                                                                                                  >  3:810924/1‑127 (MQ=255)
                                                     tATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGtaataa                                                                                                  <  4:810924/127‑1 (MQ=255)
                                                          aCCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTg                                                                                                                                >  1:767301/1‑92 (MQ=255)
                                                          aCCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTg                                                                                                                                <  2:767301/92‑1 (MQ=255)
                                                            cAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCt                                                                                                                           >  4:56983/1‑95 (MQ=255)
                                                            cAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCt                                                                                                                           <  3:56983/95‑1 (MQ=255)
                                                            cAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAg                                                                                                             <  3:441092/109‑1 (MQ=255)
                                                            cAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAg                                                                                                             >  4:441092/1‑109 (MQ=255)
                                                            cAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTg                                                                                  >  4:787266/1‑136 (MQ=255)
                                                            cAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTg                                                                                  <  3:787266/136‑1 (MQ=255)
                                                             aGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATgg                                                                                               <  2:1063566/122‑1 (MQ=255)
                                                             aGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATgg                                                                                               >  1:1063566/1‑122 (MQ=255)
                                                              gCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCtt                                                                                    <  3:451525/132‑1 (MQ=255)
                                                              gCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCtt                                                                                    >  4:451525/1‑132 (MQ=255)
                                                               cAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAGATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAg                                                                                                             <  3:793449/106‑1 (MQ=255)
                                                               cAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAg                                                                                                             >  4:793449/1‑106 (MQ=255)
                                                                   cATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGAt                                                                                                                    >  2:1019641/1‑95 (MQ=255)
                                                                   cATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGAt                                                                                                                    <  1:1019641/95‑1 (MQ=255)
                                                                   cATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGt                                                                      <  1:1000901/141‑1 (MQ=255)
                                                                   cATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGt                                                                      >  2:1000901/1‑141 (MQ=255)
                                                                   cATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGc                                                               >  3:609642/1‑148 (MQ=255)
                                                                     tGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaa                                                                                                     >  1:668415/1‑108 (MQ=255)
                                                                     tGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaa                                                                                                     <  2:668415/108‑1 (MQ=255)
                                                                     tGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGTCCAGt                                                                      <  4:419968/139‑1 (MQ=255)
                                                                     tGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGt                                                                      >  3:419968/1‑139 (MQ=255)
                                                                     tGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTg                                                                    >  2:735983/1‑141 (MQ=255)
                                                                     tGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTg                                                                    <  1:735983/141‑1 (MQ=255)
                                                                          ccGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGa                                                                                                                     <  2:728537/87‑1 (MQ=255)
                                                                          ccGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGa                                                                                                                     >  1:728537/1‑87 (MQ=255)
                                                                              tGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGcc                                                    <  4:609642/148‑1 (MQ=255)
                                                                                            cTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGc                                                                                              >  2:382998/1‑92 (MQ=255)
                                                                                            cTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGc                                                                                              <  3:523376/92‑1 (MQ=255)
                                                                                            cTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGc                                                                                              >  4:523376/1‑92 (MQ=255)
                                                                                            cTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGc                                                                                              <  1:382998/92‑1 (MQ=255)
                                                                                            cTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCg                                                                            >  3:929491/1‑110 (MQ=255)
                                                                                            cTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCg                                                                            <  4:929491/110‑1 (MQ=255)
                                                                                                    aaaTCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTgcg                                >  3:874501/1‑146 (MQ=255)
                                                                                                    aaaTCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTgcg                                <  4:874501/146‑1 (MQ=255)
                                                                                                     aaTCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGcc                                                    >  1:345820/1‑125 (MQ=255)
                                                                                                     aaTCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGcc                                                    <  2:345820/125‑1 (MQ=255)
                                                                                                      aTCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGaa                                                                                                                        <  4:373137/56‑1 (MQ=255)
                                                                                                      aTCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGaa                                                                                                                        >  3:373137/1‑56 (MQ=255)
                                                                                                       tCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTg                                                                                  <  1:1054416/93‑1 (MQ=255)
                                                                                                       tCGAACCAGCAGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTg                                                                                  >  2:1054416/1‑93 (MQ=255)
                                                                                                        cGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCttt                                                                                   >  1:774854/1‑91 (MQ=255)
                                                                                                        cGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCttt                                                                                   <  2:774854/91‑1 (MQ=255)
                                                                                                        cGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGt                                                                      >  4:377728/1‑104 (MQ=255)
                                                                                                        cGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGt                                                                      <  3:377728/104‑1 (MQ=255)
                                                                                                        cGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGc                                                     <  4:538342/121‑1 (MQ=255)
                                                                                                        cGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGc                                                     >  3:538342/1‑121 (MQ=255)
                                                                                                            ccAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCt                                                                                     >  1:69809/1‑85 (MQ=255)
                                                                                                            ccAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCt                                                                                     <  2:69809/85‑1 (MQ=255)
                                                                                                            ccAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGagca                      >  1:880458/1‑148 (MQ=255)
                                                                                                                   aGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTgg                                                                              >  4:971781/1‑85 (MQ=255)
                                                                                                                   aGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTgg                                                                              <  3:971781/85‑1 (MQ=255)
                                                                                                                         ttCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGcc         >  3:874853/1‑148 (MQ=255)
                                                                                                                         ttCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGcc         >  3:427184/1‑148 (MQ=255)
                                                                                                                           taTGCGCGGCCTGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTg                                                                 >  2:1011503/2‑90 (MQ=255)
                                                                                                                           taTGCGCGGCCTGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTg                                                                 <  1:1011503/89‑1 (MQ=255)
                                                                                                                           taTGCGCGGCCTGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGtt                                                        <  2:288703/98‑1 (MQ=255)
                                                                                                                           taTGCGCGGCCTGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGtt                                                        >  1:288703/2‑99 (MQ=255)
                                                                                                                           taTGCGCGGCCTGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTg           <  3:357121/143‑1 (MQ=255)
                                                                                                                             tGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCg                                                      <  2:143233/99‑1 (MQ=255)
                                                                                                                             tGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCg                                                      >  1:143233/1‑99 (MQ=255)
                                                                                                                                gcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAAt  <  4:427184/148‑1 (MQ=255)
                                                                                                                                      cgcgCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGtt                                                        >  2:968479/1‑88 (MQ=255)
                                                                                                                                      cgcgCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGtt                                                        <  1:968479/88‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                              
GCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAAT  >  NC_000913/75578‑75853

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: