Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,113,378 A→T 14.0% L34I (TTA→ATA) ‡ yghG ← secretin (GspDbeta) OM localization lipoprotein pilotin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,113,3780AT14.0% 199.2 / 22.4 71L34I (TTA→ATA) ‡yghGsecretin (GspDbeta) OM localization lipoprotein pilotin
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base T (7/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTT  >  NC_000913/3113231‑3113523
                                                                                                                                                   |                                                                                                                                                 
caggcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTaaa                                                                                                                                                  <  3:372040/149‑1 (MQ=255)
   gcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTg                                                                                                                                               <  2:376909/149‑1 (MQ=255)
   gcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTg                                                                                                                                               <  4:460895/149‑1 (MQ=255)
   gcgtcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTg                                                                                                                                               <  1:325114/149‑1 (MQ=255)
      tcTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCa                                                                                                                                            <  1:352678/149‑1 (MQ=255)
                     tACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctg                                                                                                                             >  3:529774/1‑149 (MQ=255)
                      aCCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAt                                                                                                                                           <  4:312410/134‑1 (MQ=255)
                      aCCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCAt                                                                                                                                           >  3:312410/1‑134 (MQ=255)
                      aCCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTAtggctgg                                                                                                                            >  1:509595/1‑149 (MQ=255)
                           cTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAg                                                                                                                       <  1:449243/149‑1 (MQ=255)
                            ttCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGTCAGTAAACTGGCATTTTCATTATGGCTGGCACAGc                                                                                                                      >  4:90884/1‑149 (MQ=255)
                            ttCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGc                                                                                                                      <  3:414272/149‑1 (MQ=255)
                                 tGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACtt                                                                                                                 >  1:227688/1‑149 (MQ=255)
                                 tGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACtt                                                                                                                 >  3:96902/1‑149 (MQ=255)
                                   aTAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaa                                                                                                               >  1:444602/1‑149 (MQ=255)
                                        ggTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGTTTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATc                                                                                                          <  4:124880/149‑1 (MQ=255)
                                         gTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATcc                                                                                                         >  3:191532/1‑149 (MQ=255)
                                          tCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaa                                                                                                               >  1:117240/1‑142 (MQ=255)
                                          tCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTtaa                                                                                                               <  2:117240/142‑1 (MQ=255)
                                             tttttACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTa                                                                                                     >  3:419893/1‑149 (MQ=255)
                                                     gTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTcaa                                                                                             <  1:490897/149‑1 (MQ=255)
                                                     gTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTcaa                                                                                             <  2:56419/149‑1 (MQ=255)
                                                     gTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTcaa                                                                                             <  3:19727/149‑1 (MQ=255)
                                                          gCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATa                                                                                        >  4:298775/1‑149 (MQ=255)
                                                          gCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATa                                                                                        <  1:497689/149‑1 (MQ=255)
                                                           ccGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAg                                                                                       >  4:251351/1‑149 (MQ=255)
                                                                 cTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGatat                                                                                 <  2:429222/149‑1 (MQ=255)
                                                                  tttGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATAtt                                                                                >  3:72155/1‑149 (MQ=255)
                                                                    tGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTaa                                                                              <  4:529774/149‑1 (MQ=255)
                                                                             aCTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAccctccc                                                                     >  1:356454/1‑149 (MQ=255)
                                                                               tAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTg                                                                   >  4:583855/1‑149 (MQ=255)
                                                                                          ccGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtttt                                                        <  2:509595/149‑1 (MQ=255)
                                                                                          ccGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtttt                                                        >  3:566594/1‑149 (MQ=255)
                                                                                          ccGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGtttt                                                        >  2:440679/1‑149 (MQ=255)
                                                                                           cGCAGATTTACTCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGGACGGCCACTTAATAATCCTGTAACGCTCCACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTa                                                       >  4:255941/1‑149 (MQ=255)
                                                                                                  ttAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACAt                                                <  4:72155/149‑1 (MQ=255)
                                                                                                    aaTCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATa                                                                            <  1:236702/119‑1 (MQ=255)
                                                                                                    aaTCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATa                                                                            >  2:236702/1‑119 (MQ=255)
                                                                                                      tCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCCCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTa                                            <  2:418914/149‑1 (MQ=255)
                                                                                                        ggCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTaaa                                          >  3:17164/1‑149 (MQ=255)
                                                                                                         gCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAaat                                         >  4:250120/1‑149 (MQ=255)
                                                                                                                   gCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCCACAATAGCGATATTAATAcc                                                                          >  2:479738/1‑106 (MQ=255)
                                                                                                                   gCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATAcc                                                                          <  1:479738/106‑1 (MQ=255)
                                                                                                                   gCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAg                               <  2:356454/149‑1 (MQ=255)
                                                                                                                     tGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACt                                                                                                                  >  2:483777/1‑64 (MQ=255)
                                                                                                                     tGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACt                                                                                                                  <  1:483777/64‑1 (MQ=255)
                                                                                                                     tGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGt                             >  4:583689/1‑149 (MQ=255)
                                                                                                                     tGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGt                             >  2:91373/1‑149 (MQ=255)
                                                                                                                      gATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGtt                            >  1:353518/1‑149 (MQ=255)
                                                                                                                        tatTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCAGTTGTCTTATCGACATGGTAAATAATATCTAAGGttat                          >  3:380570/1‑149 (MQ=255)
                                                                                                                          ttttGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCt                                                                                                        >  1:226697/1‑69 (MQ=255)
                                                                                                                          ttttGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCt                                                                                                        <  2:226697/69‑1 (MQ=255)
                                                                                                                                 ccTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAagag                 <  2:227688/149‑1 (MQ=255)
                                                                                                                                  cTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGt                >  1:331445/1‑149 (MQ=255)
                                                                                                                                    gTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGt                                             <  4:370237/118‑1 (MQ=255)
                                                                                                                                    gTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGt                                             >  3:370237/1‑118 (MQ=255)
                                                                                                                                    gTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAACAGTTa              >  4:506135/1‑149 (MQ=255)
                                                                                                                                       ttttCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTaaaa           <  4:96902/149‑1 (MQ=255)
                                                                                                                                       ttttCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAACAAGAGTTaaaa           <  3:506135/149‑1 (MQ=255)
                                                                                                                                         aatgCCAGTTTACTGGCATTTTCATTATTGCTGACACAGCCACTTAATAATCCTGTAACACTCAACAATAGCGATATTAATACCCTCCCTGACATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATg         >  4:586530/4‑149 (MQ=255)
                                                                                                                                         aatgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATg         >  3:161618/4‑149 (MQ=255)
                                                                                                                                         aatgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATg         >  2:3103/4‑149 (MQ=255)
                                                                                                                                         aatgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATg         >  1:154340/4‑149 (MQ=255)
                                                                                                                                         aatgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGAGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATg         >  1:370667/4‑149 (MQ=255)
                                                                                                                                          atgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTaaaa           <  4:387503/144‑1 (MQ=255)
                                                                                                                                          atgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTaaaa           >  3:387503/3‑146 (MQ=255)
                                                                                                                                          atgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTaaaa           <  2:247669/144‑1 (MQ=255)
                                                                                                                                          atgCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTaaaa           >  1:247669/3‑146 (MQ=255)
                                                                                                                                            gCCAGTTTACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCa      <  4:161618/149‑1 (MQ=255)
                                                                                                                                              cAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACt    <  4:419893/149‑1 (MQ=255)
                                                                                                                                              cAGTAAACTGGCATTTGCATTATGGCGGGCACAGCCACATAATAATCCTGTAACGCTCAACAATAGTGATATTAATACCCTCCCTGGCATTTGTTTTAGCGACATGGTGAATAATCTCTAAGGTTATTAATAAGAGTTAACATGTCACt    >  3:363697/1‑149 (MQ=255)
                                                                                                                                                gTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACttt  <  1:91373/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAGGCGTCTGCGTGGTTTGTGTACCCGCTTCGCTGATAATGGTCATTTTTACCGTTGTGCCGCTACTTTGCGCCAGCACTAAGGTATATCCCGCAGATTTAATCGGCAGGTTTTGGCTGATATTTTGCGCCTGTTTTTTCGCCAGTAAACTGGCATTTTCATTATGGCTGGCACAGCCACTTAATAATCCTGTAACGCTCAACAATAGCGATATTAATACCCTCCCTGGCATTTGTTTTATCGACATGGTAAATAATCTCTAAGGTTATTAATAAGAGTTAAAATGTCACTTT  >  NC_000913/3113231‑3113523

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