Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 656,969 A→G W65R (TGG→CGG)  crcB ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913656,9690AG99.1% 713.1 / ‑3.0 221W65R (TGG→CGG) crcBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (116/103);  total (118/103)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

Only 100 of 238 total aligned reads displayed.

CGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCA  >  NC_000913/656878‑657058
                                                                                           |                                                                                         
 ggCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCcgca                                                                                         >  1:4565404/1‑93 (MQ=255)
 ggCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCcgca                                                                                         >  1:874509/1‑93 (MQ=255)
   cccTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACt                                                                                       <  1:1327568/93‑1 (MQ=255)
   cccTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACt                                                                                       >  2:20385/1‑93 (MQ=255)
   cccTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACt                                                                                       <  1:2879458/93‑1 (MQ=255)
    ccTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTg                                                                                      >  1:2409644/1‑93 (MQ=255)
     ctctTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTgg                                                                                     >  1:3527413/1‑93 (MQ=255)
      tctTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGa                                                                                    >  2:2757266/1‑93 (MQ=255)
        ttGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGAta                                                                                  >  2:3744355/1‑92 (MQ=255)
        ttGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATc                                                                                  <  1:1062912/93‑1 (MQ=255)
          gTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCaa                                                                                >  1:3405013/1‑93 (MQ=255)
            aaCAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATg                                                                              >  2:3617209/1‑93 (MQ=255)
            aaCAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATg                                                                              >  1:1526704/1‑93 (MQ=255)
             aCAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGt                                                                             <  1:4232779/93‑1 (MQ=255)
             aCAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGt                                                                             >  1:412227/1‑93 (MQ=255)
              cAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGtt                                                                            >  2:1626123/1‑93 (MQ=255)
              cAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGtt                                                                            >  1:1308228/1‑93 (MQ=255)
               aaaaaCACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTc                                                                           >  1:3940860/1‑93 (MQ=255)
               aaaaaCACCACTTCGGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTc                                                                           >  1:273545/1‑93 (MQ=255)
                  aaCACCGCTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTc                                                                        <  2:2994293/93‑1 (MQ=255)
                  aaCACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTc                                                                        >  1:2172551/1‑93 (MQ=255)
                    caccacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCAt                                                                      >  1:326466/1‑93 (MQ=255)
                     accacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATc                                                                     >  2:4050160/1‑93 (MQ=255)
                      ccacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATATTTTCCGCACTGGATCAATGTTCGTCATcc                                                                    >  2:1976513/1‑93 (MQ=255)
                      ccacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATcc                                                                    <  1:2786839/93‑1 (MQ=255)
                      ccacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATcc                                                                    <  2:3772441/93‑1 (MQ=255)
                       cacTTCTGCCGAACATGTTGAGAAGGTTATAAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCTGTGGATCAATGTTTGTCATCct                                                                   <  2:830513/93‑1 (MQ=255)
                       cacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCct                                                                   <  2:3438109/93‑1 (MQ=255)
                       cacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCct                                                                   >  1:3946322/1‑93 (MQ=255)
                          ttCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCctgct                                                                >  2:1291805/1‑93 (MQ=255)
                           tCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCctgctg                                                               >  2:1202209/1‑93 (MQ=255)
                           tCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCctgctg                                                               <  1:2667682/93‑1 (MQ=255)
                             tGCCGAAAAGGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGaa                                                             <  1:2064007/93‑1 (MQ=255)
                              gCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAAc                                                            >  1:3126194/1‑93 (MQ=255)
                                cGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGATCGTCATCCTGCTGAACCa                                                          >  1:3613474/1‑93 (MQ=255)
                                 gAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCAt                                                         >  2:1660646/1‑93 (MQ=255)
                                  aaaaTGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGCTCGTCATCCTGCTGAACCATg                                                        >  1:2281439/1‑93 (MQ=255)
                                  aaaaTGTTGAGAAGGTTGTTAAACCGCCACAAACTCCGGCGGTGATTAATACTTTCCGCACTGGATCCATGTTCGTCATCCTCCTGAACCTTg                                                        >  2:1820221/1‑93 (MQ=255)
                                     aTGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGa                                                     >  2:4476120/1‑93 (MQ=255)
                                     aTGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGa                                                     >  2:4271426/1‑93 (MQ=255)
                                      tGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGaa                                                    >  2:876613/1‑93 (MQ=255)
                                      tGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGaa                                                    >  2:515374/1‑93 (MQ=255)
                                      tGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGaa                                                    >  2:4188110/1‑93 (MQ=255)
                                       gTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAt                                                   <  2:61266/93‑1 (MQ=255)
                                         tGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcc                                                 <  1:1130053/93‑1 (MQ=255)
                                          gagaAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcct                                                >  2:1399761/1‑93 (MQ=255)
                                           agaAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtccta                                               <  2:2094749/93‑1 (MQ=255)
                                             aaGGTTGTTAGACCGCCACAAAATCCGGTGGTGCTTAATACTTTCCGCACTGGATCACTGTTCGTCATCCTGCTGAACTATGCGCAtcctatt                                             >  1:2496759/1‑93 (MQ=255)
                                             aaGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctatt                                             >  1:3429799/1‑93 (MQ=255)
                                              aGGCTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattc                                            <  2:2476290/93‑1 (MQ=255)
                                               ggTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattcc                                           >  2:3968639/1‑93 (MQ=255)
                                               ggTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattcc                                           <  1:3993177/93‑1 (MQ=255)
                                                 ttgttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattccta                                         <  1:389055/93‑1 (MQ=255)
                                                   tttAGACCCCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatg                                       <  2:1316304/92‑1 (MQ=255)
                                                   gttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatg                                       <  2:3146047/93‑1 (MQ=255)
                                                    ttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatga                                      >  1:4404776/1‑93 (MQ=255)
                                                    ttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatga                                      >  1:2582937/1‑93 (MQ=255)
                                                      aacccGCCAAAAAATCCGGTTGTGATTAATACTTTCCGCACTGGATCAGTGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatgatg                                    <  1:2309996/90‑1 (MQ=255)
                                                      aGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatgatg                                    <  1:2004110/93‑1 (MQ=255)
                                                      aGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACAGGATAAATGTTCGTCATCCTGCTGATCCATGCGAATCCTATTCCTatgatg                                    >  1:3074404/1‑93 (MQ=255)
                                                       gACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTAtgatga                                   >  1:3102071/1‑93 (MQ=255)
                                                       gACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTAtgatga                                   >  1:4304574/1‑93 (MQ=255)
                                                          cGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATg                                <  1:1906829/93‑1 (MQ=255)
                                                           cgcacaAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGc                               <  1:909318/91‑1 (MQ=255)
                                                            ccACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGcc                              <  2:2863120/93‑1 (MQ=255)
                                                            ccACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGcc                              <  2:1101179/93‑1 (MQ=255)
                                                            ccACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGcc                              <  2:4369507/93‑1 (MQ=255)
                                                             cacaAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGccc                             <  2:132445/93‑1 (MQ=255)
                                                             cacaAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGccc                             <  2:388649/93‑1 (MQ=255)
                                                              acaAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGcccc                            >  2:3762844/1‑93 (MQ=255)
                                                                aaaaTCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCaa                          >  1:2232571/1‑93 (MQ=255)
                                                                aaaaTCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCaa                          >  1:2814250/1‑93 (MQ=255)
                                                                 aaaTCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAAt                         >  2:2954868/1‑93 (MQ=255)
                                                                  aaTCCGGTGTTGATTAATACTTGCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATTCGAATCCTATTCCTATGATGAATGCCCCAATc                        <  2:3273108/93‑1 (MQ=255)
                                                                   aTCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCa                       <  2:3551709/93‑1 (MQ=255)
                                                                   aTCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGCTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCa                       >  1:1364613/1‑93 (MQ=255)
                                                                    tCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAg                      <  1:4629927/93‑1 (MQ=255)
                                                                    tCCGGAGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATACTGCTGAACCATGCGAATCCTATTCCTATGATTAATGCCCCAATCAg                      >  2:3019082/1‑93 (MQ=255)
                                                                     ccGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAgg                     >  2:3092176/1‑93 (MQ=255)
                                                                      cGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGt                    >  2:1489979/1‑93 (MQ=255)
                                                                      cGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGt                    >  1:959614/1‑93 (MQ=255)
                                                                      cGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGt                    <  1:71668/93‑1 (MQ=255)
                                                                       ggtggtGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGtt                   <  2:1642842/93‑1 (MQ=255)
                                                                         tggtgATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTg                 >  1:3283725/1‑93 (MQ=255)
                                                                           gtgATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCt               >  2:3325830/1‑93 (MQ=255)
                                                                           gtgATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCt               <  1:2998216/93‑1 (MQ=255)
                                                                           gtgATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCt               <  2:1048392/93‑1 (MQ=255)
                                                                              aTTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTc            >  2:2348710/1‑93 (MQ=255)
                                                                               ttAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCa           <  1:621804/93‑1 (MQ=255)
                                                                                 aaTACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGc         <  2:3934610/93‑1 (MQ=255)
                                                                                 aaTACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGc         >  2:1714813/1‑93 (MQ=255)
                                                                                  aTACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCg        >  1:2436064/1‑93 (MQ=255)
                                                                                   tACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGt       >  1:10144/1‑93 (MQ=255)
                                                                                   tACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGt       >  2:2564838/1‑93 (MQ=255)
                                                                                    aCTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTc      <  1:3663641/93‑1 (MQ=255)
                                                                                    aCTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTc      >  1:78401/1‑93 (MQ=255)
                                                                                     cTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTcc     >  2:159347/1‑93 (MQ=255)
                                                                                     cTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTcc     >  1:500457/1‑93 (MQ=255)
                                                                                      tttCCGCACTGGATCACTGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTccc    >  2:4600013/1‑93 (MQ=255)
                                                                                      tttCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTccc    >  1:406260/1‑93 (MQ=255)
                                                                                           |                                                                                         
CGGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCA  >  NC_000913/656878‑657058

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: