Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,256 (C)8→9 pseudogene (3831/5422 nt) yneO ← AIDA‑I family autotransporter YneO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,594,2481.C98.3% 204.5 / ‑3.3 60pseudogene (3839/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (36/23);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

AGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAG  >  NC_000913/1594111‑1594377
                                                                                                                                            |                                                                                                                                 
aGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTc                                                                                                                                  >  2:968726/1‑140 (MQ=255)
  cGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                                                                                                <  1:968726/140‑1 (MQ=255)
     cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccc                                                                                                                             <  2:578910/140‑1 (MQ=255)
       cGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccc                                                                                                                           <  2:448960/140‑1 (MQ=255)
       cGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccc                                                                                                                           <  2:168581/140‑1 (MQ=255)
         cAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTcccccccccac                                                                                                                        <  2:440070/140‑4 (MQ=255)
               ccccccccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAACTTCCCCCCCCCACGatat                                                                                                                    <  1:707091/140‑1 (MQ=255)
                cccccccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATc                                                                                                                   <  1:13978/140‑1 (MQ=255)
                 ccccccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATcc                                                                                                                  >  2:568694/1‑140 (MQ=255)
                  cccccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCcg                                                                                                                 >  1:493444/1‑140 (MQ=255)
                  cccccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCcg                                                                                                                 <  1:568694/140‑1 (MQ=255)
                     ccAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCcgtcg                                                                                                              <  2:493444/140‑1 (MQ=255)
                           aTAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCt                                                                                                        >  1:24611/1‑140 (MQ=255)
                                         aCCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGaa                                                                                          >  1:466311/1‑140 (MQ=255)
                                          ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAAt                                                                                         >  2:525184/1‑140 (MQ=255)
                                          ccTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAAt                                                                                         >  2:78885/1‑140 (MQ=255)
                                           cTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATc                                                                                        >  1:958832/1‑140 (MQ=255)
                                           cTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATc                                                                                        >  1:621187/1‑140 (MQ=255)
                                               aGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg                                                                                    >  1:216953/1‑140 (MQ=255)
                                                  ggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                                                                 >  1:417747/1‑140 (MQ=255)
                                                  ggtgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                                                                 <  2:352022/140‑1 (MQ=255)
                                                    tgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAag                                                                               >  2:498941/1‑140 (MQ=255)
                                                    tgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAag                                                                               >  2:317747/1‑140 (MQ=255)
                                                    tgACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGCCGTACCCTGATACGCCGCGAAATCGTTGTTagg                                                                               >  1:617947/1‑138 (MQ=255)
                                                          ttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAg                                                                         <  2:887932/140‑1 (MQ=255)
                                                          ttCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAg                                                                         <  1:317747/140‑1 (MQ=255)
                                                                       acCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTAccc                                                            <  2:809527/140‑1 (MQ=255)
                                                                          gCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGt                                                         <  1:498941/140‑1 (MQ=255)
                                                                           cTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGtt                                                        >  2:702431/1‑140 (MQ=255)
                                                                              ttATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTg                                                     >  1:118175/1‑140 (MQ=255)
                                                                              ttATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTg                                                     >  1:494273/1‑140 (MQ=255)
                                                                                      atTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCa                                             >  1:935053/1‑140 (MQ=255)
                                                                                      atTAATGTGTTGACTTGCCATATTAATGGCAGAATCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCa                                             >  1:629172/1‑140 (MQ=255)
                                                                                        tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCacc                                           >  1:43237/1‑140 (MQ=255)
                                                                                        tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCacc                                           >  2:607629/1‑140 (MQ=255)
                                                                                        tAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTCAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCacc                                           >  1:912866/1‑140 (MQ=255)
                                                                                         aaTGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccg                                          <  1:607629/140‑1 (MQ=255)
                                                                                           tgtgTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgct                                        <  1:78885/140‑1 (MQ=255)
                                                                                           tgtgTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgct                                        <  2:466311/140‑1 (MQ=255)
                                                                                             tgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctac                                      >  2:714756/1‑140 (MQ=255)
                                                                                              gtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCaccgctacc                                     >  1:716345/1‑140 (MQ=255)
                                                                                                  aCTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTg                                 >  2:469316/1‑140 (MQ=255)
                                                                                                   cTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTCAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGa                                <  2:912866/140‑1 (MQ=255)
                                                                                                      gCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATAc                             >  2:481743/1‑140 (MQ=255)
                                                                                                      gCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATAc                             >  2:529869/1‑140 (MQ=255)
                                                                                                       ccATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATa                              <  1:719452/138‑1 (MQ=255)
                                                                                                       ccATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATa                              >  2:719452/1‑138 (MQ=255)
                                                                                                            ttAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCg                       >  2:609102/1‑140 (MQ=255)
                                                                                                               aTGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCa                    <  2:494273/140‑1 (MQ=255)
                                                                                                               aTGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCa                    <  2:216953/140‑1 (MQ=255)
                                                                                                                 ggCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATc                  >  1:719189/1‑140 (MQ=255)
                                                                                                                    agagTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGcc               <  2:417747/140‑1 (MQ=255)
                                                                                                                     gagTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGccc              >  1:702151/1‑140 (MQ=255)
                                                                                                                      agTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCa             >  1:332046/1‑140 (MQ=255)
                                                                                                                       gTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAg            >  2:111924/1‑140 (MQ=255)
                                                                                                                        tCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGc           >  2:707311/1‑140 (MQ=255)
                                                                                                                         cGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCt          <  2:169235/140‑1 (MQ=255)
                                                                                                                         cGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCt          <  2:958832/140‑1 (MQ=255)
                                                                                                                            aaCCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACg       <  1:609102/140‑1 (MQ=255)
                                                                                                                            aaCCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACg       <  2:719189/140‑1 (MQ=255)
                                                                                                                               cGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTc    <  1:714756/140‑1 (MQ=255)
                                                                                                                               cGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTc    <  2:43237/140‑1 (MQ=255)
                                                                                                                                gAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCa   >  1:496215/1‑140 (MQ=255)
                                                                                                                                 aaaGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAg  <  2:332046/140‑1 (MQ=255)
                                                                                                                                 aaaGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAg  >  2:21101/1‑140 (MQ=255)
                                                                                                                                 aaaGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAg  >  1:286915/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                 
AGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCACATCACCGCTACCGCTGATACTGTGCGTCATCGCCCAGCTACGCTCAG  >  NC_000913/1594111‑1594377

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