Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,841,014 A→G Y161H (TAC→CAC)  yicO ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,841,0140AG99.5% 1333.6 / ‑3.9 421Y161H (TAC→CAC) yicOhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (210/209);  total (212/209)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

Only 100 of 428 total aligned reads displayed.

TATTTTTTAATCCCATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTACCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCC  >  NC_000913/3840924‑3841103
                                                                                          |                                                                                         
tATTTTTTAATCCCATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGcc                                                                                         >  1:7845740/1‑93 (MQ=255)
   tttttAATCCCATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGa                                                                                      <  1:5254586/93‑1 (MQ=255)
    ttttAATCCCATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGAt                                                                                     >  2:5105580/1‑93 (MQ=255)
     tttAATCCCATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATa                                                                                    <  1:6383414/93‑1 (MQ=255)
       tAATCCCATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATACAGTGCCGGATACg                                                                                  >  1:5353820/1‑93 (MQ=255)
         aTCCCATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGaa                                                                                >  2:870779/1‑93 (MQ=255)
          tCCCATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGaaa                                                                               <  1:6789808/93‑1 (MQ=255)
          tCCCATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGaaa                                                                               <  1:3883646/93‑1 (MQ=255)
          tCCCATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGaaa                                                                               <  2:5544440/93‑1 (MQ=255)
           cccATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAAc                                                                              >  1:457457/1‑93 (MQ=255)
           cccATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAAc                                                                              <  1:6911805/93‑1 (MQ=255)
           cccATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAAc                                                                              <  1:2935869/93‑1 (MQ=255)
            ccATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACa                                                                             <  2:5881117/93‑1 (MQ=255)
            ccATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACa                                                                             <  2:4543949/93‑1 (MQ=255)
            ccATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACa                                                                             <  1:4260783/93‑1 (MQ=255)
            ccATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACa                                                                             <  1:1653428/93‑1 (MQ=255)
            ccATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACa                                                                             <  1:5058006/93‑1 (MQ=255)
              aTTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGc                                                                           >  1:1015796/1‑93 (MQ=255)
              aTTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGc                                                                           >  1:6211276/1‑93 (MQ=255)
              aTTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGc                                                                           >  1:555935/1‑93 (MQ=255)
              aTTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGc                                                                           >  1:2084861/1‑93 (MQ=255)
              aTTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGc                                                                           <  2:6759998/93‑1 (MQ=255)
              aTTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGc                                                                           <  2:6830881/93‑1 (MQ=255)
               ttAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCg                                                                          <  2:4054613/93‑1 (MQ=255)
               ttAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCg                                                                          <  2:6579600/93‑1 (MQ=255)
               ttAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCg                                                                          >  2:6641680/1‑93 (MQ=255)
               ttAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCg                                                                          >  2:7128810/1‑93 (MQ=255)
                tAAGGCGATAAATAATCTAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGt                                                                         <  2:2383882/93‑1 (MQ=255)
                tAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGt                                                                         >  1:6716792/1‑93 (MQ=255)
                tAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGt                                                                         <  1:1145568/93‑1 (MQ=255)
                 aaGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTg                                                                        >  2:5253512/1‑93 (MQ=255)
                  aGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGa                                                                       >  2:4868692/1‑93 (MQ=255)
                   ggCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAg                                                                      >  2:6386081/1‑93 (MQ=255)
                      gATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCaa                                                                   >  2:2905673/1‑93 (MQ=255)
                      gATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCaa                                                                   >  2:5352213/1‑93 (MQ=255)
                       ataaataaTCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCaaa                                                                  >  1:7226202/1‑93 (MQ=255)
                       ataaataaTCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCaaa                                                                  <  1:3793875/93‑1 (MQ=255)
                        taaataaTCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCaaaa                                                                 >  1:17803/1‑93 (MQ=255)
                          aataatCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAAt                                                               <  2:3943644/93‑1 (MQ=255)
                          aataatCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAAt                                                               >  1:7098070/1‑93 (MQ=255)
                           ataatCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATa                                                              <  2:3124571/93‑1 (MQ=255)
                           ataatCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATa                                                              >  2:2132993/1‑93 (MQ=255)
                            taatCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAg                                                             >  1:5867948/1‑93 (MQ=255)
                            taatCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAg                                                             <  2:5933661/93‑1 (MQ=255)
                            taatCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAg                                                             >  2:7589404/1‑93 (MQ=255)
                             aatCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGt                                                            <  1:6496144/93‑1 (MQ=255)
                               tCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTcc                                                          <  1:1271445/93‑1 (MQ=255)
                                ccAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCa                                                         >  2:358945/1‑93 (MQ=255)
                                 cAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCaa                                                        >  1:7156946/1‑93 (MQ=255)
                                 cAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCaa                                                        <  2:5775524/93‑1 (MQ=255)
                                  aaTTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAAc                                                       <  2:1245072/93‑1 (MQ=255)
                                  aaTTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAAc                                                       >  2:543176/1‑93 (MQ=255)
                                   aTTGCGCTGGTGATACCAATACGTAAACGTAAGGGAATGTTGGAGATCATCCAGTGCTGGATACGAAACAGCGTGAGCAAAAATAGTCCAACt                                                      <  1:1363391/93‑1 (MQ=255)
                                   aTTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACt                                                      <  2:7866745/93‑1 (MQ=255)
                                   aTTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACt                                                      <  1:1575030/93‑1 (MQ=255)
                                    ttCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTg                                                     >  2:3840555/1‑93 (MQ=255)
                                    ttCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTg                                                     >  1:5759341/1‑93 (MQ=255)
                                          tGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGcccc                                               <  2:1426790/93‑1 (MQ=255)
                                          tGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGcccc                                               >  1:6618987/1‑93 (MQ=255)
                                          tGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGcccc                                               >  1:3541887/1‑93 (MQ=255)
                                           ggTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCa                                              <  1:3106137/93‑1 (MQ=255)
                                              gATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGaa                                           <  2:3495077/93‑1 (MQ=255)
                                                tACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAata                                         <  1:2695627/93‑1 (MQ=255)
                                                tACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAata                                         >  1:4003267/1‑93 (MQ=255)
                                                tACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAata                                         >  1:3030432/1‑93 (MQ=255)
                                                 aCCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAatat                                        <  1:2515512/93‑1 (MQ=255)
                                                 aCCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAatat                                        <  1:3712581/93‑1 (MQ=255)
                                                  ccAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATAtt                                       >  1:7700551/1‑93 (MQ=255)
                                                   cAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTg                                      <  2:4466680/93‑1 (MQ=255)
                                                      tACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTgcgc                                   <  1:896157/93‑1 (MQ=255)
                                                      tACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTgcgc                                   <  2:1722287/93‑1 (MQ=255)
                                                       aCGTAAACTTAAGGGAATGTTGGAGATCATCCAGTACCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGcc                                  >  2:5027605/1‑93 (MQ=255)
                                                              cTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCg                           >  2:6482386/1‑93 (MQ=255)
                                                              cTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCg                           <  1:441287/93‑1 (MQ=255)
                                                              cTTAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCg                           <  1:3296095/93‑1 (MQ=255)
                                                               ttAAGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCgg                          >  2:1106484/1‑93 (MQ=255)
                                                                  aGGGAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCt                       >  2:6281628/1‑93 (MQ=255)
                                                                   gggAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTg                      >  2:2823723/1‑93 (MQ=255)
                                                                    ggAATGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGc                     <  1:2271201/93‑1 (MQ=255)
                                                                       aTGTTGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAg                  <  1:988269/93‑1 (MQ=255)
                                                                           tGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGgaga              <  1:6027911/93‑1 (MQ=255)
                                                                           tGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGgaga              <  2:2604165/93‑1 (MQ=255)
                                                                           tGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGgaga              <  2:3660815/93‑1 (MQ=255)
                                                                           tGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGgaga              <  1:4481934/93‑1 (MQ=255)
                                                                           tGGAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGgaga              <  1:5450796/93‑1 (MQ=255)
                                                                            ggAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGgagag             >  2:1874911/1‑92 (MQ=255)
                                                                            ggAGATCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGAt             >  2:2692446/1‑93 (MQ=255)
                                                                             gagaTCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATg            <  2:2725936/93‑1 (MQ=255)
                                                                              agaTCATCCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGc           >  1:1523345/1‑93 (MQ=255)
                                                                                 tcatcCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCa        <  2:1526691/93‑1 (MQ=255)
                                                                                 tcatcCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCa        <  2:4662090/93‑1 (MQ=255)
                                                                                 tcatcCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCa        <  1:4135740/93‑1 (MQ=255)
                                                                                 acatcCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCa        <  1:7490494/92‑1 (MQ=255)
                                                                                  catcCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCAt       >  1:4702847/1‑93 (MQ=255)
                                                                                  catcCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCAt       <  1:3364123/93‑1 (MQ=255)
                                                                                   atcCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATc      >  1:1978457/1‑93 (MQ=255)
                                                                                   atcCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATc      <  2:6184264/93‑1 (MQ=255)
                                                                                    tcCAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCTGTCTGCCAGGAGATGCCCATCg     >  2:7681216/1‑93 (MQ=255)
                                                                                     ccAGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGc    >  2:4318544/1‑93 (MQ=255)
                                                                                       aGTGCCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGccc  <  2:729239/93‑1 (MQ=255)
                                                                                          |                                                                                         
TATTTTTTAATCCCATTAAGGCGATAAATAATCCAATTCCGCTGGTGATACCAATACGTAAACTTAAGGGAATGTTGGAGATCATCCAGTACCGGATACGAAACAGCGTGAGCAAAAATAGTCCAACTGCGCCCCAGAATATTGCGCCCATCCCGGTCTGCCAGGAGATGCCCATCGCCC  >  NC_000913/3840924‑3841103

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