Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,144,732 G→C 10.4% Q246H (CAG→CAC pflD → putative formate acetyltransferase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,144,7320GC10.4% 297.4 / 21.7 96Q246H (CAG→CACpflDputative formate acetyltransferase 2
Reads supporting (aligned to +/- strand):  ref base G (45/41);  new base C (6/4);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

GAAGCCTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAAC  >  NC_000913/4144583‑4144880
                                                                                                                                                     |                                                                                                                                                    
gAAGCCTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAg                                                                                                                                                      <  1:1566955/150‑1 (MQ=255)
 aaGCCTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAgg                                                                                                                                                     <  1:945547/150‑1 (MQ=255)
 aaGCCTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAgg                                                                                                                                                     <  2:805271/150‑1 (MQ=255)
     cTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTg                                                                                                                                                 <  1:1487727/150‑1 (MQ=255)
       cGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGcc                                                                                                                                               >  1:154097/1‑150 (MQ=255)
       cGCAGAAACACATTTTGCGTTACGCCGAACAGGCGGAAACGATGGCGGCAAACTGCACACATGCCCCCCCTCGCGAAGAGCTGCTGACTAATGCAGCGATCACCCGCCAAAACACCCAACATAAGCCGCAGACGATCTGGCAGCCGAGcc                                                                                                                                               >  2:1292835/1‑150 (MQ=255)
         cAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAg                                                                                                                                             >  2:544169/1‑150 (MQ=255)
           gAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaa                                                                                                                                           <  1:1105706/150‑3 (MQ=255)
            aaaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTa                                                                                                                                          <  2:1120319/150‑1 (MQ=255)
            aaaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTa                                                                                                                                          <  1:838822/150‑1 (MQ=255)
             aaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATTCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaacg                                                                                                                                         <  1:1538964/150‑5 (MQ=255)
             aaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaacg                                                                                                                                         <  2:1582806/150‑5 (MQ=255)
             aaCACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCACGCCTGCCAgaacg                                                                                                                                         <  2:989583/150‑5 (MQ=255)
              acacATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTAtt                                                                                                                                        >  1:103190/1‑150 (MQ=255)
                  aTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTgg                                                                                                                                    <  1:248778/150‑1 (MQ=255)
                  aTTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTgg                                                                                                                                    >  1:1517218/1‑150 (MQ=255)
                    tttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTa                                                                                                                                  >  2:1160446/1‑150 (MQ=255)
                    tttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTa                                                                                                                                  <  1:810239/150‑1 (MQ=255)
                    tttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTa                                                                                                                                  <  1:616090/150‑1 (MQ=255)
                     ttGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTAc                                                                                                                                 >  1:679079/1‑150 (MQ=255)
                        cGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATg                                                                                                                              <  1:986963/150‑1 (MQ=255)
                        cGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATg                                                                                                                              <  1:449034/150‑1 (MQ=255)
                        cGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATg                                                                                                                              <  1:1152513/150‑1 (MQ=255)
                         gTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGa                                                                                                                             >  1:1377705/1‑150 (MQ=255)
                          ttACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGaa                                                                                                                            <  2:239018/150‑1 (MQ=255)
                            aCGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAca                                                                                                                          <  2:803179/150‑1 (MQ=255)
                              gCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcatc                                                                                                                        >  1:1459967/1‑150 (MQ=255)
                              gCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAAcatc                                                                                                                        <  1:1013311/150‑1 (MQ=255)
                                     tGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGc                                                                                                                 <  1:1242439/150‑1 (MQ=255)
                                     tGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGc                                                                                                                 <  1:182228/150‑1 (MQ=255)
                                      ggCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCa                                                                                                                <  2:1214083/150‑1 (MQ=255)
                                        gggAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAAt                                                                                                              <  2:738822/149‑1 (MQ=255)
                                        cGGAAACGGTGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAAt                                                                                                              >  1:605443/1‑150 (MQ=255)
                                          gAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATAc                                                                                                            <  2:457840/150‑1 (MQ=255)
                                           aaaCGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACg                                                                                                           >  1:1166804/1‑150 (MQ=255)
                                               gttGGCCGCAAACTGCACAGATCCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGATATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATc                                                                                                       <  2:1616762/148‑1 (MQ=255)
                                                     ggcAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTCCATGCACATCAGTCTGCAAGACGAATcccccgc                                                                                                 >  1:1498275/1‑145 (MQ=255)
                                                          aCTGCACAGATGCCCAGCGTCGCGAAGCGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGtt                                                                                            >  2:158568/1‑150 (MQ=255)
                                                           cTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTc                                                                                           >  2:365860/1‑150 (MQ=255)
                                                            tGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCg                                                                                          >  1:967078/1‑150 (MQ=255)
                                                             gCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGc                                                                                         >  1:1076458/1‑150 (MQ=255)
                                                                caGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTAt                                                                                      >  2:214670/1‑150 (MQ=255)
                                                                  gATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCg                                                                                    <  1:1481684/150‑1 (MQ=255)
                                                                      cccAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTgg                                                                                <  2:1236487/150‑1 (MQ=255)
                                                                       ccAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTggg                                                                               >  1:568395/1‑150 (MQ=255)
                                                                            gTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCt                                                                          >  1:786060/1‑150 (MQ=255)
                                                                              cgcgAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTc                                                                        <  1:514460/150‑1 (MQ=255)
                                                                              cgcgAACAGCTGATGACTATTGCAGAGATCTCCCGCAATAACGCGCAACATAACCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTc                                                                        >  2:1322083/1‑150 (MQ=255)
                                                                               gcgAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACTTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCg                                                                       >  2:1120726/1‑150 (MQ=255)
                                                                                 gAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGAc                                                                     <  1:50733/150‑1 (MQ=255)
                                                                                 gAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGAc                                                                     >  1:764111/1‑150 (MQ=255)
                                                                                 gAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGAc                                                                     >  1:970200/1‑150 (MQ=255)
                                                                                            gACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGtt                                                          >  2:998998/1‑150 (MQ=255)
                                                                                              cTATTGCAGAGATCTCACGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCCTTCTGCAATACGAATCCAACGCCAGTTCGCTATAGTTGGGGCGCTTCGACCAGTATATGTTGc                                                        >  1:1266884/1‑150 (MQ=255)
                                                                                               tATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGcc                                                       <  1:1120726/150‑1 (MQ=255)
                                                                                                aTTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCg                                                      >  1:1494082/1‑150 (MQ=255)
                                                                                                aTTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCg                                                      <  2:1166804/150‑1 (MQ=255)
                                                                                                   gCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCATTACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTc                                                   >  1:1514690/1‑150 (MQ=255)
                                                                                                   gCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTc                                                   <  1:365860/150‑1 (MQ=255)
                                                                                                    cAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCt                                                  >  1:1260522/1‑150 (MQ=255)
                                                                                                     agagATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTa                                                 >  1:123228/1‑150 (MQ=255)
                                                                                                      gagaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTAt                                                >  1:507090/1‑150 (MQ=255)
                                                                                                        gaTCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCa                                              <  2:171667/150‑1 (MQ=255)
                                                                                                            tcCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAca                                          >  2:1449233/1‑150 (MQ=255)
                                                                                                              ccGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcatc                                        <  1:998998/150‑1 (MQ=255)
                                                                                                              ccGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcatc                                        >  2:331520/1‑150 (MQ=255)
                                                                                                               cGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGAcatca                                       >  2:169893/1‑150 (MQ=255)
                                                                                                                    tAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAAc                                  <  2:967078/150‑1 (MQ=255)
                                                                                                                       cgcgCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCa                               >  2:38276/1‑150 (MQ=255)
                                                                                                                       cgcgCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCa                               >  1:1365092/1‑150 (MQ=255)
                                                                                                                       cgcgCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCa                               >  1:186285/1‑150 (MQ=255)
                                                                                                                       cgcgCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCa                               >  1:812408/1‑150 (MQ=255)
                                                                                                                          gcAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAggg                            <  1:1380636/150‑1 (MQ=255)
                                                                                                                              cATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGaa                        <  2:1377705/150‑1 (MQ=255)
                                                                                                                               aTAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAg                       <  1:544169/150‑1 (MQ=255)
                                                                                                                                tAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGa                      >  1:393139/1‑150 (MQ=255)
                                                                                                                                tAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGa                      >  2:231239/1‑150 (MQ=255)
                                                                                                                                   gCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgc                   <  2:1517218/150‑1 (MQ=255)
                                                                                                                                    ccGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcg                  >  1:566333/1‑150 (MQ=255)
                                                                                                                                    ccGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcg                  <  2:786060/150‑1 (MQ=255)
                                                                                                                                     cGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcgg                 <  1:1591096/150‑1 (MQ=255)
                                                                                                                                      gCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATGATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcggc                <  1:522986/150‑1 (MQ=255)
                                                                                                                                      gCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATgcggc                <  2:970200/150‑1 (MQ=255)
                                                                                                                                        agaataacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGt              >  1:1602918/8‑150 (MQ=255)
                                                                                                                                        agaataacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGt              >  2:433964/8‑150 (MQ=255)
                                                                                                                                         gaataacTGGCCCGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGCCCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGtt             >  2:1204657/7‑150 (MQ=255)
                                                                                                                                          aataacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTc            >  1:1489105/6‑150 (MQ=255)
                                                                                                                                           ggTTCTGGGAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTcc           <  1:1292835/149‑1 (MQ=255)
                                                                                                                                           cGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTcc           <  1:1230678/150‑1 (MQ=255)
                                                                                                                                           ataacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTcc           >  1:1408473/5‑150 (MQ=255)
                                                                                                                                            taacTGGCACGCCTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCTGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCt          >  2:728859/4‑150 (MQ=255)
                                                                                                                                             ttCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTg         <  2:764111/150‑1 (MQ=255)
                                                                                                                                                tGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGaaa      >  2:458793/1‑150 (MQ=255)
                                                                                                                                                tGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGaaa      >  2:757570/1‑150 (MQ=255)
                                                                                                                                                tGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGaaa      <  1:1530961/150‑1 (MQ=255)
                                                                                                                                                tGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATAATTAACCCAGGGCGAAGATGCGGCGTTCCTGaaa      >  1:370468/1‑150 (MQ=255)
                                                                                                                                                 ggCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAg     >  1:1342970/1‑150 (MQ=255)
                                                                                                                                                   cAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGaa   >  1:1107838/1‑150 (MQ=255)
                                                                                                                                                    aGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAAc  <  2:103190/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAAGCCTCGCAGAAACACATTTTGCGTTACGCCGAACTGGCGGAAACGATGGCGGCAAACTGCACAGATGCCCAGCGTCGCGAAGAGCTGCTGACTATTGCAGAGATCTCCCGCCATAACGCGCAACATAAGCCGCAGACGTTCTGGCAGGCGTGCCAGTTATTCTGGTACATGAACATCATTCTGCAATACGAATCCAACGCCAGTTCGCTATCGTTGGGGCGCTTCGACCAGTATATGTTGCCGTTCTATCAGACATCATTAACCCAGGGCGAAGATGCGGCGTTCCTGAAAGAAC  >  NC_000913/4144583‑4144880

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: