Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,624,700 T→C I807V (ATT→GTT)  yhiH ← putative ATP‑binding component of a transport system, fragment 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,624,7000TC99.7% 1215.8 / ‑5.2 388I807V (ATT→GTT) yhiHputative ATP‑binding component of a transport system, fragment 1
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (183/204);  total (183/205)
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.96e-01

Only 100 of 403 total aligned reads displayed.

GCCCTTCCAGCGAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAATCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGTCAGGAAACTG  >  NC_000913/3624610‑3624791
                                                                                          |                                                                                           
gCCCTTCCAGCGAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCt                                                                                           >  1:5097336/1‑93 (MQ=255)
 cccTTCCAGCGAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTg                                                                                          <  2:4645182/93‑1 (MQ=255)
  ccTTCCAGCGAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTgg                                                                                         <  2:4629612/93‑1 (MQ=255)
  ccTTCCAGCGAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTgg                                                                                         >  2:4049039/1‑93 (MQ=255)
    ttCCAGCGAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCt                                                                                       >  2:769641/1‑93 (MQ=255)
     tCCAGCGAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCtt                                                                                      >  2:6440010/1‑93 (MQ=255)
         gCGAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCa                                                                                  <  2:7574663/93‑1 (MQ=255)
           gAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATa                                                                                <  1:3919332/93‑1 (MQ=255)
           gAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATa                                                                                <  2:6669599/93‑1 (MQ=255)
           gAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCCGGCTTTTCATa                                                                                >  1:5191108/1‑93 (MQ=255)
           gAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCAGCAACCTGGCTTTTCATa                                                                                <  2:3982627/93‑1 (MQ=255)
            aaGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTTATaa                                                                               >  1:1479067/1‑93 (MQ=255)
            aaGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATaa                                                                               >  1:4124522/1‑93 (MQ=255)
             aGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATcaa                                                                              >  2:4202999/1‑93 (MQ=255)
             aGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATaaa                                                                              <  2:2474285/93‑1 (MQ=255)
               cTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAgg                                                                            >  2:2939942/1‑93 (MQ=255)
                  ccGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTgg                                                                         <  2:1661709/93‑1 (MQ=255)
                  ccGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTgg                                                                         <  1:4362643/93‑1 (MQ=255)
                   cGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGga                                                                        <  1:6297777/93‑1 (MQ=255)
                   cGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGga                                                                        <  2:6065122/93‑1 (MQ=255)
                     gatcgatcATCCCGGAAAACTGTGTCGCCGTGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGTTTTTCATAAAGGTGggga                                                                      <  2:6886000/93‑1 (MQ=255)
                     gatcgatcATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGgaga                                                                      >  2:7476150/1‑93 (MQ=255)
                      atcgatcATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGAt                                                                     <  1:7335960/93‑1 (MQ=255)
                      atcgatcATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGAt                                                                     >  1:3379427/1‑93 (MQ=255)
                       tcgatcATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGAta                                                                    >  1:3579183/1‑92 (MQ=255)
                       tcgatcATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATc                                                                    >  2:6274494/1‑93 (MQ=255)
                       tcgatcATCCCGGAAAACTGTGGCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATc                                                                    <  2:2759237/93‑1 (MQ=255)
                        cgatcATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATca                                                                   <  2:3444181/93‑1 (MQ=255)
                        cgatcATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATca                                                                   <  2:4380724/93‑1 (MQ=255)
                          atcatcCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATcagc                                                                 <  2:2787692/93‑1 (MQ=255)
                            catcCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATcagcag                                                               >  1:4238861/1‑93 (MQ=255)
                             atcCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCagcagc                                                              >  1:3493024/1‑93 (MQ=255)
                              tcCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGcc                                                             >  2:2644021/1‑93 (MQ=255)
                                  ggAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCAt                                                         >  2:4818252/1‑93 (MQ=255)
                                  ggAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGACCCAt                                                         >  1:75647/1‑93 (MQ=255)
                                   gAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCAtt                                                        <  1:5571944/93‑1 (MQ=255)
                                   gAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCAtt                                                        >  1:1151464/1‑93 (MQ=255)
                                   gAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCAtt                                                        >  1:6868198/1‑93 (MQ=255)
                                   gAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCAtt                                                        >  1:1548786/1‑93 (MQ=255)
                                    aaaaCTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTc                                                       <  2:3708137/93‑1 (MQ=255)
                                     aaaCTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTcc                                                      >  2:6570884/1‑93 (MQ=255)
                                     aaaCTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTcc                                                      <  1:2858388/93‑1 (MQ=255)
                                        tggtgtCAGCGAGATCAACGTGATAATCGCCCTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGt                                                   <  1:4745442/91‑1 (MQ=255)
                                        cTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATTAGCAGCCCCATTCCggg                                                   >  1:7036197/1‑92 (MQ=255)
                                        cTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGt                                                   <  1:2403867/93‑1 (MQ=255)
                                             tCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCaa                                              <  2:5441766/93‑1 (MQ=255)
                                             tCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCaa                                              <  1:722071/93‑1 (MQ=255)
                                             tCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCACCAGCCCCATCCCGTGGGCaa                                              >  1:4473317/1‑93 (MQ=255)
                                              cGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAat                                             <  1:1798348/93‑1 (MQ=255)
                                              cGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAat                                             <  2:5133946/93‑1 (MQ=255)
                                              cGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAat                                             >  1:378326/1‑93 (MQ=255)
                                               gCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGGGACAAGCAGCCCCCTCCCGGGGCCAatg                                            >  2:3106345/1‑93 (MQ=255)
                                               gCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAatg                                            >  2:3617534/1‑93 (MQ=255)
                                                ccGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAatga                                           <  1:2329127/93‑1 (MQ=255)
                                                 cGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAatgat                                          <  1:3143793/93‑1 (MQ=255)
                                                   ggATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAAtgatga                                        >  1:3303433/1‑93 (MQ=255)
                                                   ggATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAAtgatga                                        <  1:4722151/93‑1 (MQ=255)
                                                   ggATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAAtgatga                                        <  2:2092413/93‑1 (MQ=255)
                                                   ggATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCCGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAAtgatga                                        >  2:180443/1‑93 (MQ=255)
                                                     aTCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGata                                      >  1:528594/1‑93 (MQ=255)
                                                     aTCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGata                                      <  1:6415053/93‑1 (MQ=255)
                                                     aTCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGata                                      <  1:1875643/93‑1 (MQ=255)
                                                      tCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGatat                                     >  2:7148405/1‑93 (MQ=255)
                                                      tCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGatat                                     >  2:3148963/1‑93 (MQ=255)
                                                       cAACGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGAtata                                    <  2:696176/93‑1 (MQ=255)
                                                         aCGTGATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATAca                                  <  1:2265649/93‑1 (MQ=255)
                                                            tGATAATCGCCGTGCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATAcagca                               <  2:5582118/93‑1 (MQ=255)
                                                             gATAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATAcagcag                              >  1:1392080/1‑93 (MQ=255)
                                                               tAATCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAgcg                            <  1:6256927/93‑1 (MQ=255)
                                                                aaTCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAgcgc                           <  1:1021086/93‑1 (MQ=255)
                                                                 aTCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGcgcc                          >  1:3256662/1‑93 (MQ=255)
                                                                  tCGCCGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGcgccg                         <  2:4515497/93‑1 (MQ=255)
                                                                     ccGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCa                      <  1:883996/93‑1 (MQ=255)
                                                                     ccGTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATAAAGCAGCgccgccg                      >  1:2004411/1‑92 (MQ=255)
                                                                       gTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAgg                    >  2:118609/1‑93 (MQ=255)
                                                                       gTTCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAgg                    >  1:98479/1‑93 (MQ=255)
                                                                        ttCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAggg                   <  2:3843846/93‑1 (MQ=255)
                                                                        ttCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAggg                   >  2:735568/1‑93 (MQ=255)
                                                                        ttCCGAAAAGGGCGGCAACCTGCCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAggg                   <  2:5805577/93‑1 (MQ=255)
                                                                         tCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCCgggg                  >  1:7844666/1‑92 (MQ=255)
                                                                         tCCGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGt                  <  1:5914149/93‑1 (MQ=255)
                                                                          ccGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTg                 <  1:4952713/93‑1 (MQ=255)
                                                                           cGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGa                >  2:3252739/1‑93 (MQ=255)
                                                                           cGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGa                <  1:2617426/93‑1 (MQ=255)
                                                                           cGAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGa                <  2:4723276/93‑1 (MQ=255)
                                                                            gAAAATGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAg               <  1:4609069/93‑1 (MQ=255)
                                                                            gAAAATGGCGGCAACCTGGCTTGTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGAGGATAGACAGCAGCGCCGCCAGGGTGAg               <  2:7935204/93‑1 (MQ=255)
                                                                             aaaaTGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGc              <  2:2005352/93‑1 (MQ=255)
                                                                             aaaaTGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGc              >  2:1360992/1‑93 (MQ=255)
                                                                             aaaaTGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGc              >  2:6220794/1‑93 (MQ=255)
                                                                               aaTGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGt            >  1:4886541/1‑93 (MQ=255)
                                                                                aTGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGTc           <  1:7520617/93‑1 (MQ=255)
                                                                                 tGGCGGCAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGTCa          <  2:6368755/93‑1 (MQ=255)
                                                                                  ggcggcAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGTCAg         <  1:6116765/93‑1 (MQ=255)
                                                                                   gcggcAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGTCAgg        <  2:3890684/93‑1 (MQ=255)
                                                                                    cggcAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGTCAGGa       <  1:2781005/93‑1 (MQ=255)
                                                                                     ggcAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGTCAGGaa      <  1:1234935/93‑1 (MQ=255)
                                                                                      gcAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGTCAGGaaa     >  1:7448667/1‑93 (MQ=255)
                                                                                       cAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGTCAGGAAAc    <  1:6631273/93‑1 (MQ=255)
                                                                                       cAACCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGTCAGGAAAc    <  1:3676597/93‑1 (MQ=255)
                                                                                          |                                                                                           
GCCCTTCCAGCGAAGCTACCGGATCGATCATCCCGGAAAACTGTGTCGCCGGGATCAACGTGATAATCGCCGTTCCGAAAATGGCGGCAATCTGGCTTTTCATAAAGGTGGAGATCAGCAGCCCCATTCCGGTGGCAATGATGATATACAGCAGCGCCGCCAGGGTGAGCGTCAGGAAACTG  >  NC_000913/3624610‑3624791

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