Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,559,003 G→A intergenic (+85/‑385) b2443 → / → b2444 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,559,0030GA99.3% 531.4 / ‑5.7 148intergenic (+85/‑385)b2443/b2444hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (74/73);  minor base T (0/1);  total (74/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 150 total aligned reads displayed.

AAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGG  >  NC_000913/2558913‑2559092
                                                                                          |                                                                                         
aaaTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATc                                                                                         >  2:4380257/1‑93 (MQ=255)
aaaTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTACATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATc                                                                                         >  1:2410765/1‑93 (MQ=255)
  aTAACGGACCTAGGCGGCCATTGTGCCCGCCTTGTCTTTGCCTTAATTCCATGAATGGAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCaa                                                                                       <  2:3239033/93‑1 (MQ=255)
   tAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCaaa                                                                                      >  1:2197336/1‑93 (MQ=255)
     aCGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAg                                                                                    >  1:4360930/1‑93 (MQ=255)
      cGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGt                                                                                   <  1:4019715/93‑1 (MQ=255)
        gAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTg                                                                                 <  1:2830624/93‑1 (MQ=255)
         aaTCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGc                                                                                <  1:1097491/93‑1 (MQ=255)
         aaTCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGc                                                                                <  1:2024984/93‑1 (MQ=255)
         aaTCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGc                                                                                >  2:2101854/1‑93 (MQ=255)
           tCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCg                                                                              <  1:3503633/93‑1 (MQ=255)
            cAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGt                                                                             <  1:4060472/93‑1 (MQ=255)
             aGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGtt                                                                            >  1:2579795/1‑93 (MQ=255)
              ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTa                                                                           >  2:183093/1‑93 (MQ=255)
              ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTa                                                                           >  2:2520935/1‑93 (MQ=255)
               gCGGGCAATGTGCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTaa                                                                          <  2:1129959/93‑1 (MQ=255)
                cGGTCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAg                                                                         <  1:3857780/93‑1 (MQ=255)
                 gggCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAgg                                                                        >  2:1444158/1‑93 (MQ=255)
                  ggCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGt                                                                       >  1:24387/1‑93 (MQ=255)
                   gCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTg                                                                      <  2:4029552/93‑1 (MQ=255)
                   gCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTg                                                                      >  2:864440/1‑93 (MQ=255)
                    cAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTgg                                                                     <  2:3001476/93‑1 (MQ=255)
                    cAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTgg                                                                     <  2:1171658/93‑1 (MQ=255)
                    cAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTgg                                                                     <  1:1809853/93‑1 (MQ=255)
                    cAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTgg                                                                     <  2:766909/93‑1 (MQ=255)
                       tGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTTAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTc                                                                  >  2:3732124/1‑93 (MQ=255)
                          tCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAANNTTGCCGTTAAGGTGGGTCnaa                                                               >  1:2232935/1‑93 (MQ=255)
                          tCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCaaa                                                               >  1:2490704/1‑93 (MQ=255)
                           cccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                              <  1:1663512/93‑1 (MQ=255)
                            ccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGTTCAAAAGTTGCCGTTAAGGTGGGTCAAAAg                                                             <  2:1343285/93‑1 (MQ=255)
                              gCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                           <  2:40669/93‑1 (MQ=255)
                               ccTTTTCTTTGCCTTAATTCCATGNATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                          <  2:704147/93‑1 (MQ=255)
                               ccTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                          <  1:2450825/93‑1 (MQ=255)
                               ccTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                          <  1:3958121/93‑1 (MQ=255)
                                cTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                                         >  2:1544422/1‑93 (MQ=255)
                                cTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                                         <  2:2907209/93‑1 (MQ=255)
                                cTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                                         >  2:2708949/1‑93 (MQ=255)
                                cTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                                         <  1:1200847/93‑1 (MQ=255)
                                 ttttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                                        >  1:4470709/1‑93 (MQ=255)
                                  tttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCg                                                       >  2:1625766/1‑93 (MQ=255)
                                  tttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCg                                                       >  1:118210/1‑93 (MQ=255)
                                  tttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGCTAGTGGTGGATCAAAAGTTGCCGTTAAGGAGGGTCAAAAGTTGCCg                                                       >  2:334693/1‑93 (MQ=255)
                                    tCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGtt                                                     >  2:2813491/1‑93 (MQ=255)
                                      tttGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTaa                                                   >  2:3255436/1‑93 (MQ=255)
                                      tttGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTaa                                                   <  2:4249776/93‑1 (MQ=255)
                                       ttGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAg                                                  <  2:3583405/93‑1 (MQ=255)
                                         gCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                >  2:2637062/1‑93 (MQ=255)
                                         gCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                <  2:2328463/93‑1 (MQ=255)
                                         gCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                <  1:596340/93‑1 (MQ=255)
                                         gCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                <  2:3358789/93‑1 (MQ=255)
                                          ccTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTg                                               <  1:1148034/93‑1 (MQ=255)
                                            ttAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                             >  1:3167602/1‑93 (MQ=255)
                                              aaTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                           <  2:446498/93‑1 (MQ=255)
                                               aTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                          >  1:3481540/1‑93 (MQ=255)
                                                  ccATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                       >  2:2100910/1‑93 (MQ=255)
                                                  ccATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                       >  1:1082139/1‑93 (MQ=255)
                                                  ccATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                       <  1:3958120/93‑1 (MQ=255)
                                                  ccATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                       >  1:4128614/1‑93 (MQ=255)
                                                   cATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAg                                      >  1:3264146/1‑93 (MQ=255)
                                                   cATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAg                                      >  1:3594388/1‑93 (MQ=255)
                                                     tGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                    <  2:1481357/93‑1 (MQ=255)
                                                     tGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                    <  2:460871/93‑1 (MQ=255)
                                                     tGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                    >  2:826215/1‑93 (MQ=255)
                                                     tGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                    >  2:3294405/1‑93 (MQ=255)
                                                       aaTTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                  >  1:609705/1‑93 (MQ=255)
                                                       aaTTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                  <  2:4262934/93‑1 (MQ=255)
                                                          tCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGa                               >  2:3406730/1‑93 (MQ=255)
                                                            aGGTGGGTCAAAAGTTGCCGTGAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGt                             <  2:3858301/93‑1 (MQ=255)
                                                             ggTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTa                            >  1:536192/1‑93 (MQ=255)
                                                              gTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTaa                           <  1:833170/93‑1 (MQ=255)
                                                              gTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTaa                           <  2:3128584/93‑1 (MQ=255)
                                                              gTGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGCaa                           <  1:4445434/93‑1 (MQ=255)
                                                               tGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAg                          >  1:4547097/1‑93 (MQ=255)
                                                               tGGGTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAg                          <  2:2941775/93‑1 (MQ=255)
                                                                gggTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGa                         >  2:3939004/1‑93 (MQ=255)
                                                                 ggTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGaa                        <  2:2423578/93‑1 (MQ=255)
                                                                  gTCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAAt                       <  2:2538233/93‑1 (MQ=255)
                                                                   tCAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAAtt                      >  1:4466188/1‑93 (MQ=255)
                                                                    cAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTc                     >  1:3385321/1‑93 (MQ=255)
                                                                    cAAAAGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTc                     <  1:2325354/93‑1 (MQ=255)
                                                                     aaaaGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTct                    >  1:617589/1‑93 (MQ=255)
                                                                      aaaGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta                   >  1:4546802/1‑93 (MQ=255)
                                                                        aGTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctact                 >  2:4091246/1‑93 (MQ=255)
                                                                        aGTTGACGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctact                 >  1:4108108/1‑93 (MQ=255)
                                                                         gTTGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctacta                <  2:3605261/93‑1 (MQ=255)
                                                                          ttGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAt               <  2:976846/93‑1 (MQ=255)
                                                                          ttGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAt               >  2:1993383/1‑93 (MQ=255)
                                                                          ttGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAt               <  1:3432868/93‑1 (MQ=255)
                                                                          ttGCCGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGCTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAt               >  1:758765/1‑93 (MQ=255)
                                                                             ccGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTaa            <  1:165613/93‑1 (MQ=255)
                                                                              cGTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAg           <  2:1128971/93‑1 (MQ=255)
                                                                               gTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGa          >  2:601194/1‑93 (MQ=255)
                                                                               gTTAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGCGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGa          >  1:3802996/1‑93 (MQ=255)
                                                                                ttAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGaa         >  2:276587/1‑93 (MQ=255)
                                                                                 tAGTGGTGGATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAAt        <  1:1740613/93‑1 (MQ=255)
                                                                                     ggtggATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa    <  1:3379111/93‑1 (MQ=255)
                                                                                     ggtggATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa    <  1:2034995/93‑1 (MQ=255)
                                                                                     ggtggATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa    <  1:1285556/93‑1 (MQ=255)
                                                                                      gtggATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAg   <  2:4595150/93‑1 (MQ=255)
                                                                                       tggATCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg  <  1:2303374/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGG  >  NC_000913/2558913‑2559092

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