Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,875,564 C→T 6.4% G256S (GGT→AGT)  cysD ← sulfate adenylyltransferase, subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,875,5640CT 6.4% 196.5 / 2.6 62G256S (GGT→AGT) cysDsulfate adenylyltransferase, subunit 2
Reads supporting (aligned to +/- strand):  ref base C (34/24);  new base T (2/2);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

ATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGATATAG  >  NC_000913/2875418‑2875708
                                                                                                                                                  |                                                                                                                                                
aTCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCTGTCAATCTCGCGGCCCGGACGTTGACTTGTGGTGGAAACCAGCTTCTCTTCGATGATTTCCGGCTGTGTTTGTGCATTTGACTCGACCGCACCgg                                                                                                                                                <  6:132769/149‑1 (MQ=255)
aTCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCgg                                                                                                                                                >  8:24270/1‑149 (MQ=255)
aTCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGAAGTGTTTGTGCATTTGACTCcaccgcaccg                                                                                                                                                 <  8:245445/148‑1 (MQ=255)
  cTTAAAATTACGCGTGTCGTTTTTTCAGGTGCATAGACCCCGACTTGTCGCGGTCAATCACGCGTCACTGACGTTCAATGGTGGTGGAATCGAGGTTCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGTCTAaaccgcaccg                                                                                                                                                 <  4:307634/146‑1 (MQ=255)
   ttAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGt                                                                                                                                               >  4:203250/1‑147 (MQ=255)
    tAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg                                                                                                                                            >  6:250766/1‑149 (MQ=255)
               cTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGGCGCGGGCAATCACGCGGCCCTGACGTTTACTGGTTGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCagcagc                                                                                                                                 <  6:180567/149‑1 (MQ=255)
                tGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCCTCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGcc                                                                                                                                >  4:126167/1‑149 (MQ=255)
                  aCGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGcc                                                                                                                                >  1:173063/1‑147 (MQ=255)
                  aCGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGcc                                                                                                                                <  2:173063/147‑1 (MQ=255)
                    gTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGc                                                                                                                            >  5:83069/1‑149 (MQ=255)
                     ttttttCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCg                                                                                                                           >  6:228642/1‑149 (MQ=255)
                             gCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGcaa                                                                                                                   >  5:251222/1‑149 (MQ=255)
                               tCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACg                                                                                                                 >  5:235498/1‑149 (MQ=255)
                                 cATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCa                                                                                                               <  6:235498/149‑1 (MQ=255)
                                 cATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCa                                                                                                               >  4:141921/1‑149 (MQ=255)
                                  aTAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCCGCGGCCAGCAGCCCAGCGTACGGAAACGCAc                                                                                                              >  3:220532/1‑149 (MQ=255)
                                  aTAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGccc                                                                                                              >  5:126193/1‑147 (MQ=255)
                                         ccGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCa                                                                                                                              <  6:134292/126‑1 (MQ=255)
                                           gCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGttt                                                                                                     <  7:218413/149‑1 (MQ=255)
                                            ccTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGtttt                                                                                                    >  5:140483/1‑149 (MQ=255)
                                             cTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGttttt                                                                                                   <  6:140483/149‑1 (MQ=255)
                                                         aaTCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCATCCCAGCGTACGGAAACGCACCATCCGTTTTTTaa                                                                                                <  6:235463/140‑1 (MQ=255)
                                                         aaTCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTaa                                                                                                >  5:235463/1‑140 (MQ=255)
                                                             aCGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCg                                                                                    <  8:55635/148‑1 (MQ=255)
                                                              cgcgGACCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGAATTTGACTCCACCGCACCGGTCAGCGGCAAGCAGCCCAGCGTACGGAAACGCACCATCCGGATTTTAATCACTACGCCCGGt                                                                                  >  6:64133/1‑149 (MQ=255)
                                                                   cccTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAg                                                                             <  7:24270/149‑1 (MQ=255)
                                                                     cTGACGTTCACTGGGGGTGTAAACGAGCATCTCTTCGTTGATTTCCGGCAGTGTTTGGGCATTTGACTCCACCGCACCGGTCAGCGGCGAGCAGCCCAGCGTACGGAAACGGTCCATCCGTTTTTTTATCACTTCGCCCGGTTGCAggg                                                                           <  4:100643/149‑2 (MQ=255)
                                                                        aCGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCTGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGa                                                                        <  2:205549/149‑1 (MQ=255)
                                                                          gTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATa                                                                      >  6:17727/1‑149 (MQ=255)
                                                                             cACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACgg                                                                   >  8:57655/1‑149 (MQ=255)
                                                                              aCTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACgg                                                                   <  7:23958/148‑1 (MQ=255)
                                                                                  gtggtggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTg                                                                <  7:216733/147‑1 (MQ=255)
                                                                                  gtggtggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGtc                                                              <  5:17727/149‑1 (MQ=255)
                                                                                   tggtggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACTGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGtca                                                             >  6:230583/1‑149 (MQ=255)
                                                                                   tggtggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACTGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGtca                                                             >  6:235331/1‑149 (MQ=255)
                                                                                   tggtggAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGtca                                                             >  6:323654/1‑149 (MQ=255)
                                                                                        gAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCatcaat                                                        <  5:74896/149‑1 (MQ=255)
                                                                                        gAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCatcaa                                                         >  7:296667/1‑148 (MQ=255)
                                                                                        gAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCatcaa                                                         <  3:203250/148‑1 (MQ=255)
                                                                                        gAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCatcaa                                                         <  3:141921/148‑1 (MQ=255)
                                                                                         aaaCCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCatcaat                                                        >  3:121608/1‑148 (MQ=255)
                                                                                           aCCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCCTCCGGTTTTTAATCACTTCGCCCGGTTGCAGGTCGAAACGGTAGTCATCAatcat                                                     >  5:235007/1‑149 (MQ=255)
                                                                                           aCCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAatcat                                                     <  5:250766/149‑1 (MQ=255)
                                                                                            ccAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAatcatc                                                    <  5:228642/149‑1 (MQ=255)
                                                                                             cAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAatcatca                                                   >  4:174374/1‑149 (MQ=255)
                                                                                              aGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAatcatcat                                                  >  2:123579/1‑149 (MQ=255)
                                                                                                  tctcTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCaa                                               >  7:68347/1‑148 (MQ=255)
                                                                                                     ctTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATa                                           >  3:9809/1‑149 (MQ=255)
                                                                                                        cGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                          >  6:249568/1‑115 (MQ=255)
                                                                                                        cGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                          >  8:134305/1‑115 (MQ=255)
                                                                                                        cGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                          <  7:134305/115‑1 (MQ=255)
                                                                                                        cGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTc                                                                          <  5:249568/115‑1 (MQ=255)
                                                                                                           tgatTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTcg                                     >  7:105975/1‑149 (MQ=255)
                                                                                                              tttCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGt                                  >  6:89082/1‑149 (MQ=255)
                                                                                                                    gCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGa                            <  4:9809/149‑1 (MQ=255)
                                                                                                                        tgtTTGTGCATTTGACTCAACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCg                        >  3:64920/1‑149 (MQ=255)
                                                                                                                             gtgCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGt                   >  2:191742/1‑149 (MQ=255)
                                                                                                                                cATTTGACTCCACCGCACTGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTc                 <  5:230583/148‑1 (MQ=255)
                                                                                                                                cATTTGACTCCACCGCACTGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTc                 <  5:235331/148‑1 (MQ=255)
                                                                                                                                 aTTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTATCATCAATCATCATCAACATACCGTCGCGTTCCAGAAcc                         <  7:7447/139‑1 (MQ=255)
                                                                                                                                 aTTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTATCATCAATCATCATCAACATACCGTCGCGTTCCAGAAcc                         >  8:7447/1‑139 (MQ=255)
                                                                                                                                   ttGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCg               >  7:137604/1‑147 (MQ=255)
                                                                                                                                       cTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCAGCAATCATCCTCAACATCCCGTCGCGTTCCAGCACCGGACGTTCCGCAGCga         >  4:253366/1‑149 (MQ=255)
                                                                                                                                          caccgcaccgGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACgg                                                                   <  2:34466/88‑1 (MQ=255)
                                                                                                                                          caccgcaccgGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACgg                                                                   >  1:34466/1‑88 (MQ=255)
                                                                                                                                          caccgcaccgGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgag        <  1:330394/147‑1 (MQ=255)
                                                                                                                                           accgcaccgGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGata     <  3:174374/149‑1 (MQ=255)
                                                                                                                                            ccgcaccgGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGatat    <  1:191742/149‑1 (MQ=255)
                                                                                                                                              gcaccgGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGGTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGATATAg  <  3:253366/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCACTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCGGCCAGCAGCCCAGCGTACGGAAACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGATATAG  >  NC_000913/2875418‑2875708

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑