Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,448,678 C→T 6.3% P387L (CCA→CTA)  feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,448,6780CT 6.3% 182.0 / 2.2 63P387L (CCA→CTA) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base T (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGATGTCACCTGTTGCACAGATTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGG  >  NC_000913/1448529‑1448826
                                                                                                                                                     |                                                                                                                                                    
ggggATGTCACCTGTTGCACAGATTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGcc                                                                                                                                                      >  1:351421/1‑150 (MQ=255)
ggggATGTCACCTGTTGCACAGATTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGcc                                                                                                                                                      <  2:1263982/150‑1 (MQ=255)
     tGTCACCTGTTGCACAGATTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGc                                                                                                                                                 >  1:806603/1‑150 (MQ=255)
      gTCACCTGTTGCACAGATTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCTAACGCt                                                                                                                                                >  2:902671/1‑150 (MQ=255)
      gTCACCTGTTGCACAGATTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCTAACGCt                                                                                                                                                >  2:1313474/1‑150 (MQ=255)
           cTGTTGCACAGATTAACCCTGTGGTTTCTCGAGCGCACTGCGACAAAGTGTGTGCATTCCTCGACGATGCGCAGGCACAGCATGCAGTGCTGATTGGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCtggtgg                                                                                                                                           <  1:1030533/150‑1 (MQ=255)
                gCACAGATTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAAt                                                                                                                                      >  1:425669/1‑150 (MQ=255)
                gCACAGATTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAAt                                                                                                                                      >  1:375971/1‑150 (MQ=255)
                 cacaGATTACCCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGGTCATTCCTCGACTATGCGCAGGCAGATCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTTCGCCAACGCTGGTGGTAAATc                                                                                                                                     <  2:943477/150‑1 (MQ=255)
                 cacaGATTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATc                                                                                                                                     >  2:438793/1‑150 (MQ=255)
                   caGATTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCACCAACGCTGGTGGTAAATccc                                                                                                                                   >  2:11477/1‑150 (MQ=255)
                      ttttaCCCTTTGGTTTCTCGTGCGCAGTGCGACAAAGGGTGGTCATTCCTCGACGATCCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGAt                                                                                                                                <  2:490379/146‑1 (MQ=255)
                      aTTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCgct                                                                                                                                >  1:793527/1‑148 (MQ=255)
                          aCCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTa                                                                                                                            >  1:1278762/1‑150 (MQ=255)
                          aCCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGGGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTa                                                                                                                            <  1:1083290/150‑1 (MQ=255)
                           cccTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTaa                                                                                                                           >  2:1242720/1‑150 (MQ=255)
                             cTTTGGTTTCTCGTGCGCAAGGCGTAAAATTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAACCAGAGCTGATTCGCGGGGCGAATTGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAAt                                                                                                                         <  2:345292/150‑1 (MQ=255)
                              tttGGTTTCTCGTGCGCACTGCTACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCACGCAGAGCTGATTCGCGGGTCGAAAGGACCAGCCGGAGAGGGGTCTTATGTTGCGCCAACGCTGGTGGTATATCCCGATGCTAAAtt                                                                                                                        <  2:165574/150‑1 (MQ=255)
                                 ggTTTCTCGTGCGCACTGCGACAAAGTGTGTTCCTTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCGAAATTgcg                                                                                                                     >  1:835394/1‑150 (MQ=255)
                                  gTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTgcgc                                                                                                                    <  2:1278762/150‑1 (MQ=255)
                                   tttCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCt                                                                                                                   >  1:1021200/1‑150 (MQ=255)
                                          tGCGCACTGCGACAAAGTGTGTTGATTACTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCg                                                                                                            >  1:1088805/1‑150 (MQ=255)
                                          ggcgcACTGCGTCACAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCTGTGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTTTTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAATTTGCGCTTAACTCg                                                                                                            <  1:141286/149‑1 (MQ=255)
                                            cgcACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTg                                                                                                          <  2:626899/150‑1 (MQ=255)
                                            cgcACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTg                                                                                                          <  1:1046899/150‑1 (MQ=255)
                                                       aaaGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGttt                                                                                               <  1:1085612/150‑1 (MQ=255)
                                                              gtTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGCTTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCgg                                                                                        >  2:559081/1‑150 (MQ=255)
                                                                 cATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCggtgg                                                                                     <  2:134196/150‑1 (MQ=255)
                                                                   ttCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTACGGTGGTa                                                                                   >  2:151998/1‑150 (MQ=255)
                                                                    accTCGACGATGCGCAGGCACAGCAAGGAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTaa                                                                                  <  2:474042/149‑1 (MQ=255)
                                                                     ccTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTaaa                                                                                 >  2:657357/1‑150 (MQ=255)
                                                                       ccgccgaTGCGCAGTCACCGCAAGCCGAGCTGATTCGCGGGCCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTTGTAAATCCCGATGCTAACTTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAAcc                                                                               <  1:127796/146‑1 (MQ=255)
                                                                          acgaTGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTgg                                                                            <  2:1127443/150‑1 (MQ=255)
                                                                             aTGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGc                                                                         <  2:953041/150‑1 (MQ=255)
                                                                              tGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCg                                                                        >  2:599215/1‑150 (MQ=255)
                                                                                  cAGGCACAGCAAGCAGCGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTa                                                                    <  1:438793/150‑1 (MQ=255)
                                                                                        caGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGAt                                                              >  1:513276/1‑150 (MQ=255)
                                                                                          tcaagcaGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGGTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATgg                                                            <  1:897468/149‑1 (MQ=255)
                                                                                           caagcaGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGga                                                           <  1:436823/150‑1 (MQ=255)
                                                                                            aagcaGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGgag                                                          >  1:9313/1‑150 (MQ=255)
                                                                                            aagcaGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGgag                                                          >  1:1346675/1‑150 (MQ=255)
                                                                                                   gCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGc                                                   <  1:937541/150‑1 (MQ=255)
                                                                                                    cTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCg                                                  <  2:379726/150‑1 (MQ=255)
                                                                                                           gcgGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAAc                                           >  2:874962/1‑150 (MQ=255)
                                                                                                            cgGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACt                                          <  2:1002983/150‑1 (MQ=255)
                                                                                                             gggTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTg                                         >  1:23701/1‑150 (MQ=255)
                                                                                                             gggACGAATGGACCAGCCGGACAGGGCTATTATGTTGCGCCACCGCGGGTGGTAAAACCCGATCCTAAAATGCGCTTAACTCCCGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTg                                         <  2:233661/150‑1 (MQ=255)
                                                                                                               gTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGc                                       <  1:758068/150‑1 (MQ=255)
                                                                                                                  gactGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCaaa                                    <  1:533889/147‑1 (MQ=255)
                                                                                                                  gAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCaaa                                    >  2:969590/1‑150 (MQ=255)
                                                                                                                   aaTGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAAc                                   <  1:657357/150‑1 (MQ=255)
                                                                                                                       gACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGaca                               >  2:1096706/1‑150 (MQ=255)
                                                                                                                        aCCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGacac                              >  1:436280/1‑150 (MQ=255)
                                                                                                                         ccAGCCGGCGAGGGGGATTATGTTGCCCCAACGCTGGTGGTAAATCCCGATGCTAAACTGCGCTAAACTCGTGAAGAGGGGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACg                             <  2:501733/150‑1 (MQ=255)
                                                                                                                             ccGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAat                         <  2:827386/150‑1 (MQ=255)
                                                                                                                                  gagGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGc                    >  1:1327987/1‑150 (MQ=255)
                                                                                                                                   tggtgTATTATGTTGCGCTAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCt                   <  1:1313474/146‑1 (MQ=255)
                                                                                                                                   agGGGTATTATGTTGCGCTAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCt                   <  1:902671/150‑1 (MQ=255)
                                                                                                                                    ggggTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCtt                  >  2:1189732/1‑150 (MQ=255)
                                                                                                                                     gggTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTa                 >  2:558134/1‑150 (MQ=255)
                                                                                                                                         attatGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGCCAAACGACACGGAATATGGCTTAACTg             >  1:1242708/1‑150 (MQ=255)
                                                                                                                                           tatGGTGCGCCAACGCTGGTGGTAAATCCCGATCCTAAATTGCGCTTAACTCGTGACGAGGTGTTTGGTCCGGTGGTAAACCGGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGcc           <  2:207774/150‑1 (MQ=255)
                                                                                                                                             tGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAg         >  1:76966/1‑150 (MQ=255)
                                                                                                                                               ttGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGGCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAgtg       <  2:617023/150‑1 (MQ=255)
                                                                                                                                                  cgcCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCt    <  2:1272543/150‑1 (MQ=255)
                                                                                                                                                    ccAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTgg  >  1:1150545/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGGATGTCACCTGTTGCACAGATTAACCCTTTGGTTTCTCGTGCGCACTGCGACAAAGTGTGTTCATTCCTCGACGATGCGCAGGCACAGCAAGCAGAGCTGATTCGCGGGTCGAATGGACCAGCCGGAGAGGGGTATTATGTTGCGCCAACGCTGGTGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGG  >  NC_000913/1448529‑1448826

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