Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,430,144 G→C 6.2% A32G (GCC→GGC)  purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,430,1440GC 6.2% 267.8 / 6.3 80A32G (GCC→GGC) purFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (32/43);  new base C (2/3);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCAAGA  >  NC_000913/2429996‑2430291
                                                                                                                                                    |                                                                                                                                                   
ggAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCcggc                                                                                                                                                    >  1:393649/1‑150 (MQ=255)
  aGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCcggcgg                                                                                                                                                  <  1:723688/150‑1 (MQ=255)
   gCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCggcggc                                                                                                                                                 <  2:1283131/150‑1 (MQ=255)
   gCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCggcggc                                                                                                                                                 >  1:1490585/1‑150 (MQ=255)
     tGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCcgg                                                                                                                                                     <  1:237089/144‑1 (MQ=255)
     tGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCcgg                                                                                                                                                     >  2:237089/1‑144 (MQ=255)
         aGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATcatc                                                                                                                                           <  1:77113/150‑4 (MQ=255)
         aGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATcatc                                                                                                                                           <  2:284769/150‑4 (MQ=255)
          gCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGCCGGCATcatca                                                                                                                                          <  2:831208/150‑5 (MQ=255)
              tGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACcgc                                                                                                                                      <  1:455346/150‑1 (MQ=255)
               ggggTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCTCTTTACGCAAACGGAAGCAGTTATTGGTATCTATGTTGATTATGTCGGCGTCATCGTGACcgcg                                                                                                                                     <  2:146779/150‑1 (MQ=255)
               ggggTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACcgcg                                                                                                                                     >  2:675223/1‑150 (MQ=255)
                 ggTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGAt                                                                                                                                   <  1:531/150‑1 (MQ=255)
                        cacaTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAg                                                                                                                            >  2:1138948/1‑150 (MQ=255)
                        cacaTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAg                                                                                                                            >  1:454998/1‑150 (MQ=255)
                        cacaTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAg                                                                                                                            <  1:1474415/150‑1 (MQ=255)
                        cacaTGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAg                                                                                                                            <  1:1051438/150‑1 (MQ=255)
                             gACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCg                                                                                                                       >  1:727505/1‑150 (MQ=255)
                               ccAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGtt                                                                                                                     <  2:763316/150‑1 (MQ=255)
                                cAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTa                                                                                                                    <  2:495829/150‑1 (MQ=255)
                                   tGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAgg                                                                                                                 <  2:199375/150‑1 (MQ=255)
                                   tGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAgg                                                                                                                 <  2:205157/150‑1 (MQ=255)
                                     cccATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGca                                                                                                               <  1:292423/150‑1 (MQ=255)
                                      ccATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGcat                                                                                                              >  1:1608065/1‑150 (MQ=255)
                                                cTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                    >  2:1189604/1‑150 (MQ=255)
                                                cTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                    >  1:1035004/1‑150 (MQ=255)
                                                cTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                    <  2:580156/150‑1 (MQ=255)
                                                cTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                    >  1:1406485/1‑150 (MQ=255)
                                                cTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGa                                                                                                    >  2:78406/1‑150 (MQ=255)
                                                   cAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTg                                                                                                 >  2:1325121/1‑150 (MQ=255)
                                                    aaaCGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTgg                                                                                                <  2:1087436/150‑1 (MQ=255)
                                                     aaCGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGt                                                                                               <  2:546877/150‑1 (MQ=255)
                                                     aaCGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGt                                                                                               <  1:658987/150‑1 (MQ=255)
                                                       cGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCAGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGAAGCCGGCGGCATCCAGACTGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTa                                                                                             >  2:650262/1‑150 (MQ=255)
                                                        gCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGTTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTaa                                                                                            <  2:550719/150‑1 (MQ=255)
                                                           gCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCg                                                                                         <  1:1138948/150‑1 (MQ=255)
                                                              tatGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCa                                                                                      <  1:959238/150‑1 (MQ=255)
                                                               atGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCAt                                                                                     <  2:1500661/150‑1 (MQ=255)
                                                                tGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATa                                                                                    <  1:1189604/150‑1 (MQ=255)
                                                                tGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATa                                                                                    <  1:78406/150‑1 (MQ=255)
                                                                   cGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAaca                                                                                 >  2:146480/1‑150 (MQ=255)
                                                                      gCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCg                                                                              >  2:1059351/1‑150 (MQ=255)
                                                                          cAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGa                                                                          >  2:630461/1‑150 (MQ=255)
                                                                          cAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACAGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGa                                                                          >  2:939958/1‑150 (MQ=255)
                                                                            aaTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATa                                                                        <  1:1573460/150‑1 (MQ=255)
                                                                            aaTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATa                                                                        >  2:1401942/1‑150 (MQ=255)
                                                                             aTACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATAc                                                                       <  1:67725/150‑1 (MQ=255)
                                                                                cATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGa                                                                    >  1:163489/1‑150 (MQ=255)
                                                                                  tCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACa                                                                  >  1:1268742/1‑150 (MQ=255)
                                                                                         ccAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGc                                                           >  2:1340632/1‑150 (MQ=255)
                                                                                         ccAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGc                                                           <  1:1401942/150‑1 (MQ=255)
                                                                                          gcttccGTTTGCGTAAAGCGAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACTGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGca                                                          >  2:884113/5‑150 (MQ=255)
                                                                                          ccttccGTTTGCGTAAAGCGAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGca                                                          >  1:625037/5‑150 (MQ=255)
                                                                                          cAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGca                                                          <  2:454998/150‑1 (MQ=255)
                                                                                          cAGCCCGTTCGCTTTACGAAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGca                                                          <  1:1257888/150‑1 (MQ=255)
                                                                                           cttccGTTTGCGTAAAGCGAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGcac                                                         >  1:1201842/4‑150 (MQ=255)
                                                                                           cttccGTTTGCGTAAAGCGAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGcac                                                         >  1:842204/4‑150 (MQ=255)
                                                                                                   cGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCAGCATAACACCGGCGATACCGACAATACCGCACATACGTCt                                                 >  2:1267276/1‑150 (MQ=255)
                                                                                                   cGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCGGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCt                                                 <  1:1372012/150‑1 (MQ=255)
                                                                                                      tttACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGTCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATTCCGACAATACCGCACATACGTCtttt                                              <  1:855580/150‑1 (MQ=255)
                                                                                                       ttACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCttttt                                             <  1:1007873/150‑1 (MQ=255)
                                                                                                            cAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCg                                        <  1:18878/150‑1 (MQ=255)
                                                                                                            cAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCg                                        <  2:1011388/150‑1 (MQ=255)
                                                                                                            cAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCg                                        <  2:1406485/150‑1 (MQ=255)
                                                                                                              aaCGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGtt                                      <  2:1394172/150‑1 (MQ=255)
                                                                                                                      cAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCcaca                              <  1:455274/150‑1 (MQ=255)
                                                                                                                        gTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACAtc                            >  2:632043/1‑150 (MQ=255)
                                                                                                                          tATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACAtctc                          >  1:142837/1‑150 (MQ=255)
                                                                                                                          tATTGGCATCTATGGTGATCATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGGTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACAtctc                          <  1:666360/150‑1 (MQ=255)
                                                                                                                              ggCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCagag                      >  1:649410/1‑150 (MQ=255)
                                                                                                                              ggCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCagag                      <  1:146480/150‑1 (MQ=255)
                                                                                                                                cATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCa                    <  2:1035004/150‑1 (MQ=255)
                                                                                                                                  tCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACt                  <  1:548518/150‑1 (MQ=255)
                                                                                                                                   cTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACtt                 <  2:636695/150‑1 (MQ=255)
                                                                                                                                    tATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTa                <  2:28752/150‑1 (MQ=255)
                                                                                                                                    tATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTa                <  2:163489/150‑1 (MQ=255)
                                                                                                                                    tATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACAGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTa                >  2:501147/1‑150 (MQ=255)
                                                                                                                                      tGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACg              >  1:1510452/1‑150 (MQ=255)
                                                                                                                                      cggtcaggatgCCGCCGGCATCATCACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACg              >  1:689145/8‑150 (MQ=255)
                                                                                                                                          caggatgCCGCCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCACAGAGCACTTACGctct          >  2:1176139/4‑150 (MQ=255)
                                                                                                                                           atgatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTg         <  2:625037/150‑1 (MQ=255)
                                                                                                                                             gatgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTggg       <  2:842204/150‑1 (MQ=255)
                                                                                                                                               tgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCa     >  2:1593322/1‑150 (MQ=255)
                                                                                                                                               tgCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCa     >  1:1379920/1‑150 (MQ=255)
                                                                                                                                                gCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCaa    >  1:683373/1‑150 (MQ=255)
                                                                                                                                                 ccGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCAAg   <  2:649410/150‑1 (MQ=255)
                                                                                                                                                  cggcggCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCAAGa  <  2:1608065/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGAGCTGCCAGCCGTGGGGTAACGCACATGACCAATGCCCATATTGCCCTGCAAACGCTGCATATGGCGAGCTTCAAATACATCGCTCACCAGCCCGTTCGCTTTACGCAAACGGAAGCAGTTATTGGCATCTATGGTGATGATGCCGGCGGCATCCTGACCGCGATGCTGAAGCACCGTTAAGGCATCATAAATCGACTGGTTAACCGGCATAACACCGGCGATACCGACAATACCGCACATACGTCTTTTTCCTCGTTAAGCCACATCTCAGAGCACTTACGCTCTGGGCAAGA  >  NC_000913/2429996‑2430291

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: