Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,838,693 A→G I280T (ATC→ACC)  yicM ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,838,6930AG99.6% 1380.1 / ‑4.2 440I280T (ATC→ACC) yicMputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (186/252);  total (188/252)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 451 total aligned reads displayed.

GTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGATGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGC  >  NC_000913/3838603‑3838782
                                                                                          |                                                                                         
gTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTg                                                                                         >  2:1472385/1‑93 (MQ=255)
 ttAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTgg                                                                                        >  1:4890072/1‑93 (MQ=255)
    aGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGa                                                                                     <  1:525396/93‑1 (MQ=255)
     gCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGat                                                                                    >  2:6140421/1‑93 (MQ=255)
       cATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatca                                                                                  <  2:2022599/93‑1 (MQ=255)
         tCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatcatc                                                                                >  1:3319361/1‑93 (MQ=255)
           cACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATccc                                                                              <  1:4047855/93‑1 (MQ=255)
            aCGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCt                                                                             <  1:4672083/93‑1 (MQ=255)
             cGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTg                                                                            <  1:2676557/93‑1 (MQ=255)
               ccGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGcc                                                                          >  2:6906374/1‑93 (MQ=255)
                cGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCa                                                                         >  1:2104799/1‑93 (MQ=255)
                 gAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCAt                                                                        >  1:4330030/1‑93 (MQ=255)
                  aaTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATc                                                                       <  2:5982678/93‑1 (MQ=255)
                    tCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATcac                                                                     >  2:3917191/1‑93 (MQ=255)
                    tCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATcac                                                                     >  2:5822114/1‑93 (MQ=255)
                      ccGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCacac                                                                   >  1:4757291/1‑93 (MQ=255)
                       cGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACAcc                                                                  <  1:3670588/93‑1 (MQ=255)
                       cGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCAACACAcc                                                                  <  2:4168774/93‑1 (MQ=255)
                        gCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACAccc                                                                 <  2:2766379/93‑1 (MQ=255)
                        gCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACAccc                                                                 >  2:3440673/1‑93 (MQ=255)
                         ccAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCg                                                                <  2:553868/93‑1 (MQ=255)
                          cAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCgg                                                               >  1:7439342/1‑93 (MQ=255)
                           aGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGa                                                              <  2:7378802/93‑1 (MQ=255)
                           aGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGa                                                              <  1:3822206/93‑1 (MQ=255)
                           aGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGa                                                              <  1:2015466/93‑1 (MQ=255)
                            ggTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGAc                                                             >  2:3223438/1‑93 (MQ=255)
                              ttCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGt                                                           >  2:2545021/1‑93 (MQ=255)
                              ttCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGt                                                           >  2:7248869/1‑93 (MQ=255)
                               tCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGtt                                                          <  1:4953968/93‑1 (MQ=255)
                               tCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGtt                                                          <  1:5414477/93‑1 (MQ=255)
                                cATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTg                                                         >  2:671566/1‑93 (MQ=255)
                                  tataCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCNNCCATCACACCCGGACGTTNta                                                       >  1:1035071/1‑93 (MQ=255)
                                    taCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaat                                                     >  1:7849820/1‑93 (MQ=255)
                                     acacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaata                                                    <  1:6838629/93‑1 (MQ=255)
                                     acacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaata                                                    <  1:4113722/93‑1 (MQ=255)
                                      cacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaataa                                                   <  1:6418092/93‑1 (MQ=255)
                                      cacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaataa                                                   <  1:2737974/93‑1 (MQ=255)
                                      cacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaataa                                                   <  2:5671732/93‑1 (MQ=255)
                                        ccGGGCTAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGc                                                 <  2:3702744/93‑1 (MQ=255)
                                          gggCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCgg                                               <  1:1358976/93‑1 (MQ=255)
                                             cGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGaaa                                            <  1:1773882/93‑1 (MQ=255)
                                              gAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAg                                           >  1:1955086/1‑93 (MQ=255)
                                              gAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAg                                           >  1:618923/1‑93 (MQ=255)
                                              gAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAg                                           <  1:3084254/93‑1 (MQ=255)
                                               aaTATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGt                                          <  1:5495337/93‑1 (MQ=255)
                                               aaTATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGt                                          >  2:1101300/1‑93 (MQ=255)
                                                 tataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAt                                        >  1:604618/1‑93 (MQ=255)
                                                  ataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtt                                       <  1:3458728/93‑1 (MQ=255)
                                                  ataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtt                                       <  2:6068361/93‑1 (MQ=255)
                                                   taCGTGAAGAAAGCATACTGCCCGGCGAAAGACATGAAGATGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAttt                                      >  2:7823812/1‑93 (MQ=255)
                                                   taCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAttt                                      >  1:7352432/1‑93 (MQ=255)
                                                   taCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAttt                                      <  1:2486244/93‑1 (MQ=255)
                                                   taCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAttt                                      <  1:1637495/93‑1 (MQ=255)
                                                   taCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAttt                                      <  2:3172849/93‑1 (MQ=255)
                                                   taCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGGGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAttt                                      <  2:649820/93‑1 (MQ=255)
                                                    aCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtttt                                     <  1:3538779/93‑1 (MQ=255)
                                                    aCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtttt                                     <  2:5426324/93‑1 (MQ=255)
                                                    aCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtttt                                     <  2:4440434/93‑1 (MQ=255)
                                                     cGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTg                                    <  1:5079545/93‑1 (MQ=255)
                                                      gTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGt                                   <  2:5179581/93‑1 (MQ=255)
                                                      gTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGt                                   <  1:7534711/93‑1 (MQ=255)
                                                      gTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGt                                   <  2:2403322/93‑1 (MQ=255)
                                                       tGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGtt                                  <  2:6557991/93‑1 (MQ=255)
                                                       tGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGtt                                  >  2:3345790/1‑93 (MQ=255)
                                                       tGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGtt                                  >  1:1175930/1‑93 (MQ=255)
                                                       tGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGtt                                  <  2:286607/93‑1 (MQ=255)
                                                       tGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGtt                                  >  2:16703/1‑93 (MQ=255)
                                                       ggaagaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGtt                                  <  1:2261041/92‑1 (MQ=255)
                                                        gaagaaAGCAAACTGCCAGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGttt                                 <  1:7041872/93‑1 (MQ=255)
                                                          agaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCt                               <  1:5808943/93‑1 (MQ=255)
                                                           gaaAGCAAACTGCCCGGCGAAAGACATGAAGGTAGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTg                              <  2:6357203/93‑1 (MQ=255)
                                                            aaaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGa                             <  1:375245/93‑1 (MQ=255)
                                                            aaaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGa                             <  2:1812118/93‑1 (MQ=255)
                                                            aaaGCAAACCGCCCGGCGAAAGACATGAAGGTGGCGGTTATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGa                             <  2:261820/93‑1 (MQ=255)
                                                             aaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGAt                            <  2:4618223/93‑1 (MQ=255)
                                                               acAAACTGCCCGGCGAAAGACATGAAGGTGGCCATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGc                          <  1:1878394/92‑1 (MQ=255)
                                                                 aaaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGa                        <  2:6265031/93‑1 (MQ=255)
                                                                 aaaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGa                        <  1:7148354/93‑1 (MQ=255)
                                                                  aaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAg                       <  2:3093223/93‑1 (MQ=255)
                                                                  aaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAg                       <  2:6992585/93‑1 (MQ=255)
                                                                   aCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAgg                      <  2:1915441/93‑1 (MQ=255)
                                                                    cTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAggg                     <  2:2235814/93‑1 (MQ=255)
                                                                     tGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGt                    <  1:84865/93‑1 (MQ=255)
                                                                     tGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGt                    >  2:3661394/1‑93 (MQ=255)
                                                                       cccGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTc                  >  2:4295274/1‑93 (MQ=255)
                                                                        cgggCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCg                 <  1:4569514/91‑1 (MQ=255)
                                                                        ccGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCg                 >  2:5277960/1‑93 (MQ=255)
                                                                        ccGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCg                 <  2:1301683/93‑1 (MQ=255)
                                                                         cGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGc                <  2:4145286/93‑1 (MQ=255)
                                                                          ggCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGcc               <  1:5853643/93‑1 (MQ=255)
                                                                          ggCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGcc               <  2:6854443/93‑1 (MQ=255)
                                                                           gCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCt              >  1:7703053/1‑93 (MQ=255)
                                                                           gCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCt              <  2:7682320/93‑1 (MQ=255)
                                                                           gCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCt              <  1:2876082/93‑1 (MQ=255)
                                                                              aaaGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGc           <  2:4039661/93‑1 (MQ=255)
                                                                               aaGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCa          <  1:1457750/93‑1 (MQ=255)
                                                                               aaGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCa          >  2:4970100/1‑93 (MQ=255)
                                                                                   cATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGa      <  1:6562415/93‑1 (MQ=255)
                                                                                    aTGAAGGTGGCGATCATCCCTGCCATCTCACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCCGGCAGTGaa     <  2:2939383/93‑1 (MQ=255)
                                                                                       aaGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGc  <  2:3580601/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGATGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGC  >  NC_000913/3838603‑3838782

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