Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,685,527 A→G 50.6% G354G (GGT→GGC fumC ← fumarate hydratase (fumarase C),aerobic Class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,685,5270AG50.6% ‑6.3 / 112.2 75G354G (GGT→GGCfumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base A (23/14);  new base G (19/19);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

ATTGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAAT  >  NC_000913/1685380‑1685673
                                                                                                                                                   |                                                                                                                                                  
aTTGAGTAATTGATCGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTa                                                                                                                                                    <  1:779414/148‑1 (MQ=255)
 ttGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTAcc                                                                                                                                                  <  3:452972/149‑1 (MQ=255)
   gAGTAATTGATTGATTCGCTCACTATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCgg                                                                                                                                                <  3:510530/149‑1 (MQ=255)
     gTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGaa                                                                                                                                              <  3:597999/149‑1 (MQ=255)
     gTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGa                                                                                                                                               <  2:89113/148‑1 (MQ=255)
       aaTTGATTGATTCGCTCACGATTCGGTTCAACACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGCACGGAAGACGTTCAGTTCAAAGTTACCGGAAg                                                                                                                                             <  1:277773/148‑1 (MQ=255)
           gattgattCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGccc                                                                                                                                        <  4:819029/149‑1 (MQ=255)
            attgattCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                                        >  1:561111/1‑148 (MQ=255)
                attCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCAt                                                                                                                                   >  4:827099/1‑149 (MQ=255)
                 ttCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCAt                                                                                                                                   >  1:261125/1‑148 (MQ=255)
                 ttCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCAt                                                                                                                                   >  1:62435/1‑148 (MQ=255)
                 ttCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGACAGCAAGCGCACCGATCGCAATACATTGTGGATCACCATTGGACGGAAGACGTTCAGCTAAAAGTTACCGGAAGCGCCCCCCAt                                                                                                                                   >  2:372365/1‑148 (MQ=255)
                  tCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGt                                                                                                                                 <  3:380871/149‑1 (MQ=255)
                     cTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGa                                                                                                                              <  3:35665/149‑1 (MQ=255)
                      tCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGa                                                                                                                              <  2:561111/148‑1 (MQ=255)
                       cACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCCCCGATTGCAGGTAATTGTGGATCACCATTGGACGGAAGACGTTCTGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGAt                                                                                                                             >  2:321182/1‑148 (MQ=255)
                       cACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGAt                                                                                                                             >  1:542255/1‑148 (MQ=255)
                           aTTCGGTTCAATACCCACTGCGCAGTTTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCa                                                                                                                        >  3:485718/1‑149 (MQ=255)
                             tCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAc                                                                                                                       >  1:479375/1‑148 (MQ=255)
                                gTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCAcgtcg                                                                                                                   >  3:239729/1‑149 (MQ=255)
                                   cAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGtt                                                                                                                 <  2:947914/148‑1 (MQ=255)
                                     aTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                              <  4:24231/149‑1 (MQ=255)
                                      tACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                              >  1:700491/1‑148 (MQ=255)
                                      tACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTAGCCGAAAGCGCCCCCCATGTTGATCGCCACGTCGTTcccc                                                                                                             >  4:622933/1‑149 (MQ=255)
                                        cccACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                           >  3:141625/1‑149 (MQ=255)
                                          cACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCa                                                                                                         >  4:944402/1‑149 (MQ=255)
                                          cACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCAGCa                                                                                                         >  3:842334/1‑149 (MQ=255)
                                           aCTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCAc                                                                                                        >  3:662393/1‑149 (MQ=255)
                                               cgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                                                                    <  4:366557/149‑1 (MQ=255)
                                                    tgtttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATTACCTGAcagc                                                                                                <  1:718391/148‑1 (MQ=255)
                                                      tttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCaga                                                                                             <  3:468840/149‑1 (MQ=255)
                                                       ttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGa                                                                                                                                               <  2:572918/98‑1 (MQ=255)
                                                       ttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGa                                                                                                                                               >  1:572918/1‑98 (MQ=255)
                                                       ttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCaga                                                                                             >  2:106081/1‑148 (MQ=255)
                                                       ttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCaga                                                                                             >  2:134622/1‑148 (MQ=255)
                                                          ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCGCCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                         >  4:185191/1‑149 (MQ=255)
                                                          ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                         >  4:855064/1‑149 (MQ=255)
                                                           tAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                          >  1:945719/1‑147 (MQ=255)
                                                           tAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                                                                          <  2:945719/147‑1 (MQ=255)
                                                           tAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                         >  1:58005/1‑148 (MQ=255)
                                                           tAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                         >  2:305356/1‑148 (MQ=255)
                                                           tAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                                                                         >  2:761421/1‑148 (MQ=255)
                                                            aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATg                                                                                        >  2:557747/1‑148 (MQ=255)
                                                             aaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                                                                      <  3:330008/149‑1 (MQ=255)
                                                             aaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                                                                      <  4:485718/149‑1 (MQ=255)
                                                             aaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                                                                      >  4:410003/1‑149 (MQ=255)
                                                              aaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                                                                      <  1:820382/148‑1 (MQ=255)
                                                                 tttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAAt                                                                                  <  4:604092/149‑1 (MQ=255)
                                                                  ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTg                                                                                                     >  1:664811/1‑129 (MQ=255)
                                                                  ttCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTg                                                                                                     <  2:664811/129‑1 (MQ=255)
                                                                     cATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                                                              <  4:600606/149‑1 (MQ=255)
                                                                       tGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTc                                                                             >  2:476081/1‑148 (MQ=255)
                                                                              tGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACtgt                                                                      >  2:963917/1‑148 (MQ=255)
                                                                                ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTgtgt                                                                    <  1:175713/148‑1 (MQ=255)
                                                                                ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTgtgt                                                                    <  1:476081/148‑1 (MQ=255)
                                                                                      aGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGtt                                                              >  1:80323/1‑148 (MQ=255)
                                                                                            ccGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc                                                       >  3:876209/1‑149 (MQ=255)
                                                                                                  gCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                                  <  2:542255/148‑1 (MQ=255)
                                                                                                       aaaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGa                                             >  1:415597/1‑148 (MQ=255)
                                                                                                          ttGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTg                                          <  1:755511/148‑1 (MQ=255)
                                                                                                           tgtgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGa                                         <  2:821583/148‑1 (MQ=255)
                                                                                                                aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGcc                                   <  4:239729/149‑1 (MQ=255)
                                                                                                                 tCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGcc                                   >  1:204199/1‑148 (MQ=255)
                                                                                                                 tCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGcc                                   >  2:658596/1‑148 (MQ=255)
                                                                                                                    ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCgg                                >  1:406530/1‑148 (MQ=255)
                                                                                                                              gAAGACGTTCAGTTCAAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTcc                      <  1:527370/148‑1 (MQ=255)
                                                                                                                              gAAGACGCCCAGTTCCAAGTTGCCGGGAGCGCCCCCGATGTTGATCGCAACGTCGTCCCCCACCACCTGACCGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTCCCCTTTCCCCGGCATGATTCAGCTGCCCGGCTCATTTTcc                      <  2:765949/148‑1 (MQ=255)
                                                                                                                                     ttcagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTga              <  3:530476/149‑1 (MQ=255)
                                                                                                                                     ttcagttcaAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTg               <  2:808778/148‑1 (MQ=255)
                                                                                                                                      tcagttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTga              >  2:933844/1‑148 (MQ=255)
                                                                                                                                        agttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGAtt                                           <  3:79951/117‑1 (MQ=255)
                                                                                                                                        agttcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGAtt                                           >  4:79951/1‑117 (MQ=255)
                                                                                                                                           tcaAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAttt         >  2:428265/1‑148 (MQ=255)
                                                                                                                                            caAAGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTGCa       >  3:635272/1‑149 (MQ=255)
                                                                                                                                              aaGTTGCCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCAcc     >  4:839528/1‑149 (MQ=255)
                                                                                                                                                  tACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAAt  >  1:959019/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTGAGTAATTGATTGATTCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCACCAAT  >  NC_000913/1685380‑1685673

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: