Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,281,967 T→C 53.7% Q370R (CAG→CGG)  yjcF ← pentapeptide repeats protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,281,9670TC53.7% 37.3 / 281.1 222Q370R (CAG→CGG) yjcFpentapeptide repeats protein
Reads supporting (aligned to +/- strand):  ref base T (53/49);  new base C (65/53);  total (119/103)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-02

Only 100 of 227 total aligned reads displayed.

GTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGCTTAAATTGG  >  NC_000913/4281868‑4282066
                                                                                                   |                                                                                                   
  ggCACAATCTTGTTAAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTACTTAGCACGTAGAATGGTAATAACGATTCCCt                                                                                                     >  2:3465728/1‑98 (MQ=255)
  ggCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTgg                                                                                                   <  1:1335723/100‑1 (MQ=255)
     acaATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACg                                                                                                >  2:1490274/1‑100 (MQ=255)
       aaTCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCt                                                                                              >  1:2586220/1‑100 (MQ=255)
         tcttctTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAg                                                                                            <  1:750645/100‑1 (MQ=255)
         tcttctTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAg                                                                                            >  1:926502/1‑100 (MQ=255)
          cttcttCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGc                                                                                           >  1:1205831/1‑100 (MQ=255)
          cttcttCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGc                                                                                           <  1:1436009/100‑1 (MQ=255)
           ttcttcAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATCTTTTCCTCACAACGGAGAATGGTAATACCCATTCCCCGGACGCTAGcc                                                                                          <  2:2208892/100‑1 (MQ=255)
            tcttcAAAGAATTTCTGTAGTGTTACCAGGCTATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCt                                                                                         >  2:766145/1‑100 (MQ=255)
              ttcAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTcc                                                                                       <  2:2680228/100‑2 (MQ=255)
              ttcAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCa                                                                                       >  1:3414449/1‑100 (MQ=255)
              ttcAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCa                                                                                       >  2:3071130/1‑100 (MQ=255)
               tcAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCAt                                                                                      >  1:2348593/1‑100 (MQ=255)
                cAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATc                                                                                     <  1:4133733/100‑1 (MQ=255)
                 aaaGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCg                                                                                    >  1:1132559/1‑100 (MQ=255)
                 aaaGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCg                                                                                    >  1:525087/1‑100 (MQ=255)
                    gAATTTCTGTAGTGTTACCAGGATATCACTATTAATTGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGaaa                                                                                 <  2:1697767/100‑1 (MQ=255)
                    gAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGTAGAATGGTAATACCCATTCTCCGGACGCTAGCCTCATCGaa                                                                                  >  2:1393788/1‑99 (MQ=255)
                      aTTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAAtt                                                                               <  1:1590339/100‑1 (MQ=255)
                      aTTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAAtt                                                                               <  1:4101638/100‑1 (MQ=255)
                          cTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAaca                                                                           <  1:79406/100‑1 (MQ=255)
                           tGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAacac                                                                          <  1:1039787/100‑1 (MQ=255)
                           ngTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAacac                                                                          >  2:1335804/2‑100 (MQ=255)
                            gTAGTGTTACCAGGCTATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGTTAATACCCATTCCCTGTACGCTAGCCTCATCGAAATTAacac                                                                          >  2:3586012/1‑99 (MQ=255)
                            gTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACg                                                                         >  2:3824705/1‑100 (MQ=255)
                              aGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGAt                                                                       <  2:1748962/100‑1 (MQ=255)
                               gtgtTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATa                                                                      <  1:2260974/100‑1 (MQ=255)
                                tgtTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATag                                                                     <  2:1942349/100‑1 (MQ=255)
                                tcttccccggACATCATTACTAATGGAGATACTTTCCTCACCACGCCGACTAGTAATACCCATCCCCCGTACCCTACCCTCATAGAAATTAACACGAt                                                                       >  2:1829528/9‑98 (MQ=255)
                                   tACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAgaga                                                                  >  1:661705/1‑100 (MQ=255)
                                   tACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAgaga                                                                  >  1:2294992/1‑100 (MQ=255)
                                   tACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAgaga                                                                  >  2:2851825/1‑100 (MQ=255)
                                    aCCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGAACGCTAGCCTCATCGAAATTAACACGATAGAGAc                                                                 >  2:2179626/1‑100 (MQ=255)
                                     ccAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATATCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACg                                                                <  2:642260/100‑1 (MQ=255)
                                     ccAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCTCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACg                                                                <  2:3408438/100‑1 (MQ=255)
                                     ccAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACg                                                                >  1:573372/1‑100 (MQ=255)
                                     ccAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACg                                                                <  2:97455/100‑1 (MQ=255)
                                     ccAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACg                                                                >  1:1834532/1‑100 (MQ=255)
                                      cAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCACCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACgg                                                               >  2:3641809/1‑100 (MQ=255)
                                         gATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGa                                                              <  2:493623/98‑1 (MQ=255)
                                          atatCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGatat                                                           >  1:3870223/1‑100 (MQ=255)
                                              cACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTaaa                                                       >  1:3786859/1‑100 (MQ=255)
                                                cTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACa                                                     >  1:3503666/1‑100 (MQ=255)
                                                  aTTAATGTAGATATTTTCCTCACCACGGAGAATCGTAATACCCATTCCCCTGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATc                                                   <  2:4052542/100‑1 (MQ=255)
                                                      aTGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAgg                                               >  1:3719099/1‑100 (MQ=255)
                                                      aTGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAgg                                               <  1:3071149/100‑1 (MQ=255)
                                                       tGGAGATATTTGCCTCACCACGGAGAATGGTAATAACCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGc                                              <  2:3645415/100‑1 (MQ=255)
                                                        ggAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCt                                             >  1:3216138/1‑100 (MQ=255)
                                                             tatTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGt                                        <  1:2461200/100‑1 (MQ=255)
                                                             tatTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGt                                        >  1:3470712/1‑100 (MQ=255)
                                                             tatTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCAAGGCTTTTGt                                        >  2:37295/1‑100 (MQ=255)
                                                              atTTTCCTCACCACGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGtt                                       <  2:2694978/100‑1 (MQ=255)
                                                                    cTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAaaca                                 <  1:1468213/100‑1 (MQ=255)
                                                                     tCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAaacag                                >  1:148896/1‑100 (MQ=255)
                                                                     tCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAaacag                                >  1:2762004/1‑100 (MQ=255)
                                                                     tCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAaacag                                <  1:2120067/100‑1 (MQ=255)
                                                                     tCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAaacag                                >  1:1957282/1‑100 (MQ=255)
                                                                      caccacGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAaacaga                               <  2:1677301/100‑1 (MQ=255)
                                                                       accacGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAaacagaa                              <  1:390418/100‑1 (MQ=255)
                                                                       accacGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAaacagaa                              <  2:1451673/100‑1 (MQ=255)
                                                                       accacGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAaacagaa                              <  1:3429477/100‑1 (MQ=255)
                                                                        ccacGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAaacagaac                             <  2:3702741/100‑1 (MQ=255)
                                                                        ccacGGAGAATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAaacagaa                              <  2:885558/99‑1 (MQ=255)
                                                                          acttagaATGGTAATACCCAATCCCCAGAAGATAGCCTCATCGAAGTTAATACGATAGAGACGGAGATTAAACAGCCAGGCTTTTGTTAAaacagaacag                           <  2:2154753/96‑1 (MQ=255)
                                                                            gtgggctGGGAATAGCCATACTCCGGACGCGTTCCTCATCGTAAGTACCACGATAGCGACTAATATTAGACATCCAGGCTTTTGTTAAaacagaac                             <  2:3930641/90‑1 (MQ=255)
                                                                            ggAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACCGtt                         >  2:1150090/1‑100 (MQ=255)
                                                                            ggAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGtt                         >  2:1873988/1‑100 (MQ=255)
                                                                             gagaATGGTAATACCCATTGCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTg                        <  2:2102015/100‑1 (MQ=255)
                                                                             gagaATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTg                        >  1:819122/1‑100 (MQ=255)
                                                                             gagaATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATGGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTg                        <  2:2294979/100‑1 (MQ=255)
                                                                             gagaATGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTg                        <  2:776570/100‑1 (MQ=255)
                                                                               gaATGGTAATACCGATTCCCCGGACGCTAGCCTCATCAAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGaa                      >  1:2807134/1‑100 (MQ=255)
                                                                               gaATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGaa                      <  1:3364444/100‑1 (MQ=255)
                                                                                aaTGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGaaa                     <  2:4096819/100‑1 (MQ=255)
                                                                                aaTGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGaaa                     >  2:3119613/1‑100 (MQ=255)
                                                                                aaTGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGaaa                     <  1:1865208/100‑1 (MQ=255)
                                                                                  tGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATg                   >  1:3439833/1‑100 (MQ=255)
                                                                                  tGGTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATg                   <  1:2510733/100‑1 (MQ=255)
                                                                                    gTAATACCCATTCCCTGTACGCTAGCCTCATCGAAATTAACACGATAGATACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGc                 >  2:344714/1‑100 (MQ=255)
                                                                                    gTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGc                 >  1:2211400/1‑100 (MQ=255)
                                                                                    gTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGc                 >  2:614435/1‑100 (MQ=255)
                                                                                    gTAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATggn                 <  2:910274/100‑2 (MQ=255)
                                                                                     tAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCa                <  2:3503644/100‑1 (MQ=255)
                                                                                     tAATACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGc                 >  1:3586591/1‑99 (MQ=255)
                                                                                      aaTACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCaa               >  1:862578/1‑100 (MQ=255)
                                                                                      aaTACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCaa               >  1:3389408/1‑100 (MQ=255)
                                                                                      aaTACCCATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCaa               >  1:472499/1‑100 (MQ=255)
                                                                                          cccATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGc           >  2:3424799/1‑100 (MQ=255)
                                                                                          cccATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGc           <  2:1971672/100‑1 (MQ=255)
                                                                                          cccATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGc           >  1:1825571/1‑100 (MQ=255)
                                                                                          cccATTCCCTGGACGCTAGCCTCATCGAAATGAACACGATAGATACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGc           >  2:1636341/1‑100 (MQ=255)
                                                                                          cccATTCCCCGGACGCTAGCCTCATCGAAATTAACATGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGc           <  1:2905725/100‑1 (MQ=255)
                                                                                           ccATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGCt          >  2:2603692/1‑100 (MQ=255)
                                                                                            cATTCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGCtt         >  1:3171321/1‑100 (MQ=255)
                                                                                             aTTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTCTTAAAACAGAACAGTTGAAATGGCAATTGCTTa        <  2:2474048/100‑1 (MQ=255)
                                                                                              ttCCCCGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGCTTaa       <  2:3206460/100‑1 (MQ=255)
                                                                                               tCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGCTTaaa      <  2:2348580/100‑1 (MQ=255)
                                                                                                ccccGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGCTTAAAt     >  2:1088694/1‑100 (MQ=255)
                                                                                                  cTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGCTTAAATTg   <  1:3223597/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GTGGCACAATCTTCTTCAAAGAATTTCTGTAGTGTTACCAGGATATCACTATTAATGGAGATATTTTCCTCACCACGGAGAATGGTAATACCCATTCCCTGGACGCTAGCCTCATCGAAATTAACACGATAGAGACGGATATTAAACATCCAGGCTTTTGTTAAAACAGAACAGTTGAAATGGCAATTGCTTAAATTGG  >  NC_000913/4281868‑4282066

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: