Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,627,711 G→A 5.3% L65L (CTG→CTA) ‡ ydfG → 3‑hydroxy acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,627,7110GA 5.3% 239.0 / 2.1 76L65L (CTG→CTA) ‡ydfG3‑hydroxy acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (33/38);  new base A (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 = 6.81e-01

CTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATGACGCG  >  NC_000913/1627566‑1627860
                                                                                                                                                 |                                                                                                                                                     
cTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCAt                                                                                                                                                   >  1:985702/1‑150 (MQ=255)
  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:322964/1‑150 (MQ=255)
       ttttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTcc                                                                                                                                            <  1:736631/150‑1 (MQ=255)
         ttATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTg                                                                                                                                          <  1:1266111/150‑1 (MQ=255)
            ttCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGCGATGCCAGCATCTCTTcaatagc                                                                                                                                       <  2:1302262/150‑7 (MQ=255)
              caacaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                                                                                                                                     >  2:455338/1‑150 (MQ=255)
              caacaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                                                                                                                                     >  1:993247/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:1404592/150‑1 (MQ=255)
                            aGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGa                                                                                                                       <  1:78607/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                                                                                                                                   >  2:391131/1‑117 (MQ=255)
                                                 ccAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtg                                                                                                                                   <  1:391131/117‑1 (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                                                         <  2:1051857/150‑1 (MQ=255)
                                                                                          aTCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTg                                                         >  1:185338/1‑150 (MQ=255)
                                                                                             tGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCGTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAg                                                      <  2:985702/150‑1 (MQ=255)
                                                                                             tGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAg                                                      >  2:102015/1‑150 (MQ=255)
                                                                                             tGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAg                                                      >  1:794901/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)
                                                                                                                                                 |                                                                                                                                                     
CTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGGTAAATAATGCCGGCCTGGCGTTGGGCATGGAGCCTGCGCATAAAGCCAGCGTTGAAGACTGGGAAACGATGATTGATACCAACAACAAAGGCCTGGTATATATGACGCG  >  NC_000913/1627566‑1627860

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: