Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,873 T→A 20.2% intergenic (+90/+62) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,504,8730TA20.2% 164.8 / 5.1 71intergenic (+90/+62)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base T (10/42);  new base A (12/2);  total (27/44)
Fisher's exact test for biased strand distribution p-value = 8.19e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCT  >  NC_000913/3504738‑3505022
                                                                                                                                       |                                                                                                                                                     
ccGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCTAACCCTCTCCCCACAGGGGc                                                                                                                                         >  2:1521996/1‑150 (MQ=255)
 cgcgCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCg                                                                                                                                        <  2:271812/150‑1 (MQ=255)
 cgcgCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCg                                                                                                                                        <  2:271762/150‑1 (MQ=255)
  gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGa                                                                                                                                       <  1:34775/150‑1 (MQ=255)
  gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGa                                                                                                                                       >  2:237091/1‑150 (MQ=255)
   cgcCAATCGCAACACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                                                                                                                      <  2:817412/150‑1 (MQ=255)
    gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTACCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                                     <  1:221553/150‑1 (MQ=255)
    gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgg                                                                                                                                     <  1:684393/150‑1 (MQ=39)
    gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgg                                                                                                                                     <  1:1466277/150‑1 (MQ=39)
    gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCGCCCCTTTGGGGAGAgg                                                                                                                                     <  1:1360859/150‑1 (MQ=39)
     ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                                    <  1:1090931/150‑1 (MQ=255)
       aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                                                  <  2:969267/150‑1 (MQ=255)
              tCACCTTAACTATTGATAATAAAAGACATTAACCGTCGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTTAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCACTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                                           <  1:275772/150‑1 (MQ=255)
              tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                                           <  1:1287503/150‑1 (MQ=255)
                   ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCACTCCCCAAAGGGGCGAGGGGACCGTTCAGTACg                                                                                                                      >  1:260982/1‑150 (MQ=255)
                    tAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACgg                                                                                                                     <  1:1131511/150‑1 (MQ=255)
                      aCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGaa                                                                                                                   <  2:74924/150‑1 (MQ=255)
                        tATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAACCCACACCCCAAACGGGCGAGGGGACCGTTCAGTACGGAAtt                                                                                                                 >  1:1036370/1‑150 (MQ=39)
                         aTTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCGCACCCCCAAGGGGCGGGGGGGCCGGTCAGTACGGAAttt                                                                                                                >  1:1404482/1‑150 (MQ=38)
                           tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                                              <  1:816788/150‑1 (MQ=255)
                           tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTCACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                                              <  1:816783/150‑1 (MQ=255)
                            gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTAAACCAAAAACAAAAAACAAATCGGCGAGGGGACCCTTCATTAAGGAAATTAGc                                                                                                             >  1:1043218/1‑150 (MQ=21)
                            gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGCGGGGCCCGTTAAGTACGGAAttttgg                                                                                                             >  1:933335/1‑147 (MQ=39)
                             ataataAAAGACATTAACCGTCGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTACCCCAAAACCGCTCCCCCCAGGGGCGCGGGGGCCGTGCAGGACGGACATTAGCt                                                                                                            >  2:1609574/1‑150 (MQ=37)
                                ataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGa                                                                                                         <  2:874730/150‑1 (MQ=255)
                                 taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAc                                                                                                        <  1:432502/150‑1 (MQ=255)
                                  aaaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAAACCACAACCCCAAGGGGCGAGGGGACCGTTCAGTAAGGAATTTAGCTTgtag                                                                                                       >  1:1032665/1‑147 (MQ=38)
                                       aCATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTGGAATTTTACATCCCGTACGTTCCCCGGACCCAAACCCACTCCCCCAAGGGGCGAGGGGACCGGTCAGTTCCGCATTTAGCTTGGCGGTGCt                                                                                                  >  2:1499931/1‑150 (MQ=37)
                                          ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATc                                                                                               <  2:300823/150‑1 (MQ=255)
                                            aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                             <  1:805251/150‑1 (MQ=255)
                                               cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                          <  1:1365202/150‑1 (MQ=255)
                                                    gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaaa                                                                                     <  1:188167/150‑1 (MQ=255)
                                                    gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaaa                                                                                     <  2:760326/150‑1 (MQ=255)
                                                        gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAAGCCTCTCCCCAAAGGGGCGTGGGGTGCGTTCAGTACGGAATTTAGCTTGTAGGTGCTATACAGCATTAAATACt                                                                                 >  2:113465/1‑150 (MQ=255)
                                                        gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACt                                                                                 <  1:925650/150‑1 (MQ=255)
                                                         aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                <  2:1569697/150‑1 (MQ=255)
                                                         aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                <  1:1455167/150‑1 (MQ=255)
                                                          tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                               <  2:464709/150‑1 (MQ=255)
                                                          tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACAACCAAAACCAAAAACCAAAGGGGCGAGTGGAAAGATAAGTACGGAAATTATATAGACGGTTCTATACCTCATTAAATACAgg                                                                               >  2:204456/1‑150 (MQ=255)
                                                          aaagaaccgcCACCATCGCAACCGCCGATGCAGGCGCGAAGAAATTAACAACCCCTACCTTCCCCACACCCAAACCCACTCCCCAAACCGGCGAGCGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCACCATTAAATACTgg                                                                               <  2:1703762/143‑1 (MQ=255)
                                                                gcgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCGAACCCACTCCCCAAAGGGGCGAGGGGGCCGTTCAGTACGGAAGTTAGCATGGCGGTGCTATCCACCAATCAATACTGGccaaac                                                                         >  2:690959/1‑145 (MQ=255)
                                                                   ccAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCCCAACCAAACCCCCAACCCAAAAGGGCGGGGGGGCCCGTCACTAAGGGAATTATCTTGACGGGGCTAAACAGCATTAAAAACATGCaaaaaaaa                                                                      >  2:301602/1‑144 (MQ=255)
                                                                    cAGCATCGCCTCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACtctc                                                                     <  2:1486855/150‑1 (MQ=255)
                                                                    cAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACtctc                                                                     >  2:429401/1‑150 (MQ=255)
                                                                    cAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTGCCCCTCACCCGAACCCCCTCCCCCAACGGGCGAGGGGTCCGGTCACTAAGGTAATTATCGTTACGGGGCGATACACCATTTAATAAAGGCAATAAtctc                                                                     >  2:806224/1‑150 (MQ=255)
                                                                      gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTACCCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                   <  2:914989/150‑1 (MQ=255)
                                                                      gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                   <  2:914984/150‑1 (MQ=255)
                                                                                 ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAAACCACTACCCAAAGGGGCGGGGGGAACGTTCACTACGGAAATTAGCATGACGGTGCTATACAGCAGAAAATACAGGCAATACTCTCAACCAAAATCCgc                                                        >  2:1204897/1‑150 (MQ=255)
                                                                                 ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                        <  1:12455/150‑1 (MQ=255)
                                                                                      tGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAAAACTGGCAATACTCTCACGCAAAAACCGcaccac                                                   >  2:949797/1‑150 (MQ=255)
                                                                                         tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGt                                                <  2:1474176/150‑1 (MQ=255)
                                                                                         tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGt                                                <  2:790536/150‑1 (MQ=255)
                                                                                         tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGt                                                <  2:1212056/150‑1 (MQ=255)
                                                                                         tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGt                                                <  2:1150620/150‑1 (MQ=255)
                                                                                          ggCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCAAAACCACAACCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAATTATCTTGAAGGTGGTATACAGCAATAAATAATGTCAAAAATATCACGCCAAATACGCAGCCCCGtg                                               >  2:459719/1‑149 (MQ=255)
                                                                                          ggCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCGAAACCACACCCCAAAGGGGCGGGGGGACCGGTGAGTAAGGAAATTAGCTTTACGGTGCTATACAGCAATAAAAAATGGCAAAAATCTCACGCAAAAACCGCACCACCggc                                               >  2:499123/1‑148 (MQ=255)
                                                                                                 tGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCACTCCCCAAAGGGGCGAGGGGACCGTTCAGTAAGGAATTTAGCTTGACGGTGCTATACAGCAATAAAAAATGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTAAc                                        >  2:266565/1‑150 (MQ=38)
                                                                                                      tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                   <  2:101966/150‑1 (MQ=255)
                                                                                                      tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                   <  2:887516/150‑1 (MQ=255)
                                                                                                            tCCCGTACGTGCCCCACACCCGAACCCACACCCCAAAGGGGCGAGGGGAACGGTCAGTACGGAATTTATCTTGACGGTGCTATCCAGCATTAAATACAGGCAATACTAACACCCAAAATACGCAGCACCGTCTCGTCAACGCTGCGATTc                             >  2:432184/1‑150 (MQ=255)
                                                                                                              ccGTACGTTCCCCTCCCCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                           <  2:1043218/150‑1 (MQ=255)
                                                                                                              ccGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                           <  2:933335/150‑1 (MQ=255)
                                                                                                              ccGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                           <  2:194928/150‑1 (MQ=255)
                                                                                                              ccGTACGTTCCCCACACCCGAACCCGCACCCCAAAGGGGCGAGGGGCCCCGTCAGTAAGGAAATTAGCTTGACGGTGCTATCCACCATTAAATAATGGCAAAACACTAACGCAAAAACCGCAGCACCGTATCATAACCGCAGCGATTCgg                           >  2:88596/1‑150 (MQ=255)
                                                                                                                      tCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGc                   <  1:1499931/150‑1 (MQ=255)
                                                                                                                           tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc              <  1:1385349/150‑1 (MQ=255)
                                                                                                                             aCCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGc            >  1:148396/1‑150 (MQ=255)
                                                                                                                               ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGACTCGcac          >  1:308346/1‑150 (MQ=255)
                                                                                                                               ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac          >  1:405390/1‑150 (MQ=255)
                                                                                                                               ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTGCAGTACGGACTTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCACAATCCGCCGCTCCGTCTCTTCGCCGCTGCGATTCGGGTTAATGCGACTCGcac          >  2:1689344/1‑150 (MQ=255)
                                                                                                                                  caCCCTCGCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTAGCCAGCATTAACTACTGGCAATACTCTCACGCAAACTCCGCAGCCCCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACTAt       <  2:161256/149‑1 (MQ=255)
                                                                                                                                       tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt  >  2:416965/1‑150 (MQ=255)
                                                                                                                                       |                                                                                                                                                     
CCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCT  >  NC_000913/3504738‑3505022

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