Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 383,897 C→T 14.0% M13I (ATG→ATA) ‡ yaiP ← putative glucosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913383,8970CT14.0% 184.6 / 16.7 71M13I (ATG→ATA) ‡yaiPputative glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base T (4/6);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATA  >  NC_000913/383748‑384046
                                                                                                                                                     |                                                                                                                                                     
cAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAAc                                                                                                                                                       <  1:833137/150‑1 (MQ=255)
cAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAAc                                                                                                                                                       >  1:1291719/1‑150 (MQ=255)
       cccGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTa                                                                                                                                                >  2:830863/1‑150 (MQ=255)
            cAAAAATAAGGGTTACGCAGTAGATTATCCAGGGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtt                                                                                                                                           <  1:1225993/150‑3 (MQ=255)
              aaaaTAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtttt                                                                                                                                         <  1:830863/150‑5 (MQ=255)
              aaaaTAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtttt                                                                                                                                         <  1:883968/150‑5 (MQ=255)
               aaaTAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATa                                                                                                                                        >  1:295022/1‑150 (MQ=255)
                aaTAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATAc                                                                                                                                       <  2:23356/150‑1 (MQ=255)
                 aTAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtttttaa                                                                                                                                      <  2:66822/150‑8 (MQ=255)
                       gTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaaa                                                                                                                                <  2:1292386/150‑1 (MQ=255)
                        ttACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAat                                                                                                                               >  2:546396/1‑150 (MQ=255)
                           cGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATc                                                                                                                            >  2:1032080/1‑150 (MQ=255)
                            gCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                                                                           <  2:1313230/150‑1 (MQ=255)
                            gCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                                                                           >  1:1313230/1‑150 (MQ=255)
                            gCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                                                                           <  1:1378470/150‑1 (MQ=255)
                              aGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAg                                                                                                                         <  2:184874/150‑1 (MQ=255)
                                     tATCCAGTTACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATTAATATCCTGGTTTTCa                                                                                                                  <  2:1076912/150‑1 (MQ=255)
                                      aTCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCAt                                                                                                                 <  2:218595/150‑1 (MQ=255)
                                        ccAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGc                                                                                                               <  1:1032080/150‑1 (MQ=255)
                                               cTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGc                                                                                                        >  1:524639/1‑150 (MQ=255)
                                                     cAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAg                                                                                                  <  2:333043/150‑1 (MQ=255)
                                                       gACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCt                                                                                                       >  1:1074797/1‑143 (MQ=255)
                                                       gACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCt                                                                                                       <  2:1074797/143‑1 (MQ=255)
                                                        acacGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAAt                                                                                               >  1:1152991/1‑150 (MQ=255)
                                                         cacGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATc                                                                                              <  2:119281/150‑1 (MQ=255)
                                                           cGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAt                                                                                            <  1:527723/150‑1 (MQ=255)
                                                              gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtttt                                                                                         >  1:368657/1‑150 (MQ=255)
                                                                  tCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGa                                                                                     <  1:86374/150‑1 (MQ=255)
                                                                    gTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAAt                                                                                   >  1:606424/1‑150 (MQ=255)
                                                                     ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAAtt                                                                                  <  1:302067/150‑1 (MQ=255)
                                                                      tataCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTa                                                                                 >  2:1395765/1‑150 (MQ=255)
                                                                           gCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAg                                                                            <  1:862342/150‑1 (MQ=255)
                                                                             aGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGcc                                                                          >  2:771859/1‑150 (MQ=255)
                                                                                tataATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCtta                                                                       >  1:1189574/1‑149 (MQ=255)
                                                                                 ataATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCa                                                                      <  1:771859/150‑1 (MQ=255)
                                                                                     tGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCa                                                                  <  2:295022/150‑1 (MQ=255)
                                                                                      ggCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAg                                                                 >  2:618960/1‑150 (MQ=255)
                                                                                       gCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa                                                                <  1:383104/150‑1 (MQ=255)
                                                                                         gTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc                                                              >  1:1174936/1‑150 (MQ=255)
                                                                                            aaTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCt                                                           <  2:398537/150‑1 (MQ=255)
                                                                                             aTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTa                                                          >  1:1421149/1‑150 (MQ=255)
                                                                                                cAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGTTTTTCATGCAATTTGCTAAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAttt                                                       >  1:725089/1‑149 (MQ=255)
                                                                                                  gCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAc                                                     >  1:1230980/1‑150 (MQ=255)
                                                                                                   ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc                                                    <  1:1323436/150‑1 (MQ=255)
                                                                                                    ctttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCt                                                   >  1:81956/1‑150 (MQ=255)
                                                                                                        ctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCa                                                                      <  1:595111/127‑1 (MQ=255)
                                                                                                        ctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCtgttg                                               <  1:1373795/150‑1 (MQ=255)
                                                                                                        ctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGATCAAAGAATCATTTTATGAATTACAAAGCCTTCa                                                                      >  2:595111/1‑127 (MQ=255)
                                                                                                             gACTCGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCcgc                                          >  1:1020581/1‑150 (MQ=255)
                                                                                                                 gggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAGAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGaaa                                      >  1:388167/1‑150 (MQ=255)
                                                                                                                 gggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAGAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGaaa                                      >  1:394794/1‑150 (MQ=255)
                                                                                                                  ggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAAt                                     <  1:1086583/150‑1 (MQ=255)
                                                                                                                    tttACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATtg                                   >  2:1021761/1‑150 (MQ=255)
                                                                                                                      tACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATtgtg                                 >  2:1305120/1‑150 (MQ=255)
                                                                                                                       aCGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGc                                <  1:892865/150‑1 (MQ=255)
                                                                                                                         gTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCaa                              >  1:1357477/1‑150 (MQ=255)
                                                                                                                             tAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGc                          <  1:666410/150‑1 (MQ=255)
                                                                                                                                    cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAgg                   <  2:1152991/150‑1 (MQ=255)
                                                                                                                                    cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAgg                   >  1:696217/1‑150 (MQ=255)
                                                                                                                                     tgtatgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGc                  >  1:272233/8‑150 (MQ=255)
                                                                                                                                      tAAAAACCATAGCAACATTGCTATGGACATACTGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCg                 >  1:414617/1‑150 (MQ=255)
                                                                                                                                       aaaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc                                                    <  1:874037/114‑1 (MQ=255)
                                                                                                                                       aaaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc                                                    >  2:874037/1‑114 (MQ=255)
                                                                                                                                        aaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGAt               >  1:1252331/1‑150 (MQ=255)
                                                                                                                                          gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCa             >  1:997516/3‑150 (MQ=255)
                                                                                                                                             cATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGa          >  2:974652/1‑150 (MQ=255)
                                                                                                                                              aTAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc                                                              <  2:608830/97‑1 (MQ=255)
                                                                                                                                              aTAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc                                                              >  1:608830/1‑97 (MQ=255)
                                                                                                                                               tAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACg        <  2:282447/150‑1 (MQ=255)
                                                                                                                                               tAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACg        <  2:282446/150‑1 (MQ=255)
                                                                                                                                               tAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGAAACTGCACGCAGGCGATCATGACg        >  1:621502/1‑150 (MQ=255)
                                                                                                                                                 gCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCa      <  2:166931/150‑1 (MQ=255)
                                                                                                                                                    aCATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCatat   >  2:226257/1‑150 (MQ=255)
                                                                                                                                                     cATAGCTATTGACATACATATAAATATCAATTTTTTAATGCAATTTTATCAAAGAATAATTTTATGAATTACAAAGAATTCAAAAAGATCTATATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCAtata  >  2:151089/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATATA  >  NC_000913/383748‑384046

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