Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 40,966 T→C A322A (GCA→GCG caiT ← putative carnitine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091340,9660TC84.8% 586.4 / 95.1 236A322A (GCA→GCGcaiTputative carnitine transporter
Reads supporting (aligned to +/- strand):  ref base T (22/14);  new base C (104/96);  total (126/110)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

Only 100 of 247 total aligned reads displayed.

AGCACCATGCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATACATCAGCAACAT  >  NC_000913/40874‑41058
                                                                                            |                                                                                            
   aCCATGCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGcc                                                                                           <  1:2791844/93‑1 (MQ=255)
   aCCATGCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGcc                                                                                           >  2:6408257/1‑93 (MQ=255)
     cATGCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCa                                                                                         >  1:4228881/1‑93 (MQ=255)
      aTGCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAg                                                                                        >  2:2669626/1‑93 (MQ=255)
      aTGCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAg                                                                                        <  2:2944652/93‑1 (MQ=255)
      aTGCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAg                                                                                        >  2:6467810/1‑93 (MQ=255)
       tGCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGt                                                                                       >  2:4524961/1‑93 (MQ=255)
       tGCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGt                                                                                       >  1:4027307/1‑93 (MQ=255)
        gCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTa                                                                                      >  1:1382037/1‑93 (MQ=255)
        gCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTa                                                                                      <  1:5731539/93‑1 (MQ=255)
        gCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTa                                                                                      >  2:1476700/1‑93 (MQ=255)
         ccGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAg                                                                                     <  2:6628761/93‑1 (MQ=255)
         ccGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAg                                                                                     <  2:2221442/93‑1 (MQ=255)
          cGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAGa                                                                                    <  1:1634220/93‑1 (MQ=255)
          cGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGa                                                                                    <  1:5382371/93‑1 (MQ=255)
           gAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAGaa                                                                                   >  1:2164149/1‑93 (MQ=255)
            aaGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAGAac                                                                                  >  2:4321436/1‑93 (MQ=255)
             aGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCTTCTGGATAGCATAAATCACCCACCACGCCCAGTAGAaca                                                                                 <  2:4786965/93‑1 (MQ=255)
             aGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAaca                                                                                 >  1:1135117/1‑93 (MQ=255)
              gCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAacac                                                                                <  2:1860984/93‑1 (MQ=255)
              gCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAacac                                                                                >  1:3517225/1‑93 (MQ=255)
              gCACAGGTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAGAacac                                                                                <  1:5106967/93‑1 (MQ=255)
               cacaGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATGACCCACCACGCCCAGTAGGACACg                                                                               >  1:6158025/1‑93 (MQ=255)
               cacaGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACg                                                                               <  2:3542602/93‑1 (MQ=255)
                 caGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGt                                                                             >  2:385319/1‑93 (MQ=255)
                  aGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTc                                                                            >  1:4303949/1‑93 (MQ=255)
                  aGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTc                                                                            >  2:4678858/1‑93 (MQ=255)
                   gTTCACGCCCAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTcc                                                                           <  1:3442628/93‑1 (MQ=255)
                   gTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTcc                                                                           <  2:5869522/93‑1 (MQ=255)
                      cACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGc                                                                        <  2:3893056/93‑1 (MQ=255)
                       aCGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAGAACACGGTCCAGcc                                                                       <  1:4824138/93‑1 (MQ=255)
                        cGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGccc                                                                      <  1:411198/93‑1 (MQ=255)
                         gCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCt                                                                     >  2:5368467/1‑93 (MQ=255)
                         gCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCt                                                                     >  2:5087194/1‑93 (MQ=255)
                          cacaGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTg                                                                    >  1:3045472/1‑93 (MQ=255)
                          cacaGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTg                                                                    <  1:4938965/93‑1 (MQ=255)
                              gTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCggg                                                                >  2:4488969/1‑93 (MQ=255)
                              gTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGACCTGCggg                                                                >  1:326594/1‑93 (MQ=255)
                               tACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGa                                                               <  1:28495/93‑1 (MQ=255)
                                   aCCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAgcc                                                           >  2:2907691/1‑93 (MQ=255)
                                      aCGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAgccgcc                                                        >  2:3418888/1‑93 (MQ=255)
                                      aCGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAgccgcc                                                        >  2:505903/1‑93 (MQ=255)
                                       cGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCt                                                       >  2:5532089/1‑93 (MQ=255)
                                       cGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCt                                                       <  2:2431924/93‑1 (MQ=255)
                                       cGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCt                                                       <  1:644418/93‑1 (MQ=255)
                                        gggAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCtt                                                      <  2:5208582/93‑1 (MQ=255)
                                        gggAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCtt                                                      >  2:1030587/1‑93 (MQ=255)
                                        gggAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCtt                                                      <  2:379883/93‑1 (MQ=255)
                                        gggAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCtt                                                      >  1:3629445/1‑93 (MQ=255)
                                         ggAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGTCCAGTAGAACACGGTCCAGCCCTGCGGGAAgccgcactt                                                     >  2:5500194/1‑89 (MQ=255)
                                         ggAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCttt                                                     <  2:158247/93‑1 (MQ=255)
                                             aTGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCg                                                 <  1:5996779/93‑1 (MQ=255)
                                             aTGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCg                                                 >  1:1688774/1‑93 (MQ=255)
                                              tGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGa                                                <  1:1541869/93‑1 (MQ=255)
                                              tGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGa                                                <  2:4228881/93‑1 (MQ=255)
                                              tGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGa                                                <  1:377400/93‑1 (MQ=255)
                                                cGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATg                                              <  1:1806531/93‑1 (MQ=255)
                                                  ggCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATggg                                            <  1:5763113/93‑1 (MQ=255)
                                                    cGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGAt                                          <  1:502226/93‑1 (MQ=255)
                                                     gAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATc                                         <  2:304051/93‑1 (MQ=255)
                                                     gAGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATc                                         <  2:4176293/93‑1 (MQ=255)
                                                      aGGAAGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGAGAAGCCGCCTTTAGCGATGGGATCg                                        >  2:2233954/1‑93 (MQ=255)
                                                      aGGAAGATACTCATCTGGAGAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCg                                        <  2:2347125/93‑1 (MQ=255)
                                                         aaGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGta                                     >  1:4438371/1‑93 (MQ=255)
                                                          aGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGtat                                    >  1:5331208/1‑93 (MQ=255)
                                                          aGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGtat                                    >  2:3144850/1‑93 (MQ=255)
                                                          aGATACTCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGtat                                    >  1:1455723/1‑93 (MQ=255)
                                                              aCTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGaa                                >  2:632905/1‑93 (MQ=255)
                                                                tCATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGaaca                              >  2:4212498/1‑93 (MQ=255)
                                                                 cATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGaacaa                             >  2:1276879/1‑93 (MQ=255)
                                                                 cATCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGaacaa                             <  2:6126255/93‑1 (MQ=255)
                                                                   tCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAacaaca                           <  2:1226279/93‑1 (MQ=255)
                                                                   tCTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAacaaca                           >  1:2493829/1‑93 (MQ=255)
                                                                    cTGGATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACAt                          <  1:3982036/93‑1 (MQ=255)
                                                                      ggATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATgc                        >  2:4009158/1‑93 (MQ=255)
                                                                       gATAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATgcg                       <  1:240940/93‑1 (MQ=255)
                                                                        aTAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATgcgc                      <  2:3333394/93‑1 (MQ=255)
                                                                         tAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGcgcg                     <  1:997396/93‑1 (MQ=255)
                                                                         tAGCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGcgcg                     <  1:4231727/93‑1 (MQ=255)
                                                                           gCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGc                   >  2:279658/1‑93 (MQ=255)
                                                                           gCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGc                   >  1:811705/1‑93 (MQ=255)
                                                                           gCATAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGc                   >  1:466848/1‑93 (MQ=255)
                                                                            cATAAATCACCCACCATGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCa                  <  1:2565308/93‑1 (MQ=255)
                                                                             aTAAATCACCCACCATGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAg                 >  2:5669396/1‑93 (MQ=255)
                                                                              tAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGa                >  2:3602448/1‑93 (MQ=255)
                                                                              tAAATCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGa                >  1:1571156/1‑93 (MQ=255)
                                                                               aaaTCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGAt               <  1:5614643/93‑1 (MQ=255)
                                                                               aaaTCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGAt               >  1:5120078/1‑93 (MQ=255)
                                                                                aaTCACCCACCATGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGTGGCAGATa              <  1:6472689/93‑1 (MQ=255)
                                                                                aaTCACCCACCATGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATa              >  2:2875427/1‑93 (MQ=255)
                                                                                aaTCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATa              <  1:5524909/93‑1 (MQ=255)
                                                                                  tCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATACa            <  1:4718339/93‑1 (MQ=255)
                                                                                  tCACCCACCACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATACa            >  1:2722334/1‑93 (MQ=255)
                                                                                    acccaccACGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATACATc          >  2:2192990/1‑93 (MQ=255)
                                                                                     cccaccATGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATACATCa         <  2:486695/93‑1 (MQ=255)
                                                                                       caccacGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGTGATGGGATCGGTATAGAACAACATGCGCGGCAGATACATCAGc       >  1:2321156/1‑93 (MQ=255)
                                                                                       caccacGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGACGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATACATCAGc       >  2:865606/1‑93 (MQ=255)
                                                                                        cccaCGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATACATCAGCa      <  1:3091721/92‑1 (MQ=255)
                                                                                        accaCGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATACATCAGCa      >  1:4629972/1‑93 (MQ=255)
                                                                                            tGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATACATCAGCAACAt  <  1:4578160/93‑1 (MQ=255)
                                                                                            |                                                                                            
AGCACCATGCCGAAGCACAGTTCACGCACAGTACGACCACGGGAGATGCGGGCGAGGAAGATACTCATCTGGATAGCATAAATCACCCACCATGCCCAGTAGAACACGGTCCAGCCCTGCGGGAAGCCGCCTTTAGCGATGGGATCGGTATAGAACAACATGCGCGGCAGATACATCAGCAACAT  >  NC_000913/40874‑41058

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