Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,366,223 A→G F184S (TTC→TCC)  dcuA ← C4‑dicarboxylate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,366,2230AG97.5% 247.9 / ‑3.9 79F184S (TTC→TCC) dcuAC4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (42/35);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGAAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGAT  >  NC_000913/4366076‑4366369
                                                                                                                                                    |                                                                                                                                                  
cccAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCagggag                                                                                                                                                  <  1:417252/150‑1 (MQ=255)
 ccAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCagggagg                                                                                                                                                 >  2:1224788/1‑150 (MQ=255)
         aaCAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCa                                                                                                                                        <  1:784231/151‑1 (MQ=255)
          aCAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCAcc                                                                                                                                         >  1:1076132/1‑149 (MQ=255)
           cAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGa                                                                                                                                      <  1:535956/151‑1 (MQ=255)
            aGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGa                                                                                                                                      >  1:413552/1‑150 (MQ=255)
             gCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAAc                                                                                                                                    <  1:575923/151‑1 (MQ=255)
              ccAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAAcc                                                                                                                                   <  1:35175/151‑1 (MQ=255)
               cAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAGCCg                                                                                                                                  <  2:1266092/151‑1 (MQ=255)
                  acggacgTTTTGGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTCCAACCAGGCCCCCTC‑CCAGACGCTGGCGATAAATCGGAGCGTCAGAGAGTTTGGAGTGGAACAGCAAAGTAACCAGGGAGGACATCACCAAAACCGCca                                                                                                                               <  2:1228358/151‑1 (MQ=255)
                    ggacgTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCcagc                                                                                                                             >  1:625511/1‑151 (MQ=255)
                       cgTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTC‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCCCCAGAACCGCCAGCAgg                                                                                                                          <  1:236591/151‑1 (MQ=255)
                          ttttGCACCGGATTAGATTTCACTCTGCTTTTCACCGCGCAGTT‑AACCGGGCCCTTTTACAACCCGCTTGCGATAAAACGACTCGTGAGAGAGTTTGGAGTTCACCAGAATTGGGCCCAGGGAGGAAATCACCAGAACCGCCAGCAGGGt                                                                                                                        <  2:1155824/150‑1 (MQ=255)
                           tttGCACCGGATTTGATTTCATTCTGCTTTTCACCGCGCTGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGTATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGt                                                                                                                        <  2:925707/149‑1 (MQ=255)
                            ttcccccGGTTTTTATTTCAACCTGTTTTTCACCGGGCAGTCCAACCAGCCCCTCTT‑CCAGACGCTTGCGATAAATCGAATCGACAGAGAGTTGGGATTTGAACAGCATAGTGACCGGGGAGGACATCACCAGAACCGCCAGCAGGGTgg                                                                                                                      <  2:1436090/145‑1 (MQ=255)
                             tGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGa                                                                                                                     >  2:797080/1‑150 (MQ=255)
                             tGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGAc                                                                                                                    >  1:599548/1‑151 (MQ=255)
                              gCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACg                                                                                                                   >  1:885776/1‑151 (MQ=255)
                                   ggATTTGATTTCAATTTGTTTTTCCCCGCGCAGTTCAACCAGCCCCTCTT‑CCAGACGCTTGCGATAAATGGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGAACATCACCAGAACCGCCAGCAGGGTGGACGGGAt                                                                                                               <  1:746884/150‑1 (MQ=255)
                                   ggATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATg                                                                                                              >  2:230226/1‑151 (MQ=255)
                                   ggATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTGGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGAt                                                                                                               <  2:470393/150‑1 (MQ=255)
                                   ggATTTATTTTCAATCTGCTTTTCCCCGCGCAGTTCACCCAGGCCCTCTC‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTGGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATg                                                                                                              <  1:905890/151‑1 (MQ=255)
                                          atttCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACgg                                                                                                       >  1:190881/1‑151 (MQ=255)
                                           tttCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACgg                                                                                                       >  2:1066731/1‑150 (MQ=255)
                                           tttCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACgg                                                                                                       >  2:581424/1‑150 (MQ=255)
                                             tCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGAAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACgg                                                                                                       >  2:1450380/1‑148 (MQ=255)
                                              cAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGgag                                                                                                     >  1:401793/1‑149 (MQ=255)
                                                aTCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGc                                                                                                  >  1:44761/1‑150 (MQ=255)
                                                 tCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGc                                                                                                  >  2:876983/1‑149 (MQ=255)
                                                  cTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGc                                                                                                  >  1:616926/1‑148 (MQ=255)
                                                  cTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGc                                                                                                  >  2:1308767/1‑148 (MQ=255)
                                                       tttCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATgg                                                                                            >  1:578366/1‑149 (MQ=255)
                                                        ttCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAgg                                                                                         >  1:544472/1‑151 (MQ=255)
                                                           aCCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGt                                                                                        >  2:128602/1‑149 (MQ=255)
                                                                gcAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGAt                                                                                 >  1:763103/1‑151 (MQ=255)
                                                                acAGTTCACCCAGGCCTCCTG‑CCAGGGGCTTGCAATAAGTCGATTCGTGAGAGAGTTAGGGGTTTAACATCATAGTGGCCAGGGAGGACATACCCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAAAGCAGATGGAGGTAGCTg                                                                                   <  1:1339492/148‑1 (MQ=255)
                                                                 cAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGa                                                                                  <  2:1014681/149‑1 (MQ=255)
                                                                  aGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGc                                                                               >  2:1269554/1‑151 (MQ=255)
                                                                   gTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGcc                                                                              >  1:480381/1‑151 (MQ=255)
                                                                    ttCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGcc                                                                              >  2:218259/1‑150 (MQ=255)
                                                                       aaCCGGGCCCTCTT‑CCAGACGCTTGCGATAAATGGAATCGTAAGAGAGTTTGGATTTAAAAGGCATGGTGCCCGGGGAGGAAATCTCGAAAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCa                                                                             <  1:862780/148‑1 (MQ=255)
                                                                       aaCCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGGGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGCAGCTGATGCCa                                                                             >  1:204346/1‑148 (MQ=255)
                                                                       aaCCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATg                                                                           <  2:175165/150‑1 (MQ=255)
                                                                               ccTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCAt                                                                  >  1:688686/1‑151 (MQ=255)
                                                                               ccTCTC‑CCGGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGGGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCAt                                                                  <  1:1118233/151‑1 (MQ=255)
                                                                                ctctT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATc                                                                 >  1:798991/1‑151 (MQ=255)
                                                                                      ccAGACGCTTGCGATAAATCGGTGCGCCGGAGAGTTTGGATTGGACAAGCTTATTGACCGGGAAGGACACCACCAAAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGa                                                            <  1:672303/151‑1 (MQ=255)
                                                                                      ccAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACgg                                                             >  2:528674/1‑150 (MQ=255)
                                                                                      ccAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACgg                                                             >  2:591052/1‑150 (MQ=255)
                                                                                       cAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACgg                                                             <  2:851131/149‑1 (MQ=255)
                                                                                        aGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGa                                                            >  1:1177710/1‑149 (MQ=255)
                                                                                        aGACGCTGGCGATAAATCGGATCGTCAGAGAGTTGGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGaa                                                           <  2:175272/150‑1 (MQ=255)
                                                                                          aCGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGAc                                                        >  2:1036885/1‑151 (MQ=255)
                                                                                           cGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGAc                                                        <  1:1469465/150‑1 (MQ=255)
                                                                                            gCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACAt                                                      >  2:610524/1‑151 (MQ=255)
                                                                                            gCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACAt                                                      >  2:204831/1‑151 (MQ=255)
                                                                                            gCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACAt                                                      >  2:1340979/1‑151 (MQ=255)
                                                                                               tGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGt                                                    >  1:1241084/1‑150 (MQ=255)
                                                                                                gCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTa                                                   >  2:261529/1‑150 (MQ=255)
                                                                                                     aaaTCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTaaa                                                 <  1:950913/147‑1 (MQ=255)
                                                                                                     aaaTCGGATCGTAAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCGCAAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGt                                                    <  2:64179/144‑1 (MQ=255)
                                                                                                      aaTCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACt                                            >  1:1083284/1‑151 (MQ=255)
                                                                                                       aTCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTg                                           >  1:573358/1‑151 (MQ=255)
                                                                                                              cGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTg                                     >  1:576991/1‑150 (MQ=255)
                                                                                                              cGTCAGAGAGTTGGGGGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCt                                      <  2:885688/149‑1 (MQ=255)
                                                                                                               gTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTg                                     >  2:374845/1‑149 (MQ=255)
                                                                                                                  agagagTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGAtt                                <  2:52119/151‑1 (MQ=255)
                                                                                                                  agagagTTTGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGAt                                 <  1:372050/150‑1 (MQ=255)
                                                                                                                    agagTTTGGAGTTGAACAGCATGGTGACCAGGGAGAACATCCCAAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGAAAGACATGTAAACCACTGCCGCTga                                    <  1:1087312/145‑1 (MQ=255)
                                                                                                                        tttGGAGTTGAACAGCATAGTGACCAGGGAGGAAATCCCCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTg                               <  1:501131/146‑1 (MQ=255)
                                                                                                                         ttGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGAt                          <  1:186007/150‑1 (MQ=255)
                                                                                                                          tGGAGTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGc                        >  2:69141/1‑151 (MQ=255)
                                                                                                                           ggggTTGAACAGCATGGTGCCCGGGGAGGACACCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGc                        <  2:693483/147‑1 (MQ=255)
                                                                                                                              gTTGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTg                    >  1:1440553/1‑151 (MQ=255)
                                                                                                                               ttGAACAGCATGGTGCCCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGt                     <  2:926457/149‑1 (MQ=255)
                                                                                                                                tGAACAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGAt                  >  2:190959/1‑151 (MQ=255)
                                                                                                                                  aaCAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGAt                  <  1:230249/149‑1 (MQ=255)
                                                                                                                                  aaCAGCATAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGAt                  <  2:1329134/149‑1 (MQ=255)
                                                                                                                                        aTAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCt          <  2:554660/151‑1 (MQ=255)
                                                                                                                                        aTAGTGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCt          >  1:985968/1‑151 (MQ=255)
                                                                                                                                            tGACCAGGGAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCt          <  2:255124/147‑1 (MQ=255)
                                                                                                                                                 agggaggaCATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGAt  >  1:757370/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTT‑CCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGAAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATGCGGTGATCGCAATCTGCGCGGAT  >  NC_000913/4366076‑4366369

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑