Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,571,055 T→C L724L (TTA→TTG pqqL ← putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,571,0550TC88.4% 785.8 / 83.6 303L724L (TTA→TTGpqqLputative peptidase
Reads supporting (aligned to +/- strand):  ref base T (17/18);  new base C (133/135);  total (150/153)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

Only 100 of 303 total aligned reads displayed.

TTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGG  >  NC_000913/1570965‑1571145
                                                                                          |                                                                                          
 tGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAg                                                                                         >  1:1780308/1‑93 (MQ=255)
  gTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGt                                                                                        >  2:1920244/1‑93 (MQ=255)
  gTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGt                                                                                        >  2:2832368/1‑93 (MQ=255)
  gTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGt                                                                                        <  1:2823058/93‑1 (MQ=255)
     cTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAAc                                                                                     >  1:6752537/1‑93 (MQ=255)
      ttttACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAAcg                                                                                    <  2:7154986/93‑1 (MQ=255)
      ttttACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAAcg                                                                                    >  1:2456388/1‑93 (MQ=255)
       tttACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAAcgc                                                                                   >  2:2283820/1‑93 (MQ=255)
        ttACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAAcgcg                                                                                  <  2:639392/93‑1 (MQ=255)
         tACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGt                                                                                 >  1:3803534/1‑93 (MQ=255)
          aCAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTa                                                                                >  2:6880649/1‑93 (MQ=255)
          aCAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTa                                                                                >  2:1477006/1‑93 (MQ=255)
           cAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTaa                                                                               <  2:5200942/93‑1 (MQ=255)
            aGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAAt                                                                              >  2:739040/1‑93 (MQ=255)
             gTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAAtt                                                                             >  1:4531527/1‑93 (MQ=255)
             gTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAAtt                                                                             >  1:2958114/1‑93 (MQ=255)
              tAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTa                                                                            <  1:5082455/93‑1 (MQ=255)
              tAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTa                                                                            <  2:3298055/93‑1 (MQ=255)
              tAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTa                                                                            >  2:3364666/1‑93 (MQ=255)
               aaCCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTaa                                                                           <  2:4290148/93‑1 (MQ=255)
                aCCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAAc                                                                          >  1:1951776/1‑93 (MQ=255)
                 ccGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACg                                                                         <  2:3912933/93‑1 (MQ=255)
                    aTGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCa                                                                      <  1:6881236/93‑1 (MQ=255)
                    aTGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCa                                                                      >  2:824899/1‑93 (MQ=255)
                      gCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACg                                                                    >  2:151271/1‑93 (MQ=255)
                       cGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGa                                                                   <  1:438447/93‑1 (MQ=255)
                        gTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAg                                                                  <  2:6827563/93‑1 (MQ=255)
                         ttGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGt                                                                 >  1:762924/1‑93 (MQ=255)
                         ttGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGt                                                                 <  1:4063044/93‑1 (MQ=255)
                           gTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGttt                                                               <  1:3854821/93‑1 (MQ=255)
                           gTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGttt                                                               <  2:5486947/93‑1 (MQ=255)
                           gTCCGACGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGttt                                                               >  1:4698691/1‑93 (MQ=255)
                              cGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGtc                                                            <  2:7505493/93‑1 (MQ=255)
                               gTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtct                                                           >  2:2014345/1‑93 (MQ=255)
                               gTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtct                                                           >  2:5656065/1‑93 (MQ=255)
                                tCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGtctt                                                          >  1:2079721/1‑93 (MQ=255)
                                tCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGtctt                                                          >  1:6411000/1‑93 (MQ=255)
                                tCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGGGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGtctt                                                          >  1:5337011/1‑93 (MQ=255)
                                 cgcgcgAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtcttc                                                         <  1:7185294/93‑1 (MQ=255)
                                 cgcgcgAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGtcttc                                                         <  1:1017248/93‑1 (MQ=255)
                                 cgcgcgAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGtcttc                                                         <  1:4438280/93‑1 (MQ=255)
                                   cgcgAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTg                                                       >  1:2206236/1‑93 (MQ=255)
                                    gcgAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGc                                                      <  1:4594285/93‑1 (MQ=255)
                                      gAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGa                                                    >  2:379865/1‑93 (MQ=255)
                                       aGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGAc                                                   <  1:2252623/93‑1 (MQ=255)
                                        gTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACa                                                  <  1:4865029/93‑1 (MQ=255)
                                        gTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACa                                                  >  2:6549498/1‑93 (MQ=255)
                                           aaTGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTa                                               >  2:1213174/1‑93 (MQ=255)
                                            aTGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAc                                              <  2:2908791/93‑1 (MQ=255)
                                            aTGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAc                                              <  2:1085149/93‑1 (MQ=255)
                                            aTGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAc                                              >  1:1065052/1‑93 (MQ=255)
                                             tGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTAcc                                             <  1:3992357/93‑1 (MQ=255)
                                               gTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGa                                           >  2:6747475/1‑93 (MQ=255)
                                                 ttACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACa                                         <  2:4607452/93‑1 (MQ=255)
                                                 ttACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACa                                         >  2:4129607/1‑93 (MQ=255)
                                                  tACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACaa                                        <  1:5826658/93‑1 (MQ=255)
                                                  tACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACaa                                        <  1:5448597/93‑1 (MQ=255)
                                                   aCCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAAt                                       <  1:1929484/93‑1 (MQ=255)
                                                    ccTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATg                                      >  1:2290388/1‑93 (MQ=255)
                                                       gCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACa                                   >  2:4847598/1‑93 (MQ=255)
                                                       gCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACa                                   >  1:3196043/1‑93 (MQ=255)
                                                         ggCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACaaa                                 <  1:4216444/93‑1 (MQ=255)
                                                         ggCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACaaa                                 >  1:5526327/1‑93 (MQ=255)
                                                         ggCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACaaa                                 <  2:4366481/93‑1 (MQ=255)
                                                          gCTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAAc                                >  2:2456449/1‑93 (MQ=255)
                                                           cTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACg                               <  1:6589949/93‑1 (MQ=255)
                                                           cTAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACg                               >  2:5734784/1‑93 (MQ=255)
                                                            tAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGt                              >  1:2747640/1‑93 (MQ=255)
                                                            tAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGt                              <  2:1799982/93‑1 (MQ=255)
                                                            tAATGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGt                              <  1:1248636/93‑1 (MQ=255)
                                                             aaTGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTg                             >  2:1768489/1‑93 (MQ=255)
                                                              aTGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGa                            <  2:6992381/93‑1 (MQ=255)
                                                              aTGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGa                            <  1:873666/93‑1 (MQ=255)
                                                               tGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGat                           >  1:6274527/1‑93 (MQ=255)
                                                               tGGCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGat                           <  2:7248638/93‑1 (MQ=255)
                                                                ggCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGata                          <  1:6325891/93‑1 (MQ=255)
                                                                ggCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGata                          <  2:2232500/93‑1 (MQ=255)
                                                                ggCGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGata                          <  2:1563623/93‑1 (MQ=255)
                                                                  cGAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATc                        <  1:5749006/93‑1 (MQ=255)
                                                                   gAATCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATcc                       <  2:6937156/93‑1 (MQ=255)
                                                                    aaTCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCg                      >  2:5150293/1‑93 (MQ=255)
                                                                      tCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCt                    <  1:3466908/93‑1 (MQ=255)
                                                                      tCAGAGTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCt                    >  1:1519844/1‑93 (MQ=255)
                                                                        agagTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTgg                  >  1:7428612/1‑93 (MQ=255)
                                                                        agagTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTgg                  >  1:6788127/1‑93 (MQ=255)
                                                                         gagTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGga                 >  2:619160/1‑93 (MQ=255)
                                                                         gagTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGga                 <  1:1426024/93‑1 (MQ=255)
                                                                         gagTGTTTGATTGATCCCAAGTAACGCATAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCTGCTGga                 >  2:5008007/1‑93 (MQ=255)
                                                                          agTGTTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGgag                <  2:5807689/93‑1 (MQ=255)
                                                                           gtgtTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGgaga               >  1:5192971/1‑93 (MQ=255)
                                                                            tgtTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGAt              <  2:6138703/93‑1 (MQ=255)
                                                                            tgtTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGAt              >  1:6149779/1‑93 (MQ=255)
                                                                            tgtTTGATTGATCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGAt              <  1:4311808/93‑1 (MQ=255)
                                                                                tgattgatCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGaaa          >  1:6012455/1‑93 (MQ=255)
                                                                                 gattgatCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCCGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGaaaa         >  2:2110796/1‑93 (MQ=255)
                                                                                   ttgatCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACa       >  2:3496805/1‑93 (MQ=255)
                                                                                     gatCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAAt     <  2:6670885/93‑1 (MQ=255)
                                                                                      atCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAAtt    >  2:5373024/1‑93 (MQ=255)
                                                                                      atCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAAtt    >  1:4188869/1‑93 (MQ=255)
                                                                                       tCCCAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTg   <  1:4806025/93‑1 (MQ=255)
                                                                                          |                                                                                          
TTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGG  >  NC_000913/1570965‑1571145

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