Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,484,849 G→A R355H (CGC→CAC)  yjgP → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,484,8490GA99.5% 696.6 / NA 203R355H (CGC→CAC) yjgPhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (101/101);  minor base C (1/0);  total (102/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 232 total aligned reads displayed.

GACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGTCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTCACCAC  >  NC_000913/4484761‑4484938
                                                                                        |                                                                                         
gACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTg                                                                                       <  1:704750/93‑1 (MQ=255)
 gcccGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTgc                                                                                      <  2:2774145/92‑1 (MQ=255)
 aCCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTgc                                                                                      <  1:1929538/93‑1 (MQ=255)
      aCGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGcc                                                                                 >  2:1104349/1‑93 (MQ=255)
       cGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCa                                                                                <  2:1786214/93‑1 (MQ=255)
       cGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCa                                                                                >  2:2189930/1‑93 (MQ=255)
       cGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCa                                                                                <  1:2123100/93‑1 (MQ=255)
        gCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAg                                                                               <  1:481378/93‑1 (MQ=255)
        gCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAg                                                                               <  2:4373266/93‑1 (MQ=255)
        gCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAg                                                                               >  2:3884317/1‑93 (MQ=255)
         cTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGt                                                                              <  1:4406480/93‑1 (MQ=255)
          tgtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGtt                                                                             >  2:4405449/1‑93 (MQ=255)
          tgtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGtt                                                                             <  2:1031418/93‑1 (MQ=255)
           gtggatgtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGttt                                                                            >  2:1471242/1‑93 (MQ=255)
                 gtggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTcgcg                                                                      <  1:2909414/93‑1 (MQ=255)
                   ggaCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTa                                                                    >  2:3110311/1‑93 (MQ=255)
                     aCCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTaaa                                                                  <  2:1738825/93‑1 (MQ=255)
                       cGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAgg                                                                >  2:3524043/1‑93 (MQ=255)
                        gTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGa                                                               >  1:242871/1‑93 (MQ=255)
                        gTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGa                                                               <  2:2490708/93‑1 (MQ=255)
                         ttAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAg                                                              >  1:2528771/1‑93 (MQ=255)
                          tAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGc                                                             <  2:3559677/93‑1 (MQ=255)
                          tAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGc                                                             <  2:2679797/93‑1 (MQ=255)
                          tAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGc                                                             >  1:3192762/1‑93 (MQ=255)
                             ccTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCCCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCggg                                                          >  2:4320179/1‑92 (MQ=255)
                             ccTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgt                                                          <  1:281353/93‑1 (MQ=255)
                              cTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgtg                                                         >  1:1980732/1‑93 (MQ=255)
                              cTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgtg                                                         >  1:550947/1‑93 (MQ=255)
                              cTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgtg                                                         >  1:3137553/1‑93 (MQ=255)
                               tGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGgtgg                                                        >  2:1625917/1‑92 (MQ=255)
                                gATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGtgtg                                                       <  1:2833420/93‑1 (MQ=255)
                                gATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGtgtg                                                       >  2:1859939/1‑93 (MQ=255)
                                gATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGtgtg                                                       >  2:3232571/1‑93 (MQ=255)
                                gATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGtgtg                                                       >  2:3850951/1‑93 (MQ=255)
                                gATTTATCTGGATTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGTGGtgtg                                                       >  2:2382468/1‑93 (MQ=255)
                                 aTTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGa                                                      <  1:907212/93‑1 (MQ=255)
                                 aTTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGa                                                      <  2:926212/93‑1 (MQ=255)
                                  tttATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGAt                                                     <  1:4643300/93‑1 (MQ=255)
                                  tttATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGAt                                                     >  2:2969024/1‑93 (MQ=255)
                                  tttATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGAt                                                     >  1:3491258/1‑93 (MQ=255)
                                    tATCTGGCTTTAGCGATTGTTATCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGc                                                   >  1:1029305/1‑93 (MQ=255)
                                     aTCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCa                                                  >  1:4185172/1‑93 (MQ=255)
                                     aTCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCa                                                  >  1:4593706/1‑93 (MQ=255)
                                     aTCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCa                                                  >  2:4675802/1‑93 (MQ=255)
                                      tCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCaa                                                 >  1:4540604/1‑93 (MQ=255)
                                      tCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCaa                                                 >  2:4452775/1‑93 (MQ=255)
                                       cTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAAc                                                >  2:242856/1‑93 (MQ=255)
                                       cTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAAc                                                <  1:85707/93‑1 (MQ=255)
                                        tGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGTCACCTGCGCGCCAGTGTTTCGCGTAAAGGAGCGGTGTGATGCAAcc                                               >  1:651468/1‑93 (MQ=255)
                                        tGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAAcc                                               >  1:2893142/1‑93 (MQ=255)
                                         ggCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCt                                              >  2:484816/1‑93 (MQ=255)
                                          gCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtt                                             <  2:95257/93‑1 (MQ=255)
                                          gCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtt                                             <  2:3156552/93‑1 (MQ=255)
                                           cTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCttt                                            <  2:4250860/93‑1 (MQ=255)
                                           cTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCttt                                            <  1:3935760/93‑1 (MQ=255)
                                           cTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCttt                                            <  1:1247372/93‑1 (MQ=255)
                                            tttAGCGATTGTTCTCAGCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtttt                                           <  1:2208474/93‑1 (MQ=255)
                                            tttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGTAACCtttt                                           >  2:1397704/1‑93 (MQ=255)
                                            tttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCtttt                                           <  2:2758908/93‑1 (MQ=255)
                                             ttAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTg                                          <  2:2993321/93‑1 (MQ=255)
                                               aGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGc                                        <  2:3629482/93‑1 (MQ=255)
                                                gCGATTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCg                                       <  1:2156283/93‑1 (MQ=255)
                                                  gATTGTTCTCAACCTTTGGGACACCGTTCCGGCACGCTTCCCGCGGGCCAGTTTTTCGCGTAAAGGAGCGGTGTGTTGCAACCTTTTGGCGTa                                     <  2:772325/93‑1 (MQ=255)
                                                   aTTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTAc                                    <  1:4160173/93‑1 (MQ=255)
                                                   aTTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTAc                                    >  2:2311516/1‑93 (MQ=255)
                                                   aTTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTAc                                    <  1:4488125/93‑1 (MQ=255)
                                                   aTTGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTAc                                    >  2:2033422/1‑93 (MQ=255)
                                                    ttGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACt                                   >  2:3272016/1‑93 (MQ=255)
                                                     tGTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACtt                                  <  2:2808913/93‑1 (MQ=255)
                                                      gTTCTCAATCTTTGGGACACAGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTg                                 >  2:4583597/1‑93 (MQ=255)
                                                      gTTCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTg                                 >  2:120087/1‑93 (MQ=255)
                                                       ttCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGa                                <  1:143488/93‑1 (MQ=255)
                                                       ttCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGa                                >  2:1307212/1‑93 (MQ=255)
                                                       ttCTCAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGa                                <  1:3055270/93‑1 (MQ=255)
                                                        tctcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGAc                               <  2:3666565/93‑1 (MQ=255)
                                                          tcAACCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCg                             >  1:3615219/1‑93 (MQ=255)
                                                             aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                          <  1:1079015/93‑1 (MQ=255)
                                                             aCCTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCta                          <  2:2587919/93‑1 (MQ=255)
                                                              ccTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCtat                         >  2:3397492/1‑93 (MQ=255)
                                                              ccTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCtat                         <  1:4051142/93‑1 (MQ=255)
                                                               cTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCtata                        <  1:2412899/93‑1 (MQ=255)
                                                               cTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCtata                        <  2:4356263/93‑1 (MQ=255)
                                                               cTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCtata                        >  2:2845829/1‑93 (MQ=255)
                                                               cTTTGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCtata                        >  1:3291226/1‑93 (MQ=255)
                                                                tttGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTatat                       >  1:670623/1‑93 (MQ=255)
                                                                tttGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTatat                       <  1:1171248/93‑1 (MQ=255)
                                                                tttGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTatat                       <  1:3202826/93‑1 (MQ=255)
                                                                  tGGGACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCg                     >  1:3988842/1‑93 (MQ=255)
                                                                   gggACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCgg                    <  1:2286519/93‑1 (MQ=255)
                                                                    ggACACCGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGt                   <  1:3178008/93‑1 (MQ=255)
                                                                         ccGTGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAAc              >  2:3970892/1‑93 (MQ=255)
                                                                            tGCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTAt           <  1:2772452/93‑1 (MQ=255)
                                                                             gCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTAtt          <  1:183910/93‑1 (MQ=255)
                                                                             gCCGGCCCGCCACCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTAtt          <  1:2330931/93‑1 (MQ=255)
                                                                                ggcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTc       <  2:2074851/90‑1 (MQ=255)
                                                                                ggcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTc       <  1:2562844/90‑1 (MQ=255)
                                                                                ggcccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGGGTGATGCAACATTTTGGCGTACTTGACCGCTATATCGGTAAAACTAtttta       >  2:361457/4‑92 (MQ=255)
                                                                                  cccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTcac     <  1:1762863/92‑1 (MQ=255)
                                                                                   ccgccaCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTcacc    <  1:1279724/93‑1 (MQ=255)
                                                                                    cgccacCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTcacca   <  1:4220106/93‑1 (MQ=255)
                                                                                        |                                                                                         
GACCCGACGCTGTGGATGTGGACCGTTAACCTGATTTATCTGGCTTTAGCGATTGTTCTCAACCTTTGGGACACCGTGCCGGTCCGCCGCCTGCGCGCCAGTTTTTCGCGTAAAGGAGCGGTGTGATGCAACCTTTTGGCGTACTTGACCGCTATATCGGTAAAACTATTTTCACCAC  >  NC_000913/4484761‑4484938

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