Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,342,604 G→A A10A (GCC→GCT b4523 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,342,6040GA99.6% 819.4 / ‑5.2 231A10A (GCC→GCTb4523hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (145/85);  total (145/86)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 239 total aligned reads displayed.

CAAGAGTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTAAAAG  >  NC_000913/1342514‑1342692
                                                                                          |                                                                                        
 aaGAGTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCa                                                                                       >  1:6574935/1‑93 (MQ=255)
  agagTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCaa                                                                                      <  1:56530/93‑1 (MQ=255)
  agagTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCaa                                                                                      >  1:4059700/1‑93 (MQ=255)
   gagTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAAc                                                                                     >  2:611940/1‑93 (MQ=255)
    agTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAAcg                                                                                    >  2:8445087/1‑93 (MQ=255)
       tCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACgcgc                                                                                 <  1:8823962/93‑1 (MQ=255)
       tCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACgcgc                                                                                 >  1:7804582/1‑93 (MQ=255)
         tGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTg                                                                               >  2:3884324/1‑93 (MQ=255)
          gCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTgg                                                                              <  1:6155593/93‑1 (MQ=255)
          gCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATAGCTAACACAATATCAATACGTTCAGCAACCCGCTgg                                                                              >  1:8577099/1‑93 (MQ=255)
          gCTTGAAGTATCTCGATATTATGAATCGCAGAGTTGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTgg                                                                              <  1:266407/93‑1 (MQ=255)
            tttAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGc                                                                            >  2:8800175/1‑93 (MQ=255)
             ttAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCa                                                                           <  2:8062245/93‑1 (MQ=255)
               aaGTATCTCGAGATTATGGATCGCAGAGTGGTAAACGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATc                                                                         >  2:5712466/1‑93 (MQ=255)
                aGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCa                                                                        >  1:1183593/1‑93 (MQ=255)
                 gTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCaa                                                                       <  2:3525882/93‑1 (MQ=255)
                    tctcGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAAtt                                                                    >  1:7536643/1‑93 (MQ=255)
                    tctcGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAAtt                                                                    >  2:636900/1‑93 (MQ=255)
                    aatggagaTTATGGATCGAAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAAtt                                                                    <  2:3136264/89‑1 (MQ=255)
                     ctcGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTc                                                                   >  2:2387061/1‑93 (MQ=255)
                      tcGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATAAAATTCa                                                                  >  1:287962/1‑93 (MQ=255)
                      tcGAGATTATGGATCGCAGAGTGGGAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCa                                                                  >  2:6723650/1‑93 (MQ=255)
                      tcGAGATAATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATGCGTTCAGCAACGCGCTGGGCATCAAATTCa                                                                  <  2:7161646/93‑1 (MQ=255)
                       cGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAg                                                                 >  2:3167129/1‑93 (MQ=255)
                        gagaTTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGa                                                                >  1:7680525/1‑93 (MQ=255)
                        gagaTTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGa                                                                >  1:3661123/1‑93 (MQ=255)
                        gagaTTATGGATCGCAGAGGGGTAATCGCCTGCCACAAGGATACCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCACCAAATCCAGa                                                                >  1:2216139/1‑93 (MQ=255)
                         agaTTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGAc                                                               >  1:886341/1‑93 (MQ=255)
                           aTTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACAt                                                             <  2:8719970/93‑1 (MQ=255)
                            ttATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATc                                                            >  1:3690238/1‑93 (MQ=255)
                            ttATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATc                                                            <  1:6934947/93‑1 (MQ=255)
                             tATGGATCGCAGAGTGGTAATCGCCTTCCACAAGGATATCTAACACAATATCAATACGTTCAACAACGCGCTGGGCATCAAATTCAGACATCa                                                           >  1:5987764/1‑93 (MQ=255)
                             tATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCa                                                           >  2:819882/1‑93 (MQ=255)
                              aTGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATCTCTAACACAATATAAATACGATCAGCAACGAGCTGTACAACAAACTCAGAACTCaa                                                          >  2:2843282/1‑93 (MQ=255)
                                  accGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGAc                                                      <  1:2045183/91‑1 (MQ=255)
                                      cAGAGTGCTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATcc                                                  <  2:8308980/93‑1 (MQ=255)
                                         agTGGTACCCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTta                                               <  2:4991169/93‑1 (MQ=255)
                                         agTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTta                                               <  2:8538889/93‑1 (MQ=255)
                                         agTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTta                                               >  1:2343939/1‑93 (MQ=255)
                                           tGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTtata                                             <  1:1885757/93‑1 (MQ=255)
                                           tGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTtata                                             >  1:4289494/1‑93 (MQ=255)
                                            ggTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAg                                            <  1:3335680/93‑1 (MQ=255)
                                             gTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGc                                           <  2:7021098/93‑1 (MQ=255)
                                             gTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGc                                           <  1:4473582/93‑1 (MQ=255)
                                             gTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGc                                           <  2:7922787/93‑1 (MQ=255)
                                             gTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGc                                           <  1:3077286/93‑1 (MQ=255)
                                               aaTCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCaa                                         >  2:4261658/1‑93 (MQ=255)
                                               aaTCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCaa                                         >  2:1312506/1‑93 (MQ=255)
                                                aTCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCaaa                                        >  2:3409133/1‑93 (MQ=255)
                                                 tCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCaaaa                                       <  2:5186325/93‑1 (MQ=255)
                                                  cGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAg                                      <  2:4021043/93‑1 (MQ=255)
                                                   gCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAgg                                     >  1:3793887/1‑93 (MQ=255)
                                                   gCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAgg                                     <  2:5593957/93‑1 (MQ=255)
                                                    ccTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAggg                                    <  1:5213310/93‑1 (MQ=255)
                                                         cacaAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGaaa                               <  1:1122579/93‑1 (MQ=255)
                                                         cacaAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGaaa                               >  2:3933597/1‑93 (MQ=255)
                                                         cacaAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGaaa                               <  1:8413543/93‑1 (MQ=255)
                                                         cacaAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGaaa                               <  2:2646303/93‑1 (MQ=255)
                                                          acaAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGaaaa                              >  2:1241310/1‑93 (MQ=255)
                                                          acaAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGaaaa                              >  1:7867623/1‑93 (MQ=255)
                                                          acaAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGaaaa                              >  1:5396947/1‑93 (MQ=255)
                                                          acaAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGaaaa                              >  2:8646536/1‑93 (MQ=255)
                                                          acaAGGATATCTAACACAATATCAATACGTTCAGCAACCCGCTGGGCATCAAATTCAGACATCAACGACATCCTTATAGCAAAAGGGGGcaaa                              >  1:6403297/1‑93 (MQ=255)
                                                             aGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGa                           >  2:6193036/1‑93 (MQ=255)
                                                             aGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGa                           >  1:4849035/1‑93 (MQ=255)
                                                             aGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGa                           >  2:5244780/1‑93 (MQ=255)
                                                               gATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaat                         <  1:5527724/93‑1 (MQ=255)
                                                               gATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaat                         >  1:3194181/1‑93 (MQ=255)
                                                               gATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaat                         >  1:4712693/1‑93 (MQ=255)
                                                                atatCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaata                        <  1:1999882/93‑1 (MQ=255)
                                                                atatCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaata                        <  2:6490388/93‑1 (MQ=255)
                                                                 tatCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaataa                       >  2:2008983/1‑93 (MQ=255)
                                                                   tCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATg                     <  1:6697639/93‑1 (MQ=255)
                                                                    cTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGc                    >  1:1065102/1‑93 (MQ=255)
                                                                    cTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGc                    <  2:3378956/93‑1 (MQ=255)
                                                                    cTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGc                    <  2:4514776/93‑1 (MQ=255)
                                                                    cTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGc                    <  2:7263586/93‑1 (MQ=255)
                                                                     tAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCa                   >  1:6436462/1‑93 (MQ=255)
                                                                       acacAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGaaaa                 >  2:6784599/1‑93 (MQ=255)
                                                                       acacAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCaaa                 >  2:2261946/1‑93 (MQ=255)
                                                                        cacaATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCACAAGGGGGAAAATGAATAATGCaaaa                >  2:7495273/1‑93 (MQ=255)
                                                                        cacaATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCaaaa                <  2:7598070/93‑1 (MQ=255)
                                                                         acaATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCaaaaa               >  1:1063250/1‑93 (MQ=255)
                                                                            atatCAATACGTTGAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTg            >  1:5801589/1‑93 (MQ=255)
                                                                            atatCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTg            >  1:2167686/1‑93 (MQ=255)
                                                                              atCAATACGTTCAGCAACGCGCTGTGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGcc          <  2:6174075/93‑1 (MQ=255)
                                                                              atCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGcc          >  2:6884874/1‑93 (MQ=255)
                                                                               tCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCg         >  1:2472941/1‑93 (MQ=255)
                                                                                cAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGc        >  1:157174/1‑93 (MQ=255)
                                                                                 aaTACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCt       <  2:2152848/93‑1 (MQ=255)
                                                                                  aTACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTa      >  1:2984850/1‑93 (MQ=255)
                                                                                  aTACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTa      >  2:6314511/1‑93 (MQ=255)
                                                                                   tACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTaa     >  1:8598991/1‑93 (MQ=255)
                                                                                   tACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTaa     >  1:7707303/1‑93 (MQ=255)
                                                                                    aCGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTaaa    >  2:2672320/1‑93 (MQ=255)
                                                                                     cGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTaaaa   <  2:3622198/93‑1 (MQ=255)
                                                                                     cGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTaaaa   <  2:3056107/93‑1 (MQ=255)
                                                                                     cGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTaaaa   >  1:528175/1‑93 (MQ=255)
                                                                                     cGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTaaaa   <  2:2528899/93‑1 (MQ=255)
                                                                                      gTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTAAAag  >  2:1414180/1‑93 (MQ=255)
                                                                                          |                                                                                        
CAAGAGTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTAAAAG  >  NC_000913/1342514‑1342692

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