Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 384,328 C→A 5.3% D97Y (GAT→TAT)  yaiS ← putative PIG‑L family deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913384,3280CA 5.3% 285.5 / 3.8 85D97Y (GAT→TAT) yaiSputative PIG‑L family deacetylase
Reads supporting (aligned to +/- strand):  ref base C (39/41);  new base A (3/2);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

ATAGTTCGGCAGGCAACCATTGAAGCTTGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACG  >  NC_000913/384183‑384470
                                                                                                                                                 |                                                                                                                                              
aTAGTTCGGCAGGCAACCATTGAAGCTTGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAatcatatcat                                                                                                                                              <  5:154637/148‑1 (MQ=255)
    ttCGGCAGGCAACCATTGAAGCTTGATAAACAGCAAGATGCTCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAg                                                                                                                                          <  5:421996/148‑1 (MQ=255)
          aggcaACCATTGAAGCTTGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTaa                                                                                                                                    <  7:181016/148‑1 (MQ=255)
              aaCCATTGAAGCTTGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGa                                                                                                                                <  3:245939/148‑1 (MQ=255)
               aCCATTGAAGCTTGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTaaa                                                                                                                                   <  6:265264/144‑1 (MQ=255)
               aCCATTGAAGCTTGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTaaa                                                                                                                                   >  5:265264/1‑144 (MQ=255)
               aCCATTGAAGCTTGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAg                                                                                                                               >  8:365868/1‑148 (MQ=255)
               aCCATTCAAGCTTGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATAAACACGCATGATTTCAAAATCAGAAGGAATTACATTATAAATGATATATTCCAGGGCGGAAATCAAAAAATTGACCTGTAAATGAg                                                                                                                               >  2:80984/1‑148 (MQ=255)
                    tGAAGCTTGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCgg                                                                                                                          >  7:348806/1‑148 (MQ=255)
                       aGCTTGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGgtg                                                                                                                       <  5:208954/148‑1 (MQ=255)
                           tGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAg                                                                                                                   >  6:185114/1‑148 (MQ=255)
                            gATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAg                                                                                                                   >  2:321113/1‑147 (MQ=255)
                             aTAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCggg                                                                                                                         <  5:29885/140‑1 (MQ=255)
                             aTAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCggg                                                                                                                         >  6:29885/1‑140 (MQ=255)
                               gtttatcAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCgg                                                                                                                          <  8:227180/131‑1 (MQ=255)
                               gtttatcAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCgg                                                                                                                          >  7:227180/7‑137 (MQ=255)
                                     cAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCTGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAg                                                                                                         <  1:389328/148‑1 (MQ=255)
                                        gatgatCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGAt                                                                                                       >  5:256962/1‑147 (MQ=255)
                                          tgatCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCggg                                                                                                                         >  6:310855/1‑127 (MQ=255)
                                          tgatCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCggg                                                                                                                         <  5:310855/127‑1 (MQ=255)
                                           gatCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGaa                                                                                                    >  2:155018/1‑147 (MQ=255)
                                             tCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAg                                                                                                 >  5:245854/1‑148 (MQ=255)
                                              ccTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAATTCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGt                                                                                                >  2:265633/1‑148 (MQ=255)
                                              ccTGATGGCGGTCGGCATCATGCATGGTATATACACGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATTATGTCTTCCAGGGCGGAAATCATATCATTGAGCTTTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGt                                                                                                >  8:374433/1‑148 (MQ=255)
                                                  aTGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGttt                                                                                              <  6:181380/146‑1 (MQ=255)
                                                  aTGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTtg                                                                                             <  6:256962/147‑1 (MQ=255)
                                                   tGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTtggt                                                                                           <  5:80042/148‑1 (MQ=255)
                                                     gCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTtggtgg                                                                                         <  4:300096/148‑1 (MQ=255)
                                                     gCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTtggtgg                                                                                         <  1:321113/148‑1 (MQ=255)
                                                                catgcatgGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATGGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAatat                                                                              >  1:172513/1‑148 (MQ=255)
                                                                catgcatgGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATAATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGc                                                                                        >  8:37964/1‑138 (MQ=255)
                                                                catgcatgGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATAATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGc                                                                                        <  7:31238/138‑1 (MQ=255)
                                                                catgcatgGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATAATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGc                                                                                        <  7:37964/138‑1 (MQ=255)
                                                                catgcatgGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATAATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGc                                                                                        >  8:31238/1‑138 (MQ=255)
                                                                  tgcatgGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCt                                                                            <  3:183340/148‑1 (MQ=255)
                                                                  tgcatgGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATAACATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCt                                                                            >  3:111474/1‑148 (MQ=255)
                                                                    catgGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGGGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCttt                                                                          <  7:232001/148‑1 (MQ=255)
                                                                     atgGTATATACCCTCATGATTTCAACATAAGAAGGAATTTGATTTTTAATGATGTCTTCCATGGCGGAAATCATATCATTGAGCTGTAAATGAGCTCTGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTa                                                                         <  7:374433/148‑1 (MQ=255)
                                                                      tgGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAATTTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAg                                                                        >  1:358878/1‑148 (MQ=255)
                                                                       ggTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAgg                                                                       >  4:202684/1‑148 (MQ=255)
                                                                        gTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAggg                                                                      >  5:352971/1‑148 (MQ=255)
                                                                           tatcCCCGCCTTATTTCAACATCAGAAGGAATTTGATTTTAAATGATGTCTTCCCGGGCGGAAATCATAGCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGt                                                                   <  8:433693/148‑1 (MQ=255)
                                                                           tataCCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGt                                                                   <  2:172513/148‑1 (MQ=255)
                                                                              aCCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTgc                                                                >  2:287194/1‑148 (MQ=255)
                                                                                ccGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTaa                                                                                                                                    <  5:56229/78‑1 (MQ=255)
                                                                                ccGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTaa                                                                                                                                    >  6:56229/1‑78 (MQ=255)
                                                                                ccGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTgcgc                                                              <  6:245854/148‑1 (MQ=255)
                                                                                 cGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTgcgc                                                              >  1:87881/1‑147 (MQ=255)
                                                                                   cATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGTCGGAAATAATATCATTTAGCTGTAAATGATCGCGGGTGTCAGCAAAATTAAGATTAATAGTTTGTTGGCACCCTAATATCTTTAGTGCGTTGCTCGAt                                                           >  4:62121/1‑148 (MQ=255)
                                                                                      gATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTaa                                                                                                                                    >  6:275353/1‑72 (MQ=255)
                                                                                      gATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTaa                                                                                                                                    <  5:275353/72‑1 (MQ=255)
                                                                                        tttCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTg                                                                 >  8:47871/1‑137 (MQ=255)
                                                                                        tttCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTg                                                                 <  7:47871/137‑1 (MQ=255)
                                                                                        tttCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGAttcttc                                                      <  7:133642/148‑1 (MQ=255)
                                                                                           cAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCAt                                                    <  6:352971/147‑1 (MQ=255)
                                                                                            aaCATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATgg                                                  <  1:287194/148‑1 (MQ=255)
                                                                                             acaacaGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGGTTGGTGGCACCCTAATATCTTTTGGGCGTTGCGCGATTCTTCATGGc                                                 >  2:61331/1‑148 (MQ=255)
                                                                                                tCAGAAGGAATTTGATTTTTAATGCTGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGTGCGCGGGTGTCAGCAAAATTTAGa                                                                                                        <  5:210229/90‑1 (MQ=255)
                                                                                                 cAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGcgatc                                             >  4:315876/1‑148 (MQ=255)
                                                                                                 cAGAAGGAATTTGATTTTTAATGAAGTCTTCCAGGGCGGAAATCATAACATTGAGCTGTAAATGAGCGCGGGAGTCAGCAAAATTAAGAt                                                                                                       >  6:210229/1‑90 (MQ=255)
                                                                                                  aGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGcgatcg                                            <  1:265633/148‑1 (MQ=255)
                                                                                                   gAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGcgatcga                                           >  3:125360/1‑148 (MQ=255)
                                                                                                      ggAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGAttcttc                                                      <  2:90386/134‑1 (MQ=255)
                                                                                                      ggAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGAttcttc                                                      >  1:90386/1‑134 (MQ=255)
                                                                                                        aaTTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAAGATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGattatt                                      >  1:211890/1‑148 (MQ=255)
                                                                                                            tGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCa                                   <  4:125360/147‑1 (MQ=255)
                                                                                                            tGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCAt                                  >  7:327016/1‑148 (MQ=255)
                                                                                                              aTTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATc                                 <  5:185114/147‑1 (MQ=255)
                                                                                                                ttttAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAg                                                                                                         >  3:366860/1‑73 (MQ=255)
                                                                                                                ttttAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAg                                                                                                         <  4:366860/73‑1 (MQ=255)
                                                                                                                 tttAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGGTGTAAATGAGAGCGTGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTAGCGAGATTCTTCATGGCGATCGATTATTCCATCtgt                              >  1:349094/1‑147 (MQ=255)
                                                                                                                   tAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCATCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGc                            <  4:111474/147‑1 (MQ=255)
                                                                                                                    aaTGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGc                            >  2:154609/1‑146 (MQ=255)
                                                                                                                     atgatgTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCag                         <  2:87881/148‑1 (MQ=255)
                                                                                                                     atgatgTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCag                         <  6:270476/148‑1 (MQ=255)
                                                                                                                      tgatgTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCaga                        <  7:365868/148‑1 (MQ=255)
                                                                                                                          gTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGatt                                         <  7:191964/127‑1 (MQ=255)
                                                                                                                          gTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGatt                                         >  8:191964/1‑127 (MQ=255)
                                                                                                                          gTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCtgtg                             >  5:369803/1‑139 (MQ=255)
                                                                                                                          gTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCtgtg                             <  6:369803/139‑1 (MQ=255)
                                                                                                                                cAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGt               <  1:155018/147‑1 (MQ=255)
                                                                                                                                  gggCGTAAATGATATCATTGAGCTGTAAATGAGCGCGTGTGTCAGCAAAATTTAGATGAATAGTTTGGTGTCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGATATCAATTAATCCATCTGTGCCTGAGTTGCCGGTAGt            <  1:80984/148‑1 (MQ=255)
                                                                                                                                        aaaTCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTa                                                                         >  3:20835/1‑81 (MQ=255)
                                                                                                                                        aaaTCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTa                                                                         <  4:20835/81‑1 (MQ=255)
                                                                                                                                        aaaTCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATcac      >  7:347280/1‑148 (MQ=255)
                                                                                                                                        aaaTCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATcac      >  5:411609/1‑148 (MQ=255)
                                                                                                                                             atatcatTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACg  <  6:411609/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ATAGTTCGGCAGGCAACCATTGAAGCTTGATAAACAGCAAGATGATCCTGATGGCGGTCGGCATCATGCATGGTATATACCCGCATGATTTCAACATCAGAAGGAATTTGATTTTTAATGATGTCTTCCAGGGCGGAAATCATATCATTGAGCTGTAAATGAGCGCGGGTGTCAGCAAAATTAAGATGAATAGTTTGGTGGCACCCTAATATCTTTAGGGCGTTGCGCGATTCTTCATGGCGATCGATTATTCCATCTGTGCCAGAGTTGCCGGTAGTCATCACCACG  >  NC_000913/384183‑384470

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