Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,627,717 G→T 5.2% S67S (TCG→TCT ydfG → 3‑hydroxy acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,627,7170GT 5.2% 253.8 / 2.1 76S67S (TCG→TCTydfG3‑hydroxy acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (33/38);  new base T (2/2);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATGACGCGCGC  >  NC_000913/1627568‑1627863
                                                                                                                                                     |                                                                                                                                                  
cgtcgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCg                                                                                                                                                    <  1:564732/150‑1 (MQ=255)
  tcgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCt                                                                                                                                                  <  2:489726/150‑1 (MQ=255)
   cgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGCCAGCATCTCtt                                                                                                                                                 <  2:1089086/150‑1 (MQ=255)
   cgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCtt                                                                                                                                                 <  1:357455/150‑1 (MQ=255)
    gtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTc                                                                                                                                                >  2:1266111/1‑150 (MQ=255)
     ttttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTcc                                                                                                                                               <  1:736631/150‑1 (MQ=255)
     ttttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTcc                                                                                                                                               >  1:322964/1‑150 (MQ=255)
       ttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTg                                                                                                                                             <  1:1266111/150‑1 (MQ=255)
          ttCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGCCAGCATCTCTTcaatagc                                                                                                                                          <  2:1302262/150‑7 (MQ=255)
            caacaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                                                                                                                                        >  1:993247/1‑150 (MQ=255)
            caacaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                                                                                                                                        >  2:455338/1‑150 (MQ=255)
                aaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtggt                                                                                                                                    >  2:1376199/1‑150 (MQ=255)
                aaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtggt                                                                                                                                    <  1:963744/150‑1 (MQ=255)
                  gggCATAAAGTTATCGCCACTGGCCGTCGAAAGGAACTGTAGCAGGAGTTAAAAGACCAACTGGGTGATAATCTTTATATCGCCCAACTGGCCCTTCGCAACCGCGCCGCTATTCAAGAGATGATTTCATCCATTCCTTCCTAgttgttc                                                                                                                                  >  2:620167/1‑148 (MQ=255)
                      aTAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAata                                                                                                                              >  1:35653/1‑150 (MQ=255)
                        aaaGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATAtt                                                                                                                            <  2:828710/150‑1 (MQ=255)
                        aaaGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTACAATAtt                                                                                                                            <  2:743850/150‑1 (MQ=255)
                         aaGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTg                                                                                                                           <  2:98363/150‑1 (MQ=255)
                          aGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGa                                                                                                                          <  1:78607/150‑1 (MQ=255)
                          aGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGa                                                                                                                          <  1:1404592/150‑1 (MQ=255)
                            ttATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGata                                                                                                                        <  1:455338/150‑1 (MQ=255)
                            ttATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGata                                                                                                                        <  2:143353/150‑1 (MQ=255)
                            ttATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGata                                                                                                                        <  2:281109/150‑1 (MQ=255)
                            ttATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGata                                                                                                                        <  2:286276/150‑1 (MQ=255)
                             tATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGatat                                                                                                                       >  1:600161/1‑150 (MQ=255)
                                 gCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGCGCAATATTGATATCCTg                                                                                                                   >  1:672546/1‑150 (MQ=255)
                                   cACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGt                                                                                                                 >  1:279884/1‑150 (MQ=255)
                                       ggCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATAtt                                                                                                                            >  1:578130/1‑135 (MQ=255)
                                       ggCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATAtt                                                                                                                            <  2:578130/135‑1 (MQ=255)
                                       ggCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATGTCCTGGTAaat                                                                                                             >  1:75948/1‑150 (MQ=255)
                                       ggCCGTCGCCAGGAACGGTTGCAGGAATTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAaat                                                                                                             >  1:1057361/1‑150 (MQ=255)
                                               ccAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtg                                                                                                                                      <  1:391131/117‑1 (MQ=255)
                                               ccAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtg                                                                                                                                      >  2:391131/1‑117 (MQ=255)
                                                        gTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGtt                                                                                            <  2:956254/150‑1 (MQ=255)
                                                         ttGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTg                                                                                           <  2:322964/150‑1 (MQ=255)
                                                          tGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTgg                                                                                          <  1:917560/150‑1 (MQ=255)
                                                            cAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTTCAATATTGATATCCTTGTATATAATGCCGGCCTGGCGTTgggg                                                                                        <  2:739587/150‑2 (MQ=255)
                                                                aGTTAAAAGACTAAATTGGAGATAATCTTTATATCGCCCAACTGGAAGTTAGCAACGGCGCAGCTATTGAAGAGATGCTGGCATCGCTTCCTGGCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGTCCTGGCGTTGGGCATgg                                                                                    <  1:620167/150‑1 (MQ=255)
                                                                   tAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGc                                                                                 >  1:533604/1‑150 (MQ=255)
                                                                       aGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTgc                                                                             >  2:1355700/1‑150 (MQ=255)
                                                                          cGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCa                                                                          >  2:1257819/1‑150 (MQ=255)
                                                                            aaCTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATa                                                                        <  1:382534/150‑1 (MQ=255)
                                                                                gggAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGc                                                                    <  1:1221882/150‑1 (MQ=255)
                                                                                   agaTAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAg                                                                 <  1:880562/150‑1 (MQ=255)
                                                                                      tAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGt                                                              >  1:545384/1‑150 (MQ=255)
                                                                                        aTCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTg                                                            >  1:185338/1‑150 (MQ=255)
                                                                                        aTCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTg                                                            <  2:1051857/150‑1 (MQ=255)
                                                                                           tGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCGTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAg                                                         <  2:985702/150‑1 (MQ=255)
                                                                                           tGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAg                                                         >  1:794901/1‑150 (MQ=255)
                                                                                           tGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAg                                                         >  2:102015/1‑150 (MQ=255)
                                                                                               tatCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTg                                                     >  1:154810/1‑150 (MQ=255)
                                                                                               tatCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTg                                                     >  1:186292/1‑150 (MQ=255)
                                                                                                atCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTgg                                                    <  2:545384/150‑1 (MQ=255)
                                                                                                  cGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGa                                                  >  1:992618/1‑150 (MQ=255)
                                                                                                    cccAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGaaa                                                >  2:1281345/1‑150 (MQ=255)
                                                                                                     ccAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAAc                                               >  1:888615/1‑150 (MQ=255)
                                                                                                        aCTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACgat                                            <  2:228004/150‑1 (MQ=255)
                                                                                                            gACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGAtgatt                                        <  2:785398/150‑1 (MQ=255)
                                                                                                               gTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGTAAACGAtgattgat                                     <  2:185338/150‑1 (MQ=255)
                                                                                                                 tCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATAc                                   <  1:102015/150‑1 (MQ=255)
                                                                                                                  cGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATAcc                                  <  1:1376199/150‑1 (MQ=255)
                                                                                                                     aaCCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCAGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACcaaa                               >  2:928386/1‑149 (MQ=255)
                                                                                                                         gcgcCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACcaacaaca                           <  2:993247/150‑1 (MQ=255)
                                                                                                                            ccgcTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAg                        >  1:762560/1‑150 (MQ=255)
                                                                                                                             cgcTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAgg                       <  1:277942/150‑1 (MQ=255)
                                                                                                                                   tGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGt                 >  1:994594/1‑150 (MQ=255)
                                                                                                                                    gAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGta                >  2:474328/1‑150 (MQ=255)
                                                                                                                                     aaGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGtat               >  1:1243542/1‑150 (MQ=255)
                                                                                                                                      agagATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGtata              <  2:206891/150‑1 (MQ=255)
                                                                                                                                      agagATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACAAACAACAAAGGCCTGGtata              >  1:978498/1‑150 (MQ=255)
                                                                                                                                       gagaTGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATGTCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGtatat             <  2:75948/150‑1 (MQ=255)
                                                                                                                                        cgaTGCCAGCATCTCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGtatata            >  1:1086774/2‑150 (MQ=255)
                                                                                                                                         gaTGCCAGCATCTCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTatatat           >  2:700636/1‑150 (MQ=255)
                                                                                                                                          aTGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGTTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCTTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATg          <  2:35653/150‑1 (MQ=255)
                                                                                                                                           tGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATGa         >  1:589685/1‑150 (MQ=255)
                                                                                                                                               ggCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATGAcgcg     <  1:925649/150‑1 (MQ=255)
                                                                                                                                                  aTCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATGACgcgcgc  >  1:557582/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATGACGCGCGC  >  NC_000913/1627568‑1627863

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: