Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,750,144 A→T 5.0% H20L (CAT→CTT)  sgbE → L‑ribulose‑5‑phosphate 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,750,1440AT 5.0% 246.3 / 2.4 79H20L (CAT→CTT) sgbEL‑ribulose‑5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (29/46);  new base T (2/2);  total (31/48)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

TGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGG  >  NC_000913/3749998‑3750291
                                                                                                                                                  |                                                                                                                                                   
tgtgGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATc                                                                                                                                                   <  3:330005/149‑1 (MQ=255)
       cGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGAc                                                                                                                                            >  6:204770/1‑149 (MQ=255)
        gAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACg                                                                                                                                           <  1:86379/149‑1 (MQ=255)
             aGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCAc                                                                                                                                      <  4:191962/149‑1 (MQ=255)
              gCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCAcc                                                                                                                                     >  1:317185/1‑149 (MQ=255)
               ccAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGCTTGAAGCGCGTATGCAGGAGGCGGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCt                                                                                                                                    <  8:251988/149‑1 (MQ=255)
                cAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTg                                                                                                                                   <  1:5603/149‑1 (MQ=255)
                cAAAGAGCCGGTGCGGGAGAGTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTg                                                                                                                                   <  5:107552/149‑1 (MQ=255)
                 aaaGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTgg                                                                                                                                  <  6:97488/149‑1 (MQ=255)
                 aaaGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTgg                                                                                                                                  >  4:194063/1‑149 (MQ=255)
                  aaGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTggg                                                                                                                                 <  7:207894/149‑1 (MQ=255)
                  aaGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTggg                                                                                                                                 <  7:282298/149‑1 (MQ=255)
                  aaGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGAAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTggg                                                                                                                                 <  1:152033/149‑1 (MQ=255)
                   agagCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTgggg                                                                                                                                <  8:54390/149‑1 (MQ=255)
                       ccGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAAt                                                                                                                            <  8:109428/149‑1 (MQ=255)
                           tGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCa                                                                                                                        <  5:204770/149‑1 (MQ=255)
                           tGCTGGAGATTATTCAGCCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCa                                                                                                                        <  4:166198/149‑1 (MQ=255)
                             cTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAgc                                                                                                                      >  8:232660/1‑149 (MQ=255)
                                gagaTTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGcgcg                                                                                                                   >  5:167661/1‑149 (MQ=255)
                                gagaTTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGcgcg                                                                                                                   >  5:277995/1‑149 (MQ=255)
                                     tattCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCCCTCACCATCAGGTGACGTTCACATGGGGCAATGTCAGCGCGGTAGa                                                                                                              >  4:89887/1‑149 (MQ=255)
                                      attCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGa                                                                                                              <  5:86102/148‑1 (MQ=255)
                                          aGGCGCGGCGTTTGATTGAAGCGCGTTTGCTGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGaaa                                                                                                         <  3:89887/149‑1 (MQ=255)
                                            gcgcGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAAcg                                                                                                       >  7:312182/1‑149 (MQ=255)
                                              gcggcgTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAAcgcg                                                                                                     <  1:203624/149‑1 (MQ=255)
                                               cggcgTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCgg                                                                                                    <  6:277995/149‑1 (MQ=255)
                                                      ggATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAatgga                                                                                             >  2:264102/1‑149 (MQ=255)
                                                       gATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTggg                                                                                                                                 <  6:129935/112‑1 (MQ=255)
                                                       gATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTggg                                                                                                                                 >  5:129935/1‑112 (MQ=255)
                                                       gATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAatggat                                                                                            <  2:317185/149‑1 (MQ=255)
                                                       gATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGGTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAatggat                                                                                            <  4:299484/149‑1 (MQ=255)
                                                        aTTGAAGCTCGTATGCAGGAGGCTGTATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAatggatg                                                                                           <  6:251917/149‑1 (MQ=255)
                                                        aTTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAatggatg                                                                                           >  8:380474/1‑149 (MQ=255)
                                                        aTTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAatggatg                                                                                           <  1:125968/149‑1 (MQ=255)
                                                         ttGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGATACGCGGCAatggatgg                                                                                          <  8:57780/149‑1 (MQ=255)
                                                         ttGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAatggatgg                                                                                          <  6:348477/149‑1 (MQ=255)
                                                            aaGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTaa                                                                                       >  4:152858/1‑149 (MQ=255)
                                                                gcgTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCaa                                                                                   <  3:306883/149‑1 (MQ=255)
                                                                 cgTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCaaa                                                                                  <  2:197287/149‑1 (MQ=255)
                                                                      gCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCaa                                                                                                 >  8:51558/1‑129 (MQ=255)
                                                                        aggaggCTGGATTTATATGATAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCg                                                                           <  7:154552/149‑1 (MQ=255)
                                                                                ggATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAg                                                                   >  8:349032/1‑149 (MQ=255)
                                                                                 gATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAg                                                                   <  8:178358/148‑1 (MQ=255)
                                                                                      atatGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCCCCATCTGGTGACGTTCCCCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAAAGGATGGTAATCAAACCTTCAGGCGTCGAGTacgac                                                             >  6:221584/1‑149 (MQ=255)
                                                                                       tatGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTacgac                                                             <  2:308271/148‑1 (MQ=255)
                                                                                         tGTTAGATCCACTGAAAGCCGACGGGCTGGCGGCTAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACTCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGtg                                                          <  7:336779/149‑1 (MQ=255)
                                                                                         tGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGt                                                           >  4:277383/1‑148 (MQ=255)
                                                                                          gTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCTTCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGtga                                                         >  5:190119/1‑149 (MQ=255)
                                                                                          gTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGtga                                                         >  1:67886/1‑149 (MQ=255)
                                                                                          gTTAGAGCAACTGAAAGCCGACGAGCTGGCGGCGAATCTGTCGCTTCCCGCTCACCTACTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACATTCCGGCGTCGATTACGACGtga                                                         >  5:190122/1‑149 (MQ=255)
                                                                                           ttAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGtgat                                                        >  4:116157/1‑149 (MQ=255)
                                                                                           ttAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCCCTCACCATCTCGTGACCTTCACCTGGGGCAATGTCAGCCCGGTAGACGAAAAGCGGCAATGGATGGTAATCAAACATTCCGGCGTCGAGTACGACGtgat                                                        >  5:334650/1‑149 (MQ=255)
                                                                                              gagCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGtgatga                                                      <  2:67886/148‑1 (MQ=255)
                                                                                                  aaCTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGAccgcc                                                 >  3:318864/1‑149 (MQ=255)
                                                                                                       aaaGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCcgacga                                            >  8:40938/1‑149 (MQ=255)
                                                                                                       aaaGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCcgac                                              <  6:167661/147‑1 (MQ=255)
                                                                                                           ccGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGCAACGCGGCAATGGATGGTAAGCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATAtg                                        <  8:116674/149‑1 (MQ=255)
                                                                                                                gTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATAtggtggt                                   <  7:349032/149‑1 (MQ=255)
                                                                                                                    tGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGtt                               >  2:172878/1‑149 (MQ=255)
                                                                                                                      gcggcgAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTg                              <  8:43287/148‑1 (MQ=255)
                                                                                                                       cggcgAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTgag                            <  4:318864/149‑1 (MQ=255)
                                                                                                                         gcgAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTgaga                           >  5:123627/1‑148 (MQ=255)
                                                                                                                         gcgAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGAt                          >  3:140041/1‑149 (MQ=255)
                                                                                                                          cgAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATa                         <  4:140041/149‑1 (MQ=255)
                                                                                                                          cgAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATa                         <  3:194063/149‑1 (MQ=255)
                                                                                                                          cgAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATa                         <  6:123627/149‑1 (MQ=255)
                                                                                                                                tGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGc                       >  5:186725/1‑145 (MQ=255)
                                                                                                                                tGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGc                       <  6:186725/145‑1 (MQ=255)
                                                                                                                                 ggCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCg                  >  6:87025/1‑149 (MQ=255)
                                                                                                                                  gcgcTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAg                    <  1:172878/146‑1 (MQ=255)
                                                                                                                                  gcgcTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCgg                 <  1:264102/149‑1 (MQ=255)
                                                                                                                                  gcgcTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCg                  <  7:232660/148‑1 (MQ=255)
                                                                                                                                     cTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGGCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGt                >  5:356416/1‑147 (MQ=255)
                                                                                                                                        cccGCTCACGTTCTGGTTACGTTCACCTGGGGAAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCGTGCGGCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGGCAGCGGTAAg             <  6:190122/147‑1 (MQ=255)
                                                                                                                                        cccGCTCACCTTCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAg             <  6:190119/147‑1 (MQ=255)
                                                                                                                                        cccGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGTATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAgg            <  8:312182/148‑1 (MQ=255)
                                                                                                                                           gCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAggtggt        >  2:47516/1‑149 (MQ=255)
                                                                                                                                               aCCATCTGGTGACGTGCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAAGGGATGGTAATCACACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGGGGTTGAGCTAGCCAGGGGGAAGGtggtgg      >  2:238161/1‑147 (MQ=255)
                                                                                                                                                 cATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAgg  <  2:372812/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGTGGACCGAAAAAGCCAAAGAGCCGGTGCTGGAGATTATTCAGGCGCGGCGTTGGATTGAAGCGCGTATGCAGGAGGCTGGATTTATATGTTAGAGCAACTGAAAGCCGACGTGCTGGCGGCGAATCTGGCGCTTCCCGCTCACCATCTGGTGACGTTCACCTGGGGCAATGTCAGCGCGGTAGACGAAACGCGGCAATGGATGGTAATCAAACCTTCCGGCGTCGAGTACGACGTGATGACCGCCGACGATATGGTGGTGGTTGAGATAGCCAGCGGTAAGGTGGTGGAAGG  >  NC_000913/3749998‑3750291

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