Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,657,532 T→C H480H (CAT→CAC b1587 → putative oxidoreductase major subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,657,5320TC95.7% 334.6 / 4.3 116H480H (CAT→CACb1587putative oxidoreductase major subunit
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (54/57);  total (56/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

Only 100 of 118 total aligned reads displayed.

AATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGC  >  NC_000913/1657442‑1657621
                                                                                          |                                                                                         
aaTCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTc                                                                                         <  1:4091919/93‑1 (MQ=255)
  tCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTg                                                                                       >  1:3948542/1‑93 (MQ=255)
  tCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTg                                                                                       >  1:1838998/1‑93 (MQ=255)
  tCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTg                                                                                       <  1:567379/93‑1 (MQ=255)
   ccGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGa                                                                                      <  2:4092762/93‑1 (MQ=255)
   ccGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGa                                                                                      >  2:2982180/1‑93 (MQ=255)
    cgcgACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGAc                                                                                     <  2:3452263/93‑1 (MQ=255)
     gcgACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACa                                                                                    >  2:1845036/1‑93 (MQ=255)
        aCGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATca                                                                                 >  1:45260/1‑93 (MQ=255)
            cGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATcaacaa                                                                             <  2:3190898/93‑1 (MQ=255)
             gTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACaaa                                                                            <  2:2341065/93‑1 (MQ=255)
             gTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACaaa                                                                            >  2:3142588/1‑93 (MQ=255)
                  cgGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCa                                                                       <  1:346852/93‑1 (MQ=255)
                       aaGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGaaa                                                                  <  2:4531181/93‑1 (MQ=255)
                        aGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAAt                                                                 >  2:1005237/1‑93 (MQ=255)
                        aGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAAt                                                                 >  1:3345768/1‑93 (MQ=255)
                        aGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAAt                                                                 <  1:2296879/93‑1 (MQ=255)
                        aGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAAt                                                                 <  2:798491/93‑1 (MQ=255)
                         gATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAAtt                                                                >  1:96768/1‑93 (MQ=255)
                         gATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAAtt                                                                <  2:4596679/93‑1 (MQ=255)
                          aTAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTc                                                               >  2:1484776/1‑93 (MQ=255)
                          aTAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTc                                                               >  2:1885429/1‑93 (MQ=255)
                          aTAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTc                                                               >  2:3137548/1‑93 (MQ=255)
                           tAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATAGCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCt                                                              >  1:4325293/1‑93 (MQ=255)
                           tAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCt                                                              <  2:425713/93‑1 (MQ=255)
                           tAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCt                                                              >  1:1349828/1‑93 (MQ=255)
                             aaCTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGCACTTCTGc                                                            >  1:2275148/1‑93 (MQ=255)
                             aaCTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGc                                                            <  1:2104921/93‑1 (MQ=255)
                              aCTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCa                                                           <  2:1566902/93‑1 (MQ=255)
                                tGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAgg                                                         >  1:376054/1‑93 (MQ=255)
                                tGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAgg                                                         >  2:1144764/1‑93 (MQ=255)
                                tGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAgg                                                         >  1:2739522/1‑93 (MQ=255)
                                 ggATGTGCGCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGa                                                        <  1:2131110/93‑1 (MQ=255)
                                 ggATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGa                                                        <  1:3023790/93‑1 (MQ=255)
                                  gATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGAc                                                       >  1:2007263/1‑93 (MQ=255)
                                  gATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGAc                                                       >  2:526490/1‑93 (MQ=255)
                                   aTGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACg                                                      >  1:297126/1‑93 (MQ=255)
                                     gtgCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGaa                                                    <  2:3701779/93‑1 (MQ=255)
                                      tgCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAAt                                                   >  2:872584/1‑93 (MQ=255)
                                      tgCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAAt                                                   >  1:2790078/1‑93 (MQ=255)
                                       gCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATc                                                  <  1:4087849/93‑1 (MQ=255)
                                        cccATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcg                                                 <  2:3290937/93‑1 (MQ=255)
                                         ccATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcga                                                >  1:3748130/1‑93 (MQ=255)
                                          cATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaa                                               >  2:1531539/1‑93 (MQ=255)
                                          cATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAAAGCATGAAATTCTGCAGGACGAATcgaa                                               >  2:4549174/1‑93 (MQ=255)
                                           aTTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaa                                              >  2:4104049/1‑93 (MQ=255)
                                            ttAAGTTCATCTGGAACTACGCGGGAAATCCCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaat                                             <  2:4535490/93‑1 (MQ=255)
                                            ttAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaat                                             <  2:4213675/93‑1 (MQ=255)
                                            ttAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaat                                             >  2:3655186/1‑93 (MQ=255)
                                            ttAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaat                                             >  1:1758454/1‑93 (MQ=255)
                                             tAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatg                                            <  2:431815/93‑1 (MQ=255)
                                              aaGTTCATCTGGACCTACGCGGCACATGCCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatgc                                           <  2:2554203/93‑1 (MQ=255)
                                              aaGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatgc                                           <  1:4359184/93‑1 (MQ=255)
                                              aaGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatgc                                           <  2:3410586/93‑1 (MQ=255)
                                              aaGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatgc                                           >  1:1922962/1‑93 (MQ=255)
                                                gTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgacatgcga                                         >  2:2116820/1‑93 (MQ=255)
                                                   cATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGTACGAATcgaaatgcgaaat                                      <  2:1498991/93‑1 (MQ=255)
                                                   cATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatgcgaaat                                      <  1:3742646/93‑1 (MQ=255)
                                                    aTCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATcgaaatgcgaaatg                                     <  1:2721450/93‑1 (MQ=255)
                                                      cTGGAACTACGCGGGAAATCCCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGAt                                   <  1:2735228/93‑1 (MQ=255)
                                                      cTGGAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGAt                                   <  1:635563/93‑1 (MQ=255)
                                                        ggAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATtg                                 <  1:1572794/93‑1 (MQ=255)
                                                        ggAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATtg                                 >  2:3416146/1‑93 (MQ=255)
                                                         gAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATtgt                                <  2:4481599/93‑1 (MQ=255)
                                                         gAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATtgt                                <  2:3888376/93‑1 (MQ=255)
                                                         gAACTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATtgt                                <  2:740544/93‑1 (MQ=255)
                                                           aCTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTgg                              <  2:4271610/93‑1 (MQ=255)
                                                            cTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTggc                             >  2:736692/1‑92 (MQ=255)
                                                            cTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGt                             <  1:3226493/93‑1 (MQ=255)
                                                            cTACGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGt                             >  1:2574920/1‑93 (MQ=255)
                                                              aCGCGGGAAATACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCa                           <  1:1907892/93‑1 (MQ=255)
                                                                     aaaTCCCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGaaaa                    <  1:995299/93‑1 (MQ=255)
                                                                      aaTACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGaaaat                   >  2:4388526/1‑92 (MQ=255)
                                                                      aaTACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAAc                   >  1:1593547/1‑93 (MQ=255)
                                                                        tACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACtt                 >  2:2098110/1‑93 (MQ=255)
                                                                        tACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACtt                 <  2:2192972/93‑1 (MQ=255)
                                                                        tACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACtt                 >  2:4631835/1‑93 (MQ=255)
                                                                        tACCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACtt                 >  1:3441714/1‑93 (MQ=255)
                                                                         aCCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACttt                >  2:600689/1‑93 (MQ=255)
                                                                         aCCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACttt                >  2:4665519/1‑93 (MQ=255)
                                                                         aCCCTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACttt                >  1:2106084/1‑93 (MQ=255)
                                                                          cccTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTa               <  1:716135/93‑1 (MQ=255)
                                                                          cccTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTa               <  1:1016360/93‑1 (MQ=255)
                                                                           ccTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTAt              <  1:4518885/93‑1 (MQ=255)
                                                                           ccTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTAt              <  1:2923581/93‑1 (MQ=255)
                                                                            cTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAACTGCGAAATGATTGTGGTCATCGAACACTTTATg             >  1:1260923/1‑93 (MQ=255)
                                                                            cTCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATg             <  2:1291338/93‑1 (MQ=255)
                                                                             tCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGa            >  1:638280/1‑93 (MQ=255)
                                                                             tCGTTAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGa            >  2:767342/1‑93 (MQ=255)
                                                                                ttAATCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCt         <  1:206398/93‑1 (MQ=255)
                                                                                  aaTCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCa       <  1:673389/93‑1 (MQ=255)
                                                                                  aaTCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCa       <  1:4289657/93‑1 (MQ=255)
                                                                                   aTCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCAt      <  2:1753136/93‑1 (MQ=255)
                                                                                   aTCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCAt      >  2:1274821/1‑93 (MQ=255)
                                                                                    tCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATc     >  2:2986969/1‑93 (MQ=255)
                                                                                    tCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATc     <  2:2682801/93‑1 (MQ=255)
                                                                                    tCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATc     <  1:2859953/93‑1 (MQ=255)
                                                                                    tCAGCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATc     <  1:2068348/93‑1 (MQ=255)
                                                                                       gCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGc  >  2:3438981/1‑93 (MQ=255)
                                                                                       gCACTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGc  >  2:513379/1‑93 (MQ=255)
                                                                                          |                                                                                         
AATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCCATTAAGTTCATCTGGAACTACGCGGGAAATACCCTCGTTAATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGC  >  NC_000913/1657442‑1657621

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