Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,769,672 (A)8→7 pseudogene (1/111 nt) yibV → pseudogene, putative cryptic Rhs immunity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,769,6650A.98.4% 282.7 / ‑3.5 64pseudogene (322/327 nt)yibVpseudogene, putative cryptic Rhs immunity protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (34/29);  minor base T (0/1);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTAAAA  >  NC_000913/3769523‑3769802
                                                                                                                                              |                                                                                                                                         
aaTAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATa                                                                                                                                           <  2:163007/143‑1 (MQ=255)
          ggACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATaaaaaaatcat                                                                                                                                 >  2:293733/1‑139 (MQ=255)
           gACTTTAAGAATCATATACATGCTAATTATTAAATAAGGTAAAGATCATAAGACCATAACTGAAAAAACTAAATCCACTACAACAAACCGAAGATTGAATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCAtt                                                                                                                               <  1:392029/143‑1 (MQ=255)
              tttCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGat                                                                                                                            <  1:293733/143‑1 (MQ=255)
                     aTCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACATGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAAtt                                                                                                                     <  1:250864/143‑1 (MQ=255)
                      tCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTc                                                                                                                    >  2:699501/1‑143 (MQ=255)
                      tCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTc                                                                                                                    <  2:1023291/143‑1 (MQ=255)
                       ccTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCa                                                                                                                   >  2:116763/1‑143 (MQ=255)
                       acTATACATGCTATTTATTAAGTAAGAGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCa                                                                                                                   <  1:591811/142‑1 (MQ=255)
                        cTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAg                                                                                                                  >  1:372204/1‑143 (MQ=255)
                                gCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATACACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATAtt                                                                                                          >  2:99942/1‑143 (MQ=255)
                                 cTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATTAAAGATAAGAAT‑ACAAAAATCATTGATATAAATTCAGTTAATAttt                                                                                                         >  1:420323/1‑143 (MQ=255)
                                 cTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATAttt                                                                                                         >  1:417517/1‑143 (MQ=255)
                                   atttattAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGt                                                                                                       >  1:261272/1‑143 (MQ=255)
                                       attaAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGtg                                                                                                   <  2:964623/143‑1 (MQ=255)
                                        atagggaaaggAAAGATCATAAGAACATAGCAGAAAAATATAAAACAAAAACAAAAAAACACAGAGAGAATTGAAACAGAAAATAAAATAAAAGAAAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGtga                                                                                                  <  1:582146/134‑1 (MQ=255)
                                          aagtaagGGAAAGATCATAAGACCATATCTGAAAAATCTAAATGCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAGGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGtgatg                                                                                                >  1:784984/1‑143 (MQ=255)
                                          aagtaagGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAataataataggaatg                                                                                                >  1:461845/1‑131 (MQ=255)
                                           agtaagGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGtgatga                                                                                               >  1:746162/1‑143 (MQ=255)
                                                 ggAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATAGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATGAAGTTAATATTTGTAGTGATGATaaaaa                                                                                         >  1:738888/1‑143 (MQ=255)
                                                   aaaGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGCTTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATCAGAAAAata                                                                                       >  1:715353/1‑143 (MQ=255)
                                                   aaaGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAata                                                                                       >  1:275130/1‑143 (MQ=255)
                                                          aTAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAAtt                                                                                <  2:975722/143‑1 (MQ=255)
                                                                    atCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAAt                                                                      >  2:582241/1‑143 (MQ=255)
                                                                      cggataaaTCTAAAACCGGATCAGCTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAg                                                                    <  2:818435/138‑1 (MQ=255)
                                                                              gcTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACtatta                                                            <  2:268276/142‑1 (MQ=255)
                                                                                tAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATg                                                          >  2:476971/1‑143 (MQ=255)
                                                                                tAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTAAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTATTATTATg                                                          >  1:186487/1‑143 (MQ=255)
                                                                                  aaTCCACTACAACATACCGATGAATGAATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGc                                                        <  1:854414/143‑1 (MQ=255)
                                                                                    tCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCa                                                      >  1:292161/1‑143 (MQ=255)
                                                                                    tCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCa                                                      >  1:92532/1‑143 (MQ=255)
                                                                                         tCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCAATtgat                                                 >  1:866329/1‑143 (MQ=255)
                                                                                           cAACTGGCCCCCGATTGCACTGAAAGAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGTCCATtgatga                                               <  2:68698/143‑1 (MQ=255)
                                                                                            aaCTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTgatgat                                              >  2:747884/1‑143 (MQ=255)
                                                                                             gcTTGCCGCTGATTGCAGTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGAtt                                             <  2:356410/142‑1 (MQ=255)
                                                                                             aCTTGCCGCTGATTGGATTGAAACAGTAAATAAAATAAAAGAGAAGAAT‑AAAAAAATCATTGATATAAAGTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGAtt                                             <  2:227971/143‑1 (MQ=255)
                                                                                             aCTTGCCACTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGAtt                                             <  1:286543/143‑1 (MQ=255)
                                                                                               tggccatataTTGCAGTGAAACAGTAAAGAAAATAAAAGATATGAAT‑AAAAAAATGATTGATATAAATACAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCt                                           <  1:193914/134‑1 (MQ=255)
                                                                                                 gCCGCTGATTGCGTTGCCCCAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCtat                                         <  2:996399/143‑1 (MQ=255)
                                                                                                 gCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGATAATAGTACTATTATGGCCATTGATGATTCtat                                         <  2:261272/143‑1 (MQ=255)
                                                                                                  ccGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCtata                                        <  1:610707/143‑1 (MQ=255)
                                                                                                   cgatGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATaa                                       <  1:191829/140‑1 (MQ=255)
                                                                                                   cGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATaa                                       >  1:48922/1‑143 (MQ=255)
                                                                                                     cTGATTGGCTTGAAACAGGAAACAAAAGACAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAACTTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAg                                     <  1:892567/143‑1 (MQ=255)
                                                                                                     cTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAg                                     >  1:131263/1‑143 (MQ=255)
                                                                                                         ttGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCa                                 >  1:260005/1‑143 (MQ=255)
                                                                                                         ttGCATTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCa                                 >  1:948786/1‑143 (MQ=255)
                                                                                                             aTTGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATATGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGa                             >  1:208151/1‑143 (MQ=255)
                                                                                                             aTTGAAAAAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGa                             >  1:184388/1‑143 (MQ=255)
                                                                                                              ttaaaaCATAAAAAAAAATAAAAG‑TAAAAAATAAAAAAATCATTAATATAAATTCATTTAATATTTGTAGTGATTATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGat                            <  1:608214/140‑1 (MQ=255)
                                                                                                              ttGAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGat                            <  2:666677/143‑1 (MQ=255)
                                                                                                                gaaaaaagaaataaaataAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGatga                          <  1:822319/135‑1 (MQ=255)
                                                                                                                gAAACAGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGatga                          >  2:513543/1‑143 (MQ=255)
                                                                                                                     aGTAAATAAAATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGAt                     >  2:294365/1‑143 (MQ=255)
                                                                                                                       taaataaaATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGCTGCTTCTATAAAGTGCATCGATGAATGCTaa                   >  2:262977/1‑143 (MQ=255)
                                                                                                                       taaataaaATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATaa                   <  2:340488/143‑1 (MQ=255)
                                                                                                                       taaataaaATAAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATaa                   >  1:272921/1‑143 (MQ=255)
                                                                                                                        aaataaaataAAAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAAc                  <  1:726398/143‑1 (MQ=255)
                                                                                                                        aaagaaaagaaaagaTAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAAc                  <  1:690528/143‑1 (MQ=255)
                                                                                                                          ataaaataaaAGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACtt                >  1:434450/1‑143 (MQ=255)
                                                                                                                              aataaaaGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAAtt            >  2:1043159/1‑143 (MQ=255)
                                                                                                                                taaaaGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAAtttt          >  2:97595/1‑143 (MQ=255)
                                                                                                                                 aaaaGATAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATttta         <  1:116763/143‑1 (MQ=255)
                                                                                                                                       tAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTaaa   <  2:598676/143‑1 (MQ=255)
                                                                                                                                       tAAGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATCTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCAAAAAAGTGCATCGATGAATGATAACTTAATTTTATTTaaa   <  1:1039312/143‑1 (MQ=255)
                                                                                                                                        aaGAAT‑AAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTaaaa  >  2:753005/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                         
AATAGATGCTGGACTTTCCGAATCCTATACATGCTATTTATTAAGTAAGGGAAAGATCATAAGACCATATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTAAAA  >  NC_000913/3769523‑3769802

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑