Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,494,308 T→G Q33H (CAA→CAC)  idnD ← L‑idonate 5‑dehydrogenase, NAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,494,3080TG98.8% 335.6 / ‑4.0 83Q33H (CAA→CAC) idnDL‑idonate 5‑dehydrogenase, NAD‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (45/37);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GCCCTTCATGTAATTCTGATGAGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATTTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTT  >  NC_000913/4494160‑4494454
                                                                                                                                                    |                                                                                                                                                  
gCCCGTCATGTAATTCTGATGAGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTAtgtg                                                                                                                                                  <  1:207579/151‑1 (MQ=255)
gCCCGTCATGTAATTCTGATGAGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTAtgtg                                                                                                                                                  <  2:781922/151‑1 (MQ=255)
    gtCATGTAATTCTGATGAGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACt                                                                                                                                              >  2:1241231/2‑151 (MQ=255)
     tCATGTAATTCTGATGAGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACt                                                                                                                                              >  2:179267/1‑150 (MQ=255)
         gTAATTCTGATGAGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATg                                                                                                                                          >  2:1097188/1‑150 (MQ=255)
           aaTTCTGATGAGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATTTGTACTAATGtt                                                                                                                                        <  1:48408/150‑1 (MQ=255)
                 gatgaGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCattat                                                                                                                                 >  1:378991/1‑151 (MQ=255)
                  atgaGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCattatt                                                                                                                                <  2:104252/151‑1 (MQ=255)
                  atgaGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCattatt                                                                                                                                >  2:388213/1‑151 (MQ=255)
                   tgaGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCattatt                                                                                                                                >  1:20198/1‑150 (MQ=255)
                    gaGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCattattat                                                                                                                              <  1:1135371/151‑1 (MQ=255)
                        cGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTcc                                                                                                                           <  2:579203/150‑1 (MQ=255)
                        cGCTATGAATAACTTTACCGATAACTTCATGACCTAAAACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAAAGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTc                                                                                                                            >  2:1174837/1‑149 (MQ=255)
                            aTGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAAt                                                                                                                        <  2:617568/149‑1 (MQ=255)
                             tGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTCTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCt                                                                                                                      <  2:412479/150‑1 (MQ=255)
                              gAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTa                                                                                                                     <  1:775150/150‑1 (MQ=255)
                                 tAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGt                                                                                                                 >  2:1157288/1‑151 (MQ=255)
                                  aaCTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGt                                                                                                                 <  2:193412/150‑1 (MQ=255)
                                  aaCTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTc                                                                                                                >  1:1135127/1‑151 (MQ=255)
                                   aCTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCt                                                                                                               <  1:408291/151‑1 (MQ=255)
                                      ttACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTg                                                                                                              <  1:433258/149‑1 (MQ=255)
                                       tACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTc                                                                                                           >  1:1073551/1‑151 (MQ=255)
                                                cTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGc                                                                                                  >  1:346430/1‑151 (MQ=255)
                                                   cATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAAc                                                                                               >  2:1231741/1‑151 (MQ=255)
                                                    aTGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAAc                                                                                               >  1:265734/1‑150 (MQ=255)
                                                      gACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGt                                                                                            >  1:916533/1‑151 (MQ=255)
                                                       aCCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGtt                                                                                           >  2:542071/1‑151 (MQ=255)
                                                            acacCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCtt                                                                                      <  2:1250897/151‑1 (MQ=255)
                                                            acacCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCtt                                                                                      <  1:307217/151‑1 (MQ=255)
                                                             cacCATAGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCt                                                                                       >  1:1245824/1‑149 (MQ=255)
                                                              acCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGc                                                                                    >  2:80485/1‑151 (MQ=255)
                                                                       gCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGcaaca                                                                            <  2:20791/150‑1 (MQ=255)
                                                                       gCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATGCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGcaaca                                                                            <  1:303245/150‑1 (MQ=255)
                                                                           ttATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGc                                                                        <  2:1050851/150‑1 (MQ=255)
                                                                             aTCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAgg                                                                     >  2:420191/1‑151 (MQ=255)
                                                                             aTCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAg                                                                      <  2:593887/150‑1 (MQ=255)
                                                                              tCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAgg                                                                     >  1:501337/1‑150 (MQ=255)
                                                                                aTGAAATTACCTACTTTTCCTTCCTGTTAATAATGTAAATCGGAACCGCAAATTCCACCTCGGTTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACt                                                                  <  2:1036904/151‑1 (MQ=255)
                                                                                aTGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACt                                                                  >  2:614304/1‑151 (MQ=255)
                                                                                  gAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACtg                                                                 <  1:788467/150‑1 (MQ=255)
                                                                                  gAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACtg                                                                 <  1:966368/150‑1 (MQ=255)
                                                                                    aaTTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTgtgt                                                              <  2:848836/151‑1 (MQ=255)
                                                                                    aaTTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTgtgt                                                              <  2:994043/151‑1 (MQ=255)
                                                                                      ttACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGttt                                                            >  2:515046/1‑151 (MQ=255)
                                                                                        aCCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTc                                                          >  1:982030/1‑151 (MQ=255)
                                                                                         ccTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTc                                                          >  2:1057300/1‑150 (MQ=255)
                                                                                           tACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACt                                                       >  1:915633/1‑151 (MQ=255)
                                                                                              ttttCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGc                                                    >  2:1008917/1‑151 (MQ=255)
                                                                                              ttttCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGc                                                    >  2:456108/1‑151 (MQ=255)
                                                                                                ttccttccTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCAt                                                  >  2:1113738/1‑151 (MQ=255)
                                                                                                    ttcctGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAAt                                              >  1:524217/1‑151 (MQ=255)
                                                                                                     tcctGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAAtt                                             >  1:82273/1‑151 (MQ=255)
                                                                                                      cctGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAAtt                                             <  1:565717/150‑1 (MQ=255)
                                                                                                      cctGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTc                                            <  2:1152870/151‑1 (MQ=255)
                                                                                                            aataatGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAACTTCCTCACt                                      >  2:47201/1‑151 (MQ=255)
                                                                                                            aataatGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACt                                      <  1:523596/151‑1 (MQ=255)
                                                                                                            aataatGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTATACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACt                                      >  2:886322/1‑151 (MQ=255)
                                                                                                               aatGTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTgg                                    <  2:1263949/150‑1 (MQ=255)
                                                                                                                  gTAAATCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCa                                 <  1:1057945/150‑1 (MQ=255)
                                                                                                                    aaaTCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGt                              >  2:525748/1‑151 (MQ=255)
                                                                                                                    aaaTCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGt                              <  2:134452/151‑1 (MQ=255)
                                                                                                                    aaaTCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGt                              >  2:583374/1‑151 (MQ=255)
                                                                                                                      aTCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGt                              >  2:1084222/1‑149 (MQ=255)
                                                                                                                      aTCGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGt                              >  1:758978/1‑149 (MQ=255)
                                                                                                                        cGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGac                          >  2:94114/1‑151 (MQ=255)
                                                                                                                        cGGAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGa                           >  2:1284492/1‑150 (MQ=255)
                                                                                                                         ggAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGaca                         <  1:709016/151‑1 (MQ=255)
                                                                                                                          gAACCGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGacac                        >  2:814745/1‑151 (MQ=255)
                                                                                                                             ccGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTc                     >  1:487068/1‑151 (MQ=255)
                                                                                                                             ccGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTc                     <  1:411619/151‑1 (MQ=255)
                                                                                                                             ccGCAAATTCCACCGCATGTTATGTGTACTAATATTCCATTATTATTCCAATCTATCGTCTGCTCTGTAACGGCAGCAGTTTACTTGCCCGCAACAACGCAGGACTTTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGGCACCTc                     <  1:704510/151‑1 (MQ=255)
                                                                                                                              cGCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCg                    >  1:687628/1‑151 (MQ=255)
                                                                                                                               gCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCgg                   <  1:542173/151‑1 (MQ=255)
                                                                                                                               gCAAATTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCgg                   >  2:1150709/1‑151 (MQ=255)
                                                                                                                                  aaTTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGa                  >  2:916992/1‑149 (MQ=255)
                                                                                                                                  aaTTCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAg                <  2:679342/151‑1 (MQ=255)
                                                                                                                                    ttCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGc               >  1:1096256/1‑150 (MQ=255)
                                                                                                                                     tCCACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCa              <  2:803402/150‑1 (MQ=255)
                                                                                                                                      ccACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCAtt            <  2:348243/151‑1 (MQ=255)
                                                                                                                                       cACCTCGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCAttt           <  1:359690/151‑1 (MQ=255)
                                                                                                                                            cGGGTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGc       <  2:818076/150‑1 (MQ=255)
                                                                                                                                             gggTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCgg     <  1:1049077/151‑1 (MQ=255)
                                                                                                                                             gggATATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCCCGGAAGCATTTAAGc       >  2:798994/1‑149 (MQ=255)
                                                                                                                                               gTTATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGtt   <  1:56775/151‑1 (MQ=255)
                                                                                                                                                ttATGTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGttt  >  1:814492/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCCCTTCATGTAATTCTGATGAGTCGCTATGAATAACTTTACCGATAACTTCATGACCTAACACCATCGGTGCCTTTATCATGAAATTACCTACTTTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATTTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTT  >  NC_000913/4494160‑4494454

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑