Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,038,386 T→C L475P (CTG→CCG)  appC → putative third cytochrome oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,038,3860TC85.8% 619.9 / 94.0 261L475P (CTG→CCG) appCputative third cytochrome oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (17/20);  new base C (129/94);  total (146/115)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-02

Only 100 of 264 total aligned reads displayed.

CGAGTTTGGTCGTCAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAGCGCGATGCAGA  >  NC_000913/1038294‑1038478
                                                                                            |                                                                                            
cGAGTTTGGTCGTCAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACt                                                                                              >  2:5398783/1‑93 (MQ=255)
  aGTTTGGTCGTCAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGa                                                                                            <  2:1290912/93‑1 (MQ=255)
  aGTTTGGTCGTCAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGa                                                                                            <  2:4329125/93‑1 (MQ=255)
   gTTTGGTCGTCAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGAt                                                                                           <  2:7042360/93‑1 (MQ=255)
   gTTTGGTCGTCAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGAt                                                                                           >  2:159327/1‑93 (MQ=255)
       ggTCGTCAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATg                                                                                       <  1:6207299/93‑1 (MQ=255)
         tcgtcAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGATCATGAt                                                                                     <  1:7160122/93‑1 (MQ=255)
         tcgtcAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGAt                                                                                     <  1:5476458/93‑1 (MQ=255)
          cgtcAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATc                                                                                    <  2:903101/93‑1 (MQ=255)
          cgtcAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATc                                                                                    <  1:5629365/93‑1 (MQ=255)
            tcAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGATCATGATCGt                                                                                  >  1:6863922/1‑93 (MQ=255)
            tcAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGt                                                                                  <  1:2180334/93‑1 (MQ=255)
            tcAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGt                                                                                  >  1:1853863/1‑93 (MQ=255)
             cAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTa                                                                                 <  2:191939/93‑1 (MQ=255)
             cAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTa                                                                                 <  2:2876111/93‑1 (MQ=255)
             cAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTa                                                                                 <  2:1144323/93‑1 (MQ=255)
              aGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGATCATGATCGTAg                                                                                <  2:5140942/93‑1 (MQ=255)
              aGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAg                                                                                <  2:6512428/93‑1 (MQ=255)
               gCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAgg                                                                               <  2:1412956/93‑1 (MQ=255)
               gCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAgg                                                                               <  2:5731067/93‑1 (MQ=255)
                ccGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAggg                                                                              <  1:2481892/93‑1 (MQ=255)
                ccGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAggg                                                                              <  1:5073262/93‑1 (MQ=255)
                 cGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGc                                                                             <  2:3170820/93‑1 (MQ=255)
                 cGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGc                                                                             <  1:7418858/93‑1 (MQ=255)
                 cGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGc                                                                             <  2:4978739/93‑1 (MQ=255)
                  gTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCt                                                                            >  2:7129915/1‑93 (MQ=255)
                   tGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGATCATGATCGTAGGGCtt                                                                           <  2:5152274/93‑1 (MQ=255)
                   tGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGATCATGATCGTAGGGCtt                                                                           <  2:5076246/93‑1 (MQ=255)
                    gggCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGATCATGATCGTAGGGCttt                                                                          <  1:4153990/93‑1 (MQ=255)
                       cGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTaca                                                                       >  1:1528414/1‑93 (MQ=255)
                       cGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTaca                                                                       <  2:2527095/93‑1 (MQ=255)
                        gATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTacac                                                                      <  2:6605742/93‑1 (MQ=255)
                         aTACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACAcc                                                                     >  2:4630702/1‑93 (MQ=255)
                         aTACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACAcc                                                                     <  1:159983/93‑1 (MQ=255)
                            cAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTg                                                                  <  1:516995/93‑1 (MQ=255)
                            cAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTg                                                                  >  2:6225642/1‑93 (MQ=255)
                            cAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTg                                                                  >  1:3822743/1‑93 (MQ=255)
                            cAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTg                                                                  <  1:4231501/93‑1 (MQ=255)
                             aGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGt                                                                 >  2:3361361/1‑93 (MQ=255)
                              ggACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGtt                                                                >  1:1259983/1‑93 (MQ=255)
                               gACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTc                                                               <  2:4332609/93‑1 (MQ=255)
                               gACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTc                                                               <  1:4988093/93‑1 (MQ=255)
                               gACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTc                                                               <  2:7416339/93‑1 (MQ=255)
                                aCATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCt                                                              <  1:7125606/93‑1 (MQ=255)
                                 cATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCtt                                                             >  2:5552868/1‑93 (MQ=255)
                                  aTCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTa                                                            >  1:666463/1‑93 (MQ=255)
                                   tCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGATCATGATCGTAGGGCTTTACACCCTGTTCTTaa                                                           <  1:6362830/93‑1 (MQ=255)
                                   tCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTaa                                                           >  1:1928169/1‑93 (MQ=255)
                                       aCCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGc                                                       >  1:6011356/1‑93 (MQ=255)
                                        ccGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGcc                                                      >  2:5618368/1‑93 (MQ=255)
                                        ccGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGcc                                                      >  2:5294063/1‑93 (MQ=255)
                                          gACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGa                                                    >  2:3626236/1‑93 (MQ=255)
                                           aCATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGaa                                                   <  2:1776825/93‑1 (MQ=255)
                                           aCATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGaa                                                   >  2:5220977/1‑93 (MQ=255)
                                            cATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAg                                                  <  2:4455105/93‑1 (MQ=255)
                                            cATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAg                                                  >  2:2128240/1‑93 (MQ=255)
                                             aTACTCCGCGCCCTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGt                                                 <  1:3635043/93‑1 (MQ=255)
                                              tACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTc                                                <  1:3250481/93‑1 (MQ=255)
                                              tACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTc                                                >  2:1852396/1‑93 (MQ=255)
                                               aCTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCt                                               >  1:1388508/1‑93 (MQ=255)
                                               aCTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCt                                               <  1:464177/93‑1 (MQ=255)
                                                cTCCGCGCACTCCGCTTTTACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTa                                              >  2:4822631/1‑93 (MQ=255)
                                                cTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTa                                              <  1:1160355/93‑1 (MQ=255)
                                                 tCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTAc                                             >  1:1497436/1‑93 (MQ=255)
                                                    gcgcACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTg                                          >  2:5491052/1‑93 (MQ=255)
                                                    gcgcACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTg                                          >  1:4051883/1‑93 (MQ=255)
                                                     cgcACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGa                                         <  1:211271/93‑1 (MQ=255)
                                                       cACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATg                                       >  2:7440102/1‑93 (MQ=255)
                                                          tCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAg                                    >  1:2230234/1‑93 (MQ=255)
                                                          tCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAg                                    >  2:1512630/1‑93 (MQ=255)
                                                          tCCGCTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAg                                    >  1:467344/1‑93 (MQ=255)
                                                              cTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAat                                <  2:6329704/93‑1 (MQ=255)
                                                              cTTTAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAat                                >  2:2343640/1‑93 (MQ=255)
                                                               tttAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAata                               <  2:3406756/93‑1 (MQ=255)
                                                                ttAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAatat                              <  2:5952264/93‑1 (MQ=255)
                                                                 tAACCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATg                             <  2:4103692/93‑1 (MQ=255)
                                                                  aaCCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGc                            <  1:1775556/93‑1 (MQ=255)
                                                                   aCCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGcc                           <  1:7303182/93‑1 (MQ=255)
                                                                   aCCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGcc                           >  1:6460747/1‑93 (MQ=255)
                                                                   aCCACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGcc                           <  1:929478/93‑1 (MQ=255)
                                                                    ccACAGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGccc                          >  2:3237994/1‑93 (MQ=255)
                                                                     cacaGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCg                         <  1:2302426/93‑1 (MQ=255)
                                                                       caGGACAACTGGCTTTCTCACTGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTc                       <  2:5991466/93‑1 (MQ=255)
                                                                        aGGACAACTGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCt                      >  2:6664018/1‑93 (MQ=255)
                                                                              aCTGGCTTTCTCACTGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGcc                >  1:4432804/1‑93 (MQ=255)
                                                                                tGGCTTTCTCACTGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGa              >  1:1074159/1‑93 (MQ=255)
                                                                                tGGCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGa              <  2:392789/93‑1 (MQ=255)
                                                                                 ggCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAg             >  1:3526056/1‑93 (MQ=255)
                                                                                 ggCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAg             <  1:1689596/93‑1 (MQ=255)
                                                                                  gCTTTCTCACTGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAgc            <  2:4526265/93‑1 (MQ=255)
                                                                                  gCTTTCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAgc            <  1:5830213/93‑1 (MQ=255)
                                                                                   cTTTCTCACAGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAgcg           <  2:2102807/93‑1 (MQ=255)
                                                                                     ttCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAGcgcg         >  1:3009402/1‑93 (MQ=255)
                                                                                     ttCTCACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAGcgcg         >  1:7522532/1‑93 (MQ=255)
                                                                                      tctcACTGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAGCGCGa        <  2:6756618/93‑1 (MQ=255)
                                                                                      tctcACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAGCGCGa        >  2:2998936/1‑93 (MQ=255)
                                                                                      tctcACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAGCGCGa        >  1:4763572/1‑93 (MQ=255)
                                                                                       ctcACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAGCGCGAt       >  1:3032314/1‑93 (MQ=255)
                                                                                         cACCGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAGCGCGATGc     <  2:7223249/93‑1 (MQ=255)
                                                                                            tGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAGCGCGATGCaga  <  2:3436011/93‑1 (MQ=255)
                                                                                            |                                                                                            
CGAGTTTGGTCGTCAGCCGTGGGCGATACAGGACATCTTACCGACATACTCCGCGCACTCCGCTTTAACCACAGGACAACTGGCTTTCTCACTGATCATGATCGTAGGGCTTTACACCCTGTTCTTAATCGCCGAAGTCTACCTGATGCAGAAATATGCCCGTCTGGGGCCGAGCGCGATGCAGA  >  NC_000913/1038294‑1038478

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