Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,282,928 T→A 5.5% I318N (ATC→AAC) ‡ agaS → putative galactosamine‑6‑phosphate deaminase/isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,282,9280TA 5.5% 228.7 / 2.5 72I318N (ATC→AAC) ‡agaSputative galactosamine‑6‑phosphate deaminase/isomerase
Reads supporting (aligned to +/- strand):  ref base T (35/33);  new base A (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCCAGTG  >  NC_000913/3282780‑3283077
                                                                                                                                                    |                                                                                                                                                     
aaaTCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACAtc                                                                                                                                                      <  2:1474315/150‑1 (MQ=255)
     gCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGc                                                                                                                                                 >  1:1547011/1‑150 (MQ=255)
         gTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGct                                                                                                                                             <  1:843867/150‑2 (MQ=255)
          tCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGct                                                                                                                                             >  2:767216/1‑148 (MQ=255)
          tCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACGATGTCGCTGct                                                                                                                                             <  1:767216/149‑2 (MQ=255)
            gATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGACAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGt                                                                                                                                          >  1:714435/1‑150 (MQ=255)
             aTGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTc                                                                                                                                         <  1:716601/150‑1 (MQ=255)
              tGAAGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGACGGTcc                                                                                                                                        >  1:676440/1‑150 (MQ=255)
               gACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCca                                                                                                                                       <  1:1828746/150‑1 (MQ=255)
                      cGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATatca                                                                                                                                >  2:787142/1‑150 (MQ=255)
                      cGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATatca                                                                                                                                <  1:402549/150‑1 (MQ=255)
                       gCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATatcat                                                                                                                               <  1:1881902/150‑1 (MQ=255)
                         tggtggtggTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATcc                                                                                                                             <  2:175429/150‑1 (MQ=255)
                           gtggtggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTg                                                                                                                           <  1:1133140/150‑1 (MQ=255)
                            tggtggtATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGc                                                                                                                          <  2:1084374/150‑1 (MQ=255)
                                     ttGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcac                                                                                                                 >  2:1430354/1‑150 (MQ=255)
                                      tGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACggtcacg                                                                                                                <  2:329645/150‑1 (MQ=255)
                                      tGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacg                                                                                                                <  2:487168/150‑1 (MQ=255)
                                      tGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacg                                                                                                                >  2:42020/1‑150 (MQ=255)
                                      tGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGACAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacg                                                                                                                <  2:714435/150‑1 (MQ=255)
                                       gTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgt                                                                                                               >  2:547265/1‑150 (MQ=255)
                                        tctcCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGGGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgtc                                                                                                              <  2:779687/150‑1 (MQ=255)
                                            cAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACtt                                                                                                          <  1:1611128/150‑1 (MQ=255)
                                             aGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACttt                                                                                                         >  1:666851/1‑150 (MQ=255)
                                                cACCCTAACAACCGTCAGTATGATCTTGATCTGATGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCCCCATCCCCGCGGAAAGCAGCGACATAGTCGCTCCCGGTCCACATATCATACTGCCACCGTCACGTCACTTTATc                                                                                                      >  2:933238/1‑150 (MQ=255)
                                                     ttACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgtc                                                                                                              <  2:1202157/137‑1 (MQ=255)
                                                     ttACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACcgtcacgtc                                                                                                              >  1:1202157/1‑137 (MQ=255)
                                                      tACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATAATCCTGCCACCGTCACGTCACTTTATCGACGtt                                                                                                >  2:1427329/1‑150 (MQ=255)
                                                       acacCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTg                                                                                               <  1:547265/150‑1 (MQ=255)
                                                       acacCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTg                                                                                               <  1:1427329/150‑1 (MQ=255)
                                                       acacCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTg                                                                                               <  2:1183263/150‑1 (MQ=255)
                                                        cacCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGa                                                                                              <  1:976049/150‑1 (MQ=255)
                                                          cccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGc                                                                                            >  1:1863354/1‑150 (MQ=255)
                                                          cccGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGc                                                                                            >  2:746665/1‑150 (MQ=255)
                                                             gTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAgg                                                                                         >  1:888364/1‑150 (MQ=255)
                                                                    tgatcttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTg                                                                                  <  1:905041/150‑1 (MQ=255)
                                                                         ttgatctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTcc                                                                             >  1:568050/1‑150 (MQ=255)
                                                                            atctGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGa                                                                          >  1:1663664/1‑150 (MQ=255)
                                                                               tgctgGCTGAACTTCGCCGTGACAACCTGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGt                                                                       >  1:1573914/1‑150 (MQ=255)
                                                                                gctggctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTa                                                                      >  1:918546/1‑150 (MQ=255)
                                                                                 ctggctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTAc                                                                     <  1:1430354/150‑1 (MQ=255)
                                                                                     ctgAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGccc                                                                 <  2:1607286/150‑1 (MQ=255)
                                                                                          cTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACg                                                            >  1:1079999/1‑150 (MQ=255)
                                                                                          cTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACg                                                            >  1:1672229/1‑150 (MQ=255)
                                                                                          cTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACg                                                            >  1:1633810/1‑150 (MQ=255)
                                                                                            tCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGtt                                                          <  2:666851/150‑1 (MQ=255)
                                                                                                 gTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCAc                                                     <  1:787142/150‑1 (MQ=255)
                                                                                                  tGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACt                                                    <  1:1160088/150‑1 (MQ=255)
                                                                                                     cAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGAt                                                 >  2:814529/1‑150 (MQ=255)
                                                                                                     cAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGAt                                                 >  2:1754232/1‑150 (MQ=255)
                                                                                                        ccAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCa                                              <  1:933238/150‑1 (MQ=255)
                                                                                                         cAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAg                                             <  1:1619659/150‑1 (MQ=255)
                                                                                                          aGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGt                                            >  1:1024593/1‑150 (MQ=255)
                                                                                                          aGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGt                                            >  2:966596/1‑150 (MQ=255)
                                                                                                           ggCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTc                                           <  2:1079999/150‑1 (MQ=255)
                                                                                                              aaTGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCt                                        <  2:1547011/150‑1 (MQ=255)
                                                                                                                tGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGc                                      <  1:1202673/150‑1 (MQ=255)
                                                                                                                 gCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGca                                     <  2:1024593/150‑1 (MQ=255)
                                                                                                                     gtAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATg                                 >  2:345586/1‑150 (MQ=255)
                                                                                                                      tAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATgg                                >  2:130429/1‑150 (MQ=255)
                                                                                                                       aaTCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATggg                               <  1:42020/150‑1 (MQ=255)
                                                                                                                        atcgccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGc                              >  2:402849/1‑150 (MQ=255)
                                                                                                                           gccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACAt                                  >  2:382432/1‑143 (MQ=255)
                                                                                                                           gccatcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACAt                                  <  1:382432/143‑1 (MQ=255)
                                                                                                                             catcgccGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATAc                         >  2:269081/1‑150 (MQ=255)
                                                                                                                                  ccgcGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCgg                    >  1:952704/1‑150 (MQ=255)
                                                                                                                                  ccgcGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCgg                    >  2:946629/1‑150 (MQ=255)
                                                                                                                                    gcgGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGAt                  <  2:1213833/150‑1 (MQ=255)
                                                                                                                                          ggcagcGACGATGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCa            >  1:1177675/2‑150 (MQ=255)
                                                                                                                                              gcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAg        >  2:1005540/1‑150 (MQ=255)
                                                                                                                                              gcGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAg        >  1:887429/1‑150 (MQ=255)
                                                                                                                                                gACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGcc      <  2:1672229/150‑1 (MQ=255)
                                                                                                                                                  cATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCCAg    >  1:647241/1‑150 (MQ=255)
                                                                                                                                                    tcgtcgCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCCAGTg  <  2:139936/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAATCGCTGGTCGATGACGAAACGCTGGTGGTGGTATTTGTCTCCAGCCACCCTTACACCCGTCAGTATGATCTTGATCTGCTGGCTGAACTTCGCCGTGACAACCAGGCAATGCGTGTAATCGCCATCGCCGCGGAAAGCAGCGACATCGTCGCTGCCGGTCCACATATCATCCTGCCACCGTCACGTCACTTTATCGACGTTGAGCAGGCATTTTGCTTCCTGATGTACGCCCAGACGTTTGCACTGATGCAGTCGCTGCACATGGGCAATACGCCGGATACCCCATCAGCCAGTG  >  NC_000913/3282780‑3283077

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: