Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,434,956 C→T 6.3% L227L (TTG→TTA)  truA ← tRNA pseudouridine(38‑40) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,9560CT 6.3% 206.7 / 3.3 63L227L (TTG→TTA) truAtRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base C (18/41);  new base T (2/2);  total (20/43)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

GCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAA  >  NC_000913/2434812‑2435095
                                                                                                                                                |                                                                                                                                           
gCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAAct                                                                                                                                         >  3:327969/1‑149 (MQ=255)
   ttattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTg                                                                                                                                      <  1:481435/149‑1 (MQ=255)
         tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGTAACTCTGCTATcc                                                                                                                                >  1:42542/1‑149 (MQ=255)
         tcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGTAACTCTGCTATcc                                                                                                                                >  1:217747/1‑149 (MQ=255)
          cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCa                                                                                                                               <  4:285276/149‑1 (MQ=255)
                     cAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCGGCTCTCCGGCt                                                                                                                    <  4:257518/149‑1 (MQ=255)
                         actAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGtt                                                                                                                <  3:462131/147‑1 (MQ=255)
                           tAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgt                                                                                                              <  4:18599/149‑1 (MQ=255)
                           tAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgt                                                                                                              <  1:129488/149‑1 (MQ=255)
                           tAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgt                                                                                                              <  2:57543/149‑1 (MQ=255)
                            aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTtgtg                                                                                                             <  2:86756/149‑1 (MQ=255)
                             gCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                                                                                                            <  1:291323/149‑1 (MQ=255)
                             gCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                                                                                                            <  4:276963/149‑1 (MQ=255)
                               gggCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGc                                                                                                          <  3:72222/149‑1 (MQ=255)
                                          cGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCa                                                                                               >  2:359313/1‑149 (MQ=255)
                                               ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg                                                                                          <  3:264256/149‑1 (MQ=255)
                                                  gAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTg                                                                                       >  1:164283/1‑149 (MQ=255)
                                                     gATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCg                                                                                    >  4:6897/1‑149 (MQ=255)
                                                      aTCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGa                                                                                   <  4:87062/149‑1 (MQ=255)
                                                        cATACCGGTCAGGGCAATCCAACGCTACCAGCTACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTCCGGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCACTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                                                                                 <  3:350103/149‑1 (MQ=255)
                                                                tCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCTTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCAGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAAGATTCCt                                                                         <  1:187632/149‑1 (MQ=255)
                                                                tCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt                                                                         <  3:341686/149‑1 (MQ=255)
                                                                tCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt                                                                         >  4:443230/1‑149 (MQ=255)
                                                                        aTCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt                                                                         >  4:81572/1‑141 (MQ=255)
                                                                        aTCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt                                                                         <  3:81572/141‑1 (MQ=255)
                                                                        aTCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCatat                                                                 <  1:494509/149‑1 (MQ=255)
                                                                         tCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATatg                                                                >  4:185561/1‑149 (MQ=255)
                                                                         tCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATatg                                                                >  2:398626/1‑149 (MQ=255)
                                                                          ccACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATatga                                                               >  1:159856/1‑149 (MQ=255)
                                                                           cACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATatgat                                                              >  1:513380/1‑149 (MQ=255)
                                                                              cgcgACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTa                                                           <  4:165800/149‑1 (MQ=255)
                                                                                   ccAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACaaaa                                                      <  2:312684/149‑1 (MQ=255)
                                                                                   ccAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACaaaa                                                      <  3:6897/149‑1 (MQ=255)
                                                                                   ccAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACaaaa                                                      <  1:456651/149‑1 (MQ=255)
                                                                                   ccAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACaaaa                                                      <  1:344033/149‑1 (MQ=255)
                                                                                   ccAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACaaaa                                                      <  1:21391/149‑1 (MQ=255)
                                                                                   ccAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACAATATGATGTACaaaa                                                      <  3:159560/149‑1 (MQ=255)
                                                                                   ccAGATACAGCCCTGCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAACAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACaaaa                                                      <  1:417095/149‑1 (MQ=255)
                                                                                    cAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAg                                                     <  3:14210/149‑1 (MQ=255)
                                                                                        tACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCAtt                                                 >  2:562080/1‑149 (MQ=255)
                                                                                           aGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGc                                              >  2:266379/1‑149 (MQ=255)
                                                                                           aGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGc                                              >  2:266776/1‑149 (MQ=255)
                                                                                                ttCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTaa                                         <  4:327969/149‑1 (MQ=255)
                                                                                                 tCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGTAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAat                                        <  2:217747/149‑1 (MQ=255)
                                                                                                 tCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGTAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAat                                        <  2:42542/149‑1 (MQ=255)
                                                                                                    gCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGc                                              >  1:1297/1‑140 (MQ=255)
                                                                                                    gCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGc                                              >  1:361649/1‑140 (MQ=255)
                                                                                                    gCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGc                                              <  2:1297/140‑1 (MQ=255)
                                                                                                    gCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGc                                              <  2:361649/140‑1 (MQ=255)
                                                                                                    gCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATc                                     >  2:530424/1‑149 (MQ=255)
                                                                                                     cTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCt                                    <  3:423141/149‑1 (MQ=255)
                                                                                                       tttGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTac                                  <  4:281922/149‑1 (MQ=255)
                                                                                                       tttGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTac                                  <  1:359313/149‑1 (MQ=255)
                                                                                                       tttGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTac                                  <  1:258032/149‑1 (MQ=255)
                                                                                                             gTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTGGCTTTAATATCTACcaccac                            >  3:465555/1‑149 (MQ=255)
                                                                                                                  tgcCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAg                       >  3:137800/1‑149 (MQ=255)
                                                                                                                    ccgccAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGa                     <  2:349301/149‑1 (MQ=255)
                                                                                                                     cgccAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGAc                    <  4:139018/149‑1 (MQ=255)
                                                                                                                        cAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGt                 <  4:435610/149‑1 (MQ=255)
                                                                                                                          gCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGa               <  3:522581/149‑1 (MQ=255)
                                                                                                                               cTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGt          <  1:37764/149‑1 (MQ=255)
                                                                                                                                      ttGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTa   <  3:204741/149‑1 (MQ=255)
                                                                                                                                       tGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTaa  <  1:398626/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
GCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAA  >  NC_000913/2434812‑2435095

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