Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 882,856 C→A 18.8% R178L (CGT→CTT)  ybjG ← undecaprenyl pyrophosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913882,8560CA18.8% 359.4 / 70.8 150R178L (CGT→CTT) ybjGundecaprenyl pyrophosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (67/54);  new base A (14/14);  total (82/68)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 152 total aligned reads displayed.

AGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCGGCCA  >  NC_000913/882757‑882953
                                                                                                   |                                                                                                 
aGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAc                                                                                                   <  1:2865991/100‑1 (MQ=255)
aGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAc                                                                                                   >  1:947794/1‑100 (MQ=255)
   gAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGtt                                                                                                <  1:3508660/100‑1 (MQ=255)
   gAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGtn                                                                                                <  2:1495107/100‑2 (MQ=255)
     gCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGa                                                                                              >  1:2989912/1‑100 (MQ=255)
        gCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATag                                                                                           >  1:155912/1‑100 (MQ=255)
        gCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATag                                                                                           >  1:2731704/1‑100 (MQ=255)
        ccTCAATACTACACTTTTTAGCATAGATCAGTCACGCACCCAGCCTTTGCGGATAGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATag                                                                                           <  2:3344461/99‑1 (MQ=255)
           cAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGt                                                                                        >  1:698852/1‑100 (MQ=255)
           cAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGt                                                                                        >  1:2781148/1‑100 (MQ=255)
           cAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGt                                                                                        >  1:3444217/1‑100 (MQ=255)
              tactacACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGACTGCAGACGTTGATAGAGTTTa                                                                                     <  2:359002/100‑1 (MQ=255)
               actacACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTAt                                                                                    <  1:3989717/100‑1 (MQ=255)
                ctacACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATg                                                                                   >  1:2952980/1‑100 (MQ=255)
                 tacACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGc                                                                                  <  1:1522741/100‑1 (MQ=255)
                 tacACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGc                                                                                  >  2:607934/1‑100 (MQ=255)
                  acacTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGcc                                                                                 >  1:1510586/1‑100 (MQ=255)
                  acacTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGcc                                                                                 >  2:1805786/1‑100 (MQ=255)
                  acacTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGcc                                                                                 >  1:3135791/1‑100 (MQ=255)
                  acacTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGcc                                                                                 >  1:3103362/1‑100 (MQ=255)
                   cacTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGTTGGTCGGTAATGCAACACAGACGCGACACCACGATTGCAGACGTTGATATAGTTTATGccc                                                                                >  2:2714701/1‑100 (MQ=255)
                   cacTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGccc                                                                                >  1:2531213/1‑100 (MQ=255)
                     cTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCa                                                                              >  1:1045324/1‑100 (MQ=255)
                     cTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCa                                                                              >  2:3944935/1‑100 (MQ=255)
                     cTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCa                                                                              >  2:3435996/1‑100 (MQ=255)
                       ttttAGCAGAGATCAGTCACGCACCCAGCTTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATc                                                                            <  2:1001422/100‑1 (MQ=255)
                       ttttAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATc                                                                            >  1:1114689/1‑100 (MQ=255)
                       ttttAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATc                                                                            <  1:3510365/100‑1 (MQ=255)
                        tttAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCg                                                                           >  1:744978/1‑100 (MQ=255)
                            gCAGAGATCAGTCACGCACCCAGCCTTTGCGGGTCGGTAATGCAAAACAGACGAGATCCCACGATTGCAGAAGTTCATAGAGGTTATGCCCCCACGCtt                                                                        >  2:4031515/1‑99 (MQ=255)
                            gCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCTTCGCTTg                                                                       <  1:2193254/100‑1 (MQ=255)
                            gCAGAGATCAGTCACGCACCCAGCCTTTGCGGAGCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTg                                                                       >  1:2473784/1‑100 (MQ=255)
                             cAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGc                                                                      >  1:4145130/1‑100 (MQ=255)
                              agagATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGcc                                                                     <  2:1472933/100‑1 (MQ=255)
                              agagATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGcc                                                                     >  2:750359/1‑100 (MQ=255)
                              agagATCAGTCACGCACCCAGCCTTTGCGGAGCGGTAATGCAAAACAGACCCGCTACCACGATTGCAGACGTGCATAGAGTTTATGCCCCa                                                                              >  2:2487789/1‑91 (MQ=255)
                               gagaTCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCa                                                                    <  1:375992/100‑1 (MQ=255)
                               gagaTCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCa                                                                    <  1:425362/100‑1 (MQ=255)
                                 gaTCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGa                                                                  >  1:3809783/1‑100 (MQ=255)
                                   tCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATa                                                                >  2:1512082/1‑100 (MQ=255)
                                     aGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGcagatcgga                                                              >  1:1621482/1‑92 (MQ=255)
                                     aGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAAt                                                              <  1:1559186/100‑1 (MQ=255)
                                      gTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATc                                                             <  1:1107488/100‑1 (MQ=255)
                                      gTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATc                                                             >  1:3938589/1‑100 (MQ=255)
                                       tCACGCACCCAGTCTTGGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATActct                                                            <  2:3853829/100‑1 (MQ=255)
                                        cACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTg                                                           >  2:1893124/1‑100 (MQ=255)
                                         aCGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTgg                                                          <  1:3455078/100‑1 (MQ=255)
                                         aCGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTgg                                                          <  2:1579263/100‑1 (MQ=255)
                                          cGCACCCAGCCTTTGCGGATCGGTAATGCAAAACCGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCTGATAATCTggg                                                         <  2:2187616/100‑1 (MQ=255)
                                          cGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTggg                                                         >  2:3184963/1‑100 (MQ=255)
                                          cGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTggg                                                         >  1:2313639/1‑100 (MQ=255)
                                           gCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGc                                                        <  1:174757/100‑1 (MQ=255)
                                           gCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGc                                                        <  1:2553552/100‑1 (MQ=255)
                                           gCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGc                                                        <  2:1510577/100‑1 (MQ=255)
                                           gCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGc                                                        <  2:202986/100‑1 (MQ=255)
                                              cccAGCCTTTGCGGATCGGTAATGCAACACAGACGCGATACCACGATTGCAGACGTTGATAGAGAATATGCCCCATCGCTTGCCAGATAATCTGGGCAc                                                      >  2:4135278/1‑99 (MQ=255)
                                              cccAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTGGCCAGATAATCTGGGCa                                                       >  2:47587/1‑98 (MQ=255)
                                               ccAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTa                                                    >  2:2255720/1‑100 (MQ=255)
                                                cAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTaa                                                   >  1:1321885/1‑100 (MQ=255)
                                                cAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTaa                                                   <  2:3278700/100‑1 (MQ=255)
                                                  gCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAgg                                                 <  1:2260599/100‑1 (MQ=255)
                                                   ccTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATCGGGTTTATGCCCCATCGCTTGCCAGATATTCTGGGCACTAAGGc                                                <  2:1779314/100‑1 (MQ=255)
                                                   ccTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGc                                                >  2:419081/1‑100 (MQ=255)
                                                       tGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGcc                                            <  1:2610788/100‑1 (MQ=255)
                                                          ggATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAAt                                         <  1:277973/100‑1 (MQ=255)
                                                           gATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATc                                        >  2:1910046/1‑100 (MQ=255)
                                                            aTCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCa                                       <  1:1859405/100‑1 (MQ=255)
                                                            aTCGGTAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCa                                       <  1:2879674/100‑1 (MQ=255)
                                                               ggTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATAc                                    <  2:1152761/100‑1 (MQ=255)
                                                               ggTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATAc                                    >  2:3561174/1‑100 (MQ=255)
                                                                gTAATGCAAAACAGACGCGATACCACGATTGCATACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGGCAATCATAcc                                   >  2:2736139/1‑100 (MQ=255)
                                                                gTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATAcc                                   >  2:550916/1‑100 (MQ=255)
                                                                 tAATGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCt                                  <  1:3646671/100‑1 (MQ=255)
                                                                  aaTGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTg                                 >  1:3713910/1‑100 (MQ=255)
                                                                  aaTGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTg                                 >  1:1224589/1‑100 (MQ=255)
                                                                    tGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGcc                               >  2:583475/1‑100 (MQ=255)
                                                                    tGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGcc                               <  2:2473769/100‑1 (MQ=255)
                                                                    tGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGcc                               >  1:1618838/1‑100 (MQ=255)
                                                                    tGCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGcc                               <  2:4090967/100‑1 (MQ=255)
                                                                     gCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCa                              >  2:985470/1‑100 (MQ=255)
                                                                     gCAAAACAGACGCGATACCACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCa                              >  1:1359550/1‑100 (MQ=255)
                                                                      cAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAg                             <  2:3237399/100‑1 (MQ=255)
                                                                      cAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAg                             >  2:2436802/1‑100 (MQ=255)
                                                                       aaaaCAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGc                            <  2:2029932/100‑1 (MQ=255)
                                                                          aCAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAAt                         <  1:3086284/100‑1 (MQ=255)
                                                                             gACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATGTGGGCACGAAAGCAGCCAATCATACCTGCCCGCAATCCa                      >  2:1117570/1‑100 (MQ=255)
                                                                                gtgATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCg                   <  2:2534229/98‑1 (MQ=255)
                                                                                 cgATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGa                  >  1:1017248/1‑100 (MQ=255)
                                                                                      ccACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATCATCTGGGCACTAACGCAGCCAATCACCCCTGCCAGCAATCAACCAAGCata             >  2:3414103/1‑100 (MQ=255)
                                                                                      ccACGATTGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCata             >  1:481590/1‑100 (MQ=255)
                                                                                       cACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCatat            >  1:2055602/1‑100 (MQ=255)
                                                                                             tgaaGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCg      <  2:3519055/97‑1 (MQ=255)
                                                                                             tGCAGAAGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCg      <  1:2643140/100‑1 (MQ=255)
                                                                                              gCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCgg     >  2:756166/1‑100 (MQ=255)
                                                                                              gCAGACGTTGATAGAGTTTATGCCCCATCGCTTACCAGATAATCTGGGCACTAAGGCAGCCAATCATACCAGCCAGCAATCCACCGAGCATATCCAGCgg     >  2:2500407/1‑100 (MQ=255)
                                                                                              gCAGACGTTGATAGAGTTTATGCCACATCGCATGCCAGATAATGTGGGCACTAAGGAAGCCAATCATACTTGCCAGCAATCCACCGAGCATATCCAGc       >  2:1615106/1‑98 (MQ=255)
                                                                                                aGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCGGcc   <  2:2952963/100‑1 (MQ=255)
                                                                                                aGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCGGcc   <  1:2095827/100‑1 (MQ=255)
                                                                                                aGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCAAGCAATCCACCGAGCATCTCCAGCgg     >  1:1384656/1‑98 (MQ=255)
                                                                                                 gACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCGGCCa  >  1:3240871/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
AGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGATCGGTAATGCAAAACAGACGCGATACCACGATTGCAGACGTTGATAGAGTTTATGCCCCATCGCTTGCCAGATAATCTGGGCACTAAGGCAGCCAATCATACCTGCCAGCAATCCACCGAGCATATCCAGCGGCCA  >  NC_000913/882757‑882953

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